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ABSTRACT 

Fin  experimental  apparatus  to  measure  the  particle  density 
distribution  in  the  freeboard  of  an  atmospheric  fluidized  bed  uas 
designed  and  constructed.   The  density  versus  height  measured  by 
the  sampling  apparatus  gives  a  similar  exponential  decrease  as 
previous  investigations  have  found. 

A  particle  trajectory  model  is  developed  uhich  calculates  the 
height  and  particle  density  distributions  above  the  bed  surface  of 
an  atmospheric  fluidized  bed.   The  parameters  input  to  the  model 
are  the  superficial  velocity,  initial  partical  velocity,  gas  jet 
velocity  and  duration,  and  the  particle  size  distribution  of  the 
bed  mass.   The  model  uas  evaluated  using  the  experimental  data  for 
jet  velocity,  duration,  and  particle  size.   The  predicted  slope  of 
the  particle  density  versus  height  in  the  freeboard  agrees  uith 
the  experimentally  measured  slope  uithin  Z0%. 

A  sensitvity  analysis  using  the  trajectory  model,  resulted  in 
a  determination  of  the  particle  distributions  in  the  freeboard  of 
a  fluidized  bed  as  affected  by  varying  the  input  parameters  to  the 
trajectory  model.   The  most  significant  effects  uere  achieved  when 
the  jet  velocity  or  duration  uas  altered. 
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CHAPTER   I 


INTRODUCTION 

Fluidized  beds  have  been  used  in  industry  for  many  years. 
They  have   been  used  to  nix  and  dry  particulate  materials  and  are 
the  principle  process  in  catalytic  cracking  plants.   In  the  past 
decade  or  so,  the  use  of  fluidized  bed  combustors  for  power 
generation  has  become  a  source  of  major  interest.   Prototype  coal 
burning  beds  have  already  been  built  uhich  are  comparable  to 
existing  coal  plants.   Fluidized  bed  combustors  have  the  added 
benifit  of  lou  NOx ,  SO2  and  hydrocarbon  emissions  and  the 
flexibility  of  being  able  to  burn  a  uide  range  of  fuels  ranging 
from  refuse  and  high  sulfur  content  coal  to  high  grade  fuels. 

A  fluidized  bed  [Fig  1]  is  composed  of  a  distributor  through 

uhich  an  air  flou  is  introduced  through  thousands  of  small 

orfices.   This  air  then  passes  through  the  dense  zone  of  the  bed 

uhich  is  comprised  of  a  mass  of  particles.   The  air  velocity 

through  the  dense  zone  is  maintained  above  the  minimum 

fluidization  velocity  ( Umf )  during  normal  operations.   fit 

velocities  equal  to  or  greater  than  Umf,  the  frictional  force 

(Drag)  of  the  air  flouing  past  a  particle  is  equal  to  the  weight 

of  the  particle.   Under  these  conditions,  the  particle  mass 

behaves  very  much  like  a  fluid.   It  will  maintain  a  horizontal 

surface  if  the  container  is  tilted,  flou  out  of  holes  in  the 
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container,  and  has  a  pressure  drop  across  any  section  of  the  bed 
apprxamately  equal  to  the  weight  of  air  and  particles  in  the 
section  [11.   At  velocities  above  Umf ,  the  excess  air  will 
coalesce  and  forn  snail  voids  or  bubbles  as  it  procedes  towards 
the  surface  of  the  fluidized  bed. 

As  bubbles  break  at  the  surface  of  the  bed,  the  solid 
particles  are  thrown  up  above  the  bed  surface  and  are  entrained  by 
the  upward  flowing  gas  stream.   This  zone  above  the  bed  surface  is 
the  freeboard  zone.   In  the  freeboard,  some  particles  are  carried 
by  the  gas  flow  far  above  the  bed  surface  and  are  removed  from  the 
fluidized  bed  (elutriated)  ,  while  the  remainder  fall  back  to  the 
bed.   In  general  the  amount  of  bed  solids  suspended  in  the 
freeboard  ( entrainment )  decreases  exponentially  along  the  freeboard 
height.   This  distribution  is  similar  to  that  of  the 
Maxwell-Boltzmann  distribution  for  the  case  of  discrete  energy 
states  as  it  applies  to  the  Law  of  Atmospheres  IZ1 . 


N(z)  =  No  exp<-mgz/kT>  (1) 


Research  in  the  area  of  entrainment  by  Lewis  et  al .  C 3]  ,  Zenz 
and  Ueil  [4]  and  others  has  resulted  in  the  following  correlation 
for  entrainment  as  a  function  of  gas  velocity  and  freeboard  height 
for  small-particle  beds  [11. 


1G 


z 

=   B  exp<  -  [(b/Uo)   +  a  H  1  > 


uhere- 


flt  Uo  (Z) 


F  =  Entrainnent 

At  =  Area  of  bed 

Uo  =  Superficial  velocity  through  bed 
8  =  Particle  dependent  constant 
b  =  Particle  dependent  constant 
a  =  Particle  dependent  constant 
H  =  Height  of  freeboard 


fti  some  point  above  the  bed  surface,  the  quantity  of 
entrained  particles  becomes  constant.   At  this  point,  the  free 
fall  velocity  of  the  remaining  particles  is  equal  to  or  less  than 
the  uniform  superficial  operating  velocity.   The  height  at  uhich 
the  entrainment  becomes  constant  is  called  the  transport 
disengagement  height  (TDH)  [11. 

The  particles  that  are  throun  above  the  bed  are  affected  in 
the  freeboard  region  by  hydrodynamic  parameters  such  as^  bubble 
size,  bubble  frequency,  fluidizing  velocity,  height  above  the  bed, 
particle  size,  particle  density,  column  diameter  [53,  and 
baffles  [61.   The  intermittent  high  velocity  bursts  of  gas  uhich 
occurs  uhen  a  bubble  bursts,  imposes  a  flucuating  and  highly 
irregular  time  dependent  velocity  profile  over  the  cross  section 
of  the  bed  surface.   At  successively  higher  levels  above  the  bed 
surface,  this  velocity  profile  becomes  more  and  more  uniform  until 
at  the  TDH,  the  flou  is  at  the  uniform  superficial  operating 
velocity  ( Uo  )  [ 4] . 
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Until  recently,  little  attention  uas  paid  to  the 
understanding  of  the  freeboard  reactions  for  large  particle  beds. 
Houever,  due  to  recent  work  in  fluidized  bed  combustors,  the 
extent  of  particle  and  fines  loading  in  the  freeboard  has  been 
shoun  to  significantly  affect  the  S02  absorbtion,  NOx  reduction,  CO 
emission.   In  general,  the  last  5  to  10X  of  the  combustibles  will 
burn  in  the  freeboard.   It  uas  shoun  that  the  fine  sorbent 
particles  entrained  into  the  freeboard  uill  enhance  sulphur 
capture  and  that  entrained  char  particles  uill  react  uith  NOx  and 
reduce  its  emission  [71.   Fines  reinfection  has  been  shoun  to 
significantly,  increases  the  fine  particle  concentration  in  the  bed 
and  in  the  freeboard  uith  the  consequence  of  further  enhancing 
char  oxidation.   Houever,  this  can  result  in  overheating  .in  the 
freeboard  region  and  excess  SO2  and  NOx  emission  E71.   The  end 
result  is  that  the  potential  for  fluidized  bed  combustion  power 
plants  is  enhanced  by  their  ability  to  burn  high  sulphur  content 
fuels  and  maintain  lou  SO-,  emissions.   Further  research  in  the 
area  directly  above  the  bed  surface  is  required  to  properly  model 
the  reactions  occuring  uithin  the  freeboard. 

Extensive  studies  on  entrainment  rate  and  elutriation  have 

been  made  uith  numerous  correlations,  some  of  uhich  are  proposed 

in  CI,  3,  4,  5,  B,  7,  81.   Houever,  most  of  the  reported  work  on 

entrainment  from  fluidized  beds  has  been  carried  out  uith  either 

a  closely  sized  fraction  of  particles  or  a  mixture  of  tuo  such 

fractions.   Virtually  all  of  this  uork  has  been  conducted  on  bench 

scale  or  catalytic  cracking  fluidized  beds.   The  results  of  these 

18 


studies  have  been  shown  to  produce  much  louer  entrainment  rates 
than  full  sized  beds  or  tend  to  operate  in  the  slugging  condition. 
Therefore,  the  entrainment  rates  and  transport  disengagenent 
heights  (TDH)  for  fluidized  beds  are  generally  estimated  fropi 
empirical  or  semiempirical  correlations  obtained  from  this  data 
and  most  of  them  shou  extreme  discrepancies  betueen  different 
experimental  results.   Extrapolation  of  these  empirical 
correlations  usually  leads  to  strange  results  [51  with 
discrepancies  uhich  can  vary  by  two  orders  of  magnitude. 


The  lack  of  good  correlations  stems  mainly  from  the 


difficulty  in  obtaining  accurate  entrainment  rate  data.   Most  of 
the  data  is  based  on  pressure  measurements  at  incremental  heights 
in  fluidized  cracking  plant  type  beds  [111.   The  effect  of  uall 
loading  by  particles  and  the  actual  relation  betueen  pressure  and 
particle  concentration  is  considered  to  be  major  sources  of  error 
uhen  using  this  method  uith  large  particles.   As  a  result,  none  of 
the  correlations  are  uidely  accepted  as  giving  accurate 
predictions  [ 8] . 

n  complete  model  of  the  entrainment  process  from  fluidized 
beds  must  take  into  account  all  the  mechanisims  involved  uithin 
the  process.   The  arrival  of  bubbles  at  the  bed  surface,  ejection 
of  particles  from  the  dense-phased  bed  into  the  freeboard  region 
as  the  bubbles  erupt,  particle-particle  interactions,  and  the 
trajectories  of  ejected  particles  are    all  important  [8].  Much 

uork  has  been  done  concerning  bubble  grouth,  velocity,  volume, 
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etc-  and  their  behavior  is  fairly  well  understood. 

The  mechanism  of  solids  ejection  at  the  bubbling  bed  surface 
is  still  not  uell  understood.   The  origin  of  ejected  particles  is 
reported  to  be  primarily  due  to  tuo  sources.   The  particles  uhich 
have  been  lifted  by  the  bubble  uake  and  thrown  upuards  follouing 
the  bubble  burst  at  the  surfaceis  the-first  source.   This  theory 
is  supported  by  uork  done  by  George  and  Grace  [81  uho  performed 
experiments  uhich  concluded  that  the  vast  majority  of  the  ejected 
particles  did  not  originate  from  the  surface  layers  but  from 
bubble  uake  pick  up.   Uork  done  by. .Page  and  Harrison  [SI  also 
appears  to  agree  uith  this.   The  second  theory  suggests  that  the 
ejected  particles  originate  at  the  nose  of  the  bursting  bubbles 
and  are  throun  outuard  uhen  the  bubble  breaks.   Research  by  Roue 
and  Partridge  [8]  and  Glicksman  et  al  [12]  have  shoun  this  second 
mechanisim  as  being  the  dominate  particle  ejection  source  and  thus 
supporting  this  second  theory.   Their  uork  has  also  shoun  that 
under  the  conditions  in  uhich  2  bubbles  coalesce  just  belou  the 
surface  of  the  bed,  the  jet  of  gas  produced  can  result  in  a 
significant  amount  of  particles  being  ejected  from  the  uake  of 
the  first  bubble. 

The  effect  of  mul t iparticle  interactions  have  been  for  the 

most  part  ignored  except  by  Peters  and  Prybylouski  [13].   The 

motion  of  any  individual  particle  is  influenced  by  the  presence 

of  other  particles,  i.e.,  through  direct  particle-particle 

interactions  and  deviations  in  the  fluid  drag  force.   The  major 

20 


drawback  of  their  uork  is  that  the  paper  compares  their  theory 
uith  only  a  single  set  of  experimental  results  t 3] . 

Studies  to  model  the  trajectories  of  particles  in  the 
freeboard  have  been  conducted  several  times.   The  uork  of  Ualsh 
et  al  [71,  George  and  Grace  [81,  and  Peters  and  Prybylouski  [131 
name  just  a  feu  of  the  latest  efforts.   All  of  these  studies 
relied  upon  experimental  data  to  develope  their  theories. 
However,  to  check  the  accuracy  of  their  theories,  more 
experimental  data  is  required. 

As  of  yet,  none  of  the  entrainment  models  available  can  be 
incorporated  into  fluidized  bed  combustion  models  uith  sufficient 
accuracy  to  uarrent  their  use.   This  is  due  to  a  lack  of 
experimental  information  on  entrainment  rate  as  a  function  of  the 
complete  f luidization  parameters  of  the  bed  to  test  the  models 
uith.   As  a  result,  the  purpose  of  this  study  has  been  to  obtain 
particle  density  distributions  above  a  cold  atmospheric  fluidized 
bed  containing  a  continuous  particle  size  distribution  [Appendix 
El. 
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CHAPTER   II 

PARTICLE  SAMPLING  APPARATUS 
Desian  Alternatives 

There  are  nany  methods  available  for  determining  particle 
distributions  in  fluid  flous.   The  more  commonly  used  methods  an 


1)  Catching  mechanisms 

2)  Trapping  mechanisms 

3)  Radiation  attenuation  measurements 

4)  Optical  measurements 

5)  Capacitance  and  Inductance  measurements 


Catching  mechanisms  are  passive  devices.   That  is,  particles 
are  captured  merely  by  the  presence  of  the  catching  mechanism  in 
the  fluid  flou  containg  the  particles  to  be  sampled.   The  data 
obtained  using  this  method  is  position  dependent  and  produces 
average  values  for  the  particle  flux  loadings.   These  catching 
mechanisims  art   also  limited  in  that  they  can  only  catch  particles 
uith  particle  fluxes  traveling  in  a  single  direction.   The  device 
used  by  Ualsh  et  al  [10]  only  captured  falling  particles  uhiie  the 
device  used  by  George  and  Grace  [8]  required  the  uouard  moving 
particles  to  deflect  off  of  a  baffle  surface  and  fall  into  a 
collecting  trough. 
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Trapping  mechanisims,  unlike  catching  mechanisms,  are  active 
particle  samplers.   Their  operation  involves  the  trapping  and 
isolation  of  a  finite  volune  of  the  fluid  flou  at  a  specific 
period  in  tine.   This  sampling  technique  produces  time  dependent 
as  uell  as  position  dependent  data.   This  uill  allou  correlations 
between  bubble  eruption  and  particle  density  to  be  made  using 
multiple  bubble  conditions  rather  than  -single  bubble  capture.   As 
the  numPer  of  random  samples  taken  by  this  method  increases,  the 
average  value  of  this  data  uill  approach  that  of  the  catching 
mechanism.   Trapping  mechanisims  also  capture  particle  fluxes 


traveling  in  multiple  directions.-  This  ability  reduces  the  error 
inherent  in  measuring  only  the  dounuard  or  only  the  upuard 
particle  flux.   The  apparatus  used  in  this  paper  is  a  trapping 
mechanism. 

Attenuation  of  nuclear  particles  from  a  radioactve  source  can 
be  used  to  give  average  particle  density  distributions  across  a 
suspension.   Houever,  this  method  is  not  adaptable  to  density 
determinations  at  a  point.   An  average  time  dependent  density 
determination  can  be  achieved  uith  this  method.   Another  drawback 
of  this  method  is  the  radiation  hazards  involved  uith  the  use  of 
nuclear  particles. 

Optical  density  determinations  consists  of  tuo  seperate 
methods.   The  first  method  uses  a  very  small  light  beam  uhich  is 
eclipsed  by  the  transition  of  a  particle  through  it.   A  related 
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method  uses  the  absorbtion  and  scattering  of  a  someuhat  larger 
light  bean  to  correlate  the  change  in  light  intensity  uith 
particle  density.   This  method  has  been  used  frequently  in  the 
study  of  aerosols  but  requires  complicated  and  intricate  equipment 
[141.   The  second  method  involves  high  speed  photographs  of  a 
small  volume  of  space.   This  method  cannot  be  used  uhen  the 
particle  density^is  so  large  that  multiple  particles  eclipse  each 
other  frequently  enough  to  produce  unacceptable  error.   This  is 
the  case  uhen  the  probe  height  above  the  bed  is  less  than  7-15  cm 
(3-S  in). 


Measurements  at  a  point  can  also  be  made  by  inserting  either 


a  toroidial  inductor  or  a  parallel  plate  capacitor  in  the  flou. 
The  presence  of  the  particles  changes  the  permeability  and  thus 
the  inductance  of  the  inductor,  or  the  dielectric  strength  and 
thus  the  capacitance  of  the  capacitor.   The  drau-backs  of  these 
methods  involves  the  unknoun  effects  of  particle  velocity  and 
external  particles  on  the  inductor  and  charge  transfer  to 
particles  from  the  capacitor  [141. 

Apparatus  Requirements 

The  goal  of  this  study  uas  to  determine  the  density 
distribution  of  particles  above  an  atmospheric  fluidized  bed  uith 
particle  velocities  of  up  to  13  meters  per  second.   The  particle 
size  distribution  of  material  ejected  from  bubbles  is  required  for 
particle  trajectory  calculations.   A  correlation  betueen  the 


average  density  and  the  density  present  immediately  after  a  bubble 
bursts  from  the  bed  surface  uas  also  of  interest.   These 
requirements  dictated  that  the  method  used  for  measuring  densities 
have  the  follouing  capabilities^ 

1)  Measure  densities  uith  good  spatial  resolution. 

Z)  Measure  densities  at  specific  moments  in  time. 

3)  Obtain  particle  size  information. 

4)  Operate  under  extremely  dirty  conditions. 

5)  Easy  sample  removal  from  bed. 
S)  Remote  operation  of  sampler. 


The  radiation  attenuation  and  inductance/capacitance  methods 


can  not  determine  particle  sizes.   Therefore,  these  methods  uere 
no  longer  considered  as  possible  measurement  alternatives. 
Because  the  optical  methods  are  either  not  reliable  at  small 
heights  above  the  bed  or  their  use  is  too  complex,  they  uere  not 
used.   Catching  devices,  although  simple  to  use,  do  not  have  the 
ability  to  measure  data  at  specific  points  in  time  and  determine 
particle  density  loading  in  space.   As  a  result,  the  determination 
to  use  a  trapping  mechanism  as  the  method  of  measuring  particle 
densities  uas  made. 


ZS 


Apparatus  Design 


General  Design  Criteria 


The  follouing  criteria  uas  used  in  determining  the  design  of 

the  trapping  device. 


1)  The  closure  time  of  the  trap  uas  chosen  to  be  equal 
to  the  time  for  a  particle  with  a  velocity  of  10 
m/s  to  transit  1/10  the  length  of  the  sample 
container.   This  velocity  is  considered  to  be  the 
upper  limit  of  the  particle  velocity  distribution 
present  in  the  test  bed,  based  on  the  uork  of 
George  and  Grace  [81. 

2)  The  apparatus  must  be  capable  of  frequent  sampling 
uithout  requiring  access  to  the  sampling  device 

i tsel f . 

3)  The  samples  trapped,  must  be  easily  accessable 
from  outside  the  fluidized  bed  uithout  interrupting 
the  bed  conditions. 

4)  The  apparatus  must  be  able  to  operate  in  the  high 
particle  flux  environment  of  the  fluidized  bed. 


5)    The  actuation  of  the  trapping  device  must  be  able 
to  be  accurately  determined  to  allow  correlation 
uith  other  tine  resolved  measurements. 


General  Design 

The  ac-paratus  is  shoun  in  Figs.   Z  and  3.   A  description  and 
list  ot  all  components  is  given  in  Appendix  H.   The  sample 
container  is  mounted  on  an  extension  arm  to  minimize  the 
disturbance  to  the  air  flou  around  the  sample  trap  caused  by  the 
rest  of  the  mechanism.   The  sample  container  is  closed  using  tuo 
(Z)  paddle  arms,  one  above  and  the  other  belou.   These  paddle  arms 
are  attached  to  aluminum  interfaces  which  are  used  to  connect  them 
to  a  rotary  solenoid.   The  solenoid  is  used  to  suing  the  paddles 
over  the  sample  trap  and  shut  it.   Not  shoun  in  these  figures  are 
the  pouer  supply  for  the  solenoid,  the  vacuum  system  used  to 
remove  the  particles  from  the  sample  trap  and  the  uater-proof 
nylon  shell  used  to  keep  the  particles  from  interfering  uith  the 
operation  of  the  solenoid.   All  of  these  systems  are  described  in 
greater  detail  in  the  follouing  sections. 

Sample  Container 

To  ensure  that  the  sampling  device  had  minimal  effect  on  the 
fluid  flou,  the  cross  sectional  area  presented  to  the  flou  had  to 

be  minimized.   This  constraint  required  that  the  sample  container 
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be  seperated  from  the  rest  of  the  apparatus.   This  uas  benificial 
in  the  final  design  because  it  helped  to  reduce  the  apparatus 
closure  tine. 

To  ensure  a  good  seal  uas  achieved  uhen  the  trap  uas  shut, 
felt  uas  used  as  a  gasket  betueen  the  sample  container  and  the 
closure  paddles.   Fig.  4: shous  the  f inal  design  used  for  the 
sampling  mechanism.   It  is  constructed  of  1/1G  inch  aluminum  uith 
1/4  inc.-  square  stock  used  for  the  frame  and  mounting  structure. 
Epoxy  is  used  to  seal  the  sides  of  the  container. 

Closure  Paddles 

ftfter  several  iterations  on  paddle  designs,  it  uas  determined 
that  the  paddle  construction  uhich  offered  the  greatest  stiffness 
for  the  least  ueight  uas  a  composite  laminate.   The  paddle,  shoun 
in  Fig.   5,  is  made  using  a  0.4  in  thick  foam  core  uith  1/32  inch 
thick  Bassuood  laminations  on  the  exterior.   Hardwood  (Maple)  end 
pieces  uere  used  to  provide  a  noncompressi ve  connection  betueen 
the  foam  paddles  and  the  aluminum  interfaces.   The  aluminum 
interfaces  couple  the  solenoid  shaft  to  the  paddles.   Epoxy  uas 
used  to  join  the  laminate  materials. 

Actuator 


ft  rotating  mechanism  utilizing  a  rotary  solenoid  uas  chosen 

to  shut  the  sample  trap.   A  rotary  solenoid  uas  selected  because 
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of  the  problems  of  maintaining  lou  friction  surfaces  for  sliding 
mechanisms  in  the  presence  of  the  particulate  matter  being 
sampled.   This  precluded  the  use  of  any  sliding  mechanism  to  shut 
the  sample  trap.   The  calculations  found  in  Appendix  A  and 
Appendix  C  determined  the  size  of  the  solenoid  required  to  achieve 
the  desired  closure  time. 

The  45  degree  stroke  solenoid  was  chosen  to  place  the  paddle 
arms  far  enough  auay  from  the  sample  trap,  such  that  uhen 
de-energized  it  prevents  interference  uith  the  particle  flou. 
This  stroke  also  minimized  the  area  uhich  must  be  clear  of 
-obstructions  to  the  travel  of  the  paddles.   The  solenoid  operates 
at  a  1/10  duty  cycle  power  rating  uhen  initially  actuated, 
providing  the  torque  output  shoun  in  Fig.  S.   After  the  solenoid 
has  shut  the  sample  trap,  the  solenoid  is  operated  at  a  louer 
pouer  rating,  providing  a  holding  torque  of  5.5  in-lbs.   This 
decreased  rating  is  necessary  to  prevent  overheating  of  the 
solenoid.   This  assembly  is  encased  uithin  a  nylon  shell. 

Pouer  Supply 

Fig.   7  is  a  schematic  of  the  electrical  system  used  to  pouer 
the  solenoid.   Appendix  H  contains  a  list  of  all  components  used 
in  the  pouer  supply.   The  pouer  supply  plugs  directly  into  a 
standard  115  volt  AC  line  source.   Suitch  SI  is  used  to  apply 
pouer  to  the  solenoid  Ml.   The  full  uave  bridge  rectifier  assembly 

converts  the  AC  line  voltage  to  DC.   The  rectifier  assembly  is 
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Fig.  G      Torque  Output  Curve  for 
Rotary  Solenoid. 
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protected  by  the  arc  suppresser  01  from  the  large  voltage  spike 
induced  in  the  solenoid  coil  uhen  the  solenoid  is  actuated. 
Resistor  R3  is  used  to  reduce  the  current  through  the  solenoid 
after  the  solenoid  has  been  closed  to  prevent  it  from  overheating. 
To  initially  shut  the  sample  trap,  full  current  is  applied  to  the 
solenoid  by  shutting  suitch  SZ ,  uhich  bypasses  resistor  R3. 
Resistors  Rl  and  R2  form  a  voltage  divider  network  to  provide  a 
lou  voltage  (5.5  volts  AC)  trigger  source  for  an  oscilloscope. 
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Vacuum  System 

To  remove  the  particles  in  the  sample  trap,  a  vacuum  system 
uas  developed.   This  system  is  shoun  in  Fig.   8  with  a  list  of  the 
components  given  in  Appendix  H.   ft  vacuum  is  produced  by  allouing 
air  from  a  100  psi  air  source  to  flou  through  valve  VI  into  the 
venturi  eductor  PI.   The  vacuum  places  a  suction  on  the  sample 
container  CZ  via  a  fine  mesh  screen.   The  purpose  of  this  screen 
is  to  prevent  particles  from  escaping  the  sample  container.   The 
suction  is  applied  to  the  sample  trap  CI  through  1/4  inch  polyflou 
tubing.   It  is  through  this  tube  that  the  particles  are  removed 
from  the  sample  trap  and  collected  in  the  sample  container.   ftn 
equalization  and  agitaion  line  is  connected  to  the  opposite  side 
of  the  sample  trap.   This  line  serves  tuo  purposes.   First,  it 
ensures  that  the  vacuum  system  does  not  pull  in  particles  from 
outside  of  the  sample  trap.   Second,  it  allows  a  flou  of  air  to  be 
introduced  uhich  stirs  up  the  particles  trapped  inside.   This 
helps  push  them  into  the  suction  line  and  reduce  the  remaining 
particles  to  a  mininum. 

Apparatus  Test i no 

Tuo  tests  uere  run  to  determine  the  effectiveness  of  the 
system.   The  first  test  determined  the  closing  time  of  the  sample 
trap.   The  second  evaluated  the  error  from  the  loss  of  particles 
uhich  uere  left  in  the  sample  trap  by  the  vacuum  system. 
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The  procedure  and  method  used  to  determine  closure  tine  for 
the  sample  trap  is  given  in  Appendix  C.   From  the  results  of  these 
tests  given  in  Appendix  C,  the  closure  time  uas  determined  to  have 
an  average  value  of  1.44  milliseconds.   This  is  equal  to  a  10  m/s 
particle  traveling  1.44  cm  or  approximately  IS  Z  of  the  sample 
trap  length.   The  average  particle  velocity  uill  be  less  than  Z 
m/s  and  uill  therefore  introduce  an  average  error  of  less  than  3  Z. 

Using  the  sample  trap  closure  time,  a  dynamic  analysis 
determined  that  the  total  time  from  initially  applying  pouer  to 
the  solenoid  until  it  shut  the  sample  trap  is  42. 5  milliseconds. 
These  calculations  are  given  in  Appendix  D.   This  actuation  time 
is  important  for  determining  uhether  a  specific  bubble's  debris 
uas  uithin  the  vicinity  of  the  trap  at  the  time  of  closure. 

The  procedure  and  results  for  the  testing  of  the  sample 
removal  vacuum  system  are  given  in  Appendix  8.   The  average  amount 
of  particles  lost  by  the  vacuum  system  uas  determined  to  be  0.52  Z 
of  the  initial  sample  placed  in  the  trap.   The  maximum  error  uas 
0.93  Z.   The  error  from  the  vacuum  system  is  therefore  considered 
to  have  an  insignificant  effect  on  the  data  obtained. 
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CHAPTER   III 

EXPERT MENTAL  PROCEDURE 

Fluidized  3ed  Configuration 

The  M.I.T.   atmospheric  fluidized  bed,  in  uhich  the  sampling 
device  uas  used,  is  a  model  of  the  20  MU  atmospheric  fluidized  bed 
combustor  prototype,  jointly  sponsored  by  the  Tennessee  Valley 
Authority  and  the  Electric  Pouer  Research  Institute-   The 
fluidized  bed  model  is  described  in  Lord  et  al  [15]  and  Jones  et 
al  [1S1.   Uhile  using  the  sampling  apparatus,  a  different  heat 
exchanger  tube  bundle  configuration  uas  used  than  is  described  in 
Jones  et  al  [IS].   The  heat  exchanger  configuration  used  is  shoun 
in  Figs.   9  and  10. 

The  heat  exchanger  used  during  this  work  is  made  of  l.ZS  cm 
(0.5  in)  O.D.   tubing  arranged  in  4  rous  of  ZZ  pipes  each.   The 
tubes  are  aligned  as  shoun  in  Fig.   10.   Each  pipe  is  spaced  uith 
a  vertical  center  to  center  distance  of  5.08  cm  (Z  in)  and  a 
horizontal  center  to  center  distance  of  3.91  cm  (1.5  in).   The 
distance  from  the  distributor  to  the  center  of  the  upper  most  tube 
is  Z7.5Z  cm  (10.875  in).   A  5.08  cm  (Z  in)  spacing  separates  the 
front  and  back  ualls  of  the  fluidized  bed  from  the  end  tubes  of 
the  bundel.   The  cross  sectional  area  oi    the  bed  is  1.079  sq  m 
(11.61  Sq  Ft).   The  particulate  material  used  in  the  bed  is  a 
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Fig.  9      Heat  Exchanger  Tube  Design  Showing  The 
Four  (4)  Rous  of  ZZ  Tubes. 
All  Dimensions  in  en. 
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Fig.  10     Hea:  Excnangsr  Tube  Design. 
All  Olnensions  in  en. 


nixture  of  steel  grit  abrasive  having  a  specific  gravity  of  8.1. 
Appendix  E  lists  the  size  distribution  of  the  steel  grit  used 
during  the  sampling  operation.   The  bed  uas  operated  uithout 
recycling  the  fines  captured  in  the  cyclones.   The  static  bed 
height  of  the  material  uas  22.54  cm  (8.875  in)  throughout  the  data 
collection  period. 

Equipment  Set-ua 

Fig.  11  shous  the  placement  of  the  sampling  device  inside 
the  fiuidized  bed.  The  sample  trap  uas  positioned  such  that  it 
uas  directly  over  a  spacing  betueen  tubes  [22.86  cm  (9  in)  from 
the  center  of  the  front-most  tube]  and  41.91  cm  (IS. 5  in)  from  the 
left  uall.  The  height  of  the  sample  trap  above  the  distributor 
uas  varied  during  the  sampling  process  as  discussed  in  the  section 
on  sampling  procedure. 

Fig.   12  shous  the  placement  of  the  bubble  probe  and  the 
anemometer  probe  uith  respect  to  the  sample  trap.   The  bubble 
probe  uas  placed  directly  belou  the  sample  trap  and  2G.G7  cm  (10.5 
in)  above  the  distributor.   The  probe  extension  uas  placed  at  an 
angle  so  as  not  to  interfere  uith  the  sample  trap  operation.   To 
protect  the  anemometer  uire  from  particles  impacting  it,  a  special 
shield  consisting  of  3  323  nesh  screen  and  an  aluminum  frame  uas 
placed  around  it.   The  anemometer  probe  uas  attached  to  the  bubble 
probe  extension  uith  the  entrance  to  the  anemometer  probe  29.21  cm 
(11.5  in)  above  the  distributor.   This  placed  the  entrance  to  the 
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Fig.  11     Position  of  Sample  Trap,  Bubble  Probe, 

and  flnenonetar  Probe  Above  Heat  Exchanger 
Tubes.   All  Dimensions  in  cm. 
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Fig.  12     Position  of  Bubble  and  Anemometer 
Probe  uuth  Respect  to  the  Sample 
Trap  and  Distributor.   All  Dimensions 
are  in  en. 

45 


anemometer  probe  Z.S4  en  (1.0  in)  auay  from  the  center  of  the 
sample  trap  and  0.S4  cm  (0.25  in)  outside  the  area  directly  belou 
the  sample  trap.   As  a  result,  the  probe  has  a  minimum  effect  on 
the  air  flou  from  the  bed  to  the  sample  trap,  but  uill  only 
measure  the  gas  velocity  at  the  edge  of  the  sample  trap  perimeter. 
The  positioning  of  the  tuo  probes  above  the  distributor  remained 
constant  throughout  the  sampling  evolution- 
Fig.   13  is  a  block  diagram  showing  the  equipment  used  during 
the  sampling  operations  and  their  interconnections  .   Table  1  is  a 
listing  of  the  equipment  used.   The  oscilloscope  time  base  uas  set 
for  MANUAL  TRIGGER,  SINGLE  SUEEP  mode  and  a  sueep  time  of  50 
ms/div.   The  channels  of  the  dual  trace  amplifier  uere  set  at  5 
volts  per  division  for  the  bubble  probe  and  2  volts  per  division 
for  the  sample  trap  inputs.   On  the  differential  amplifier,  one 
channel  uas  not  used  and  the  second  channel  uas  set  at  1  volt  per 
division  for  the  anemometer  probe  input. 
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Fig.  13     Block  Diagram  of  the  Equipment  Used 
During  the  Sampling  Operations. 
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TABLE  1 


Oscilloscope  TEKTRONIX  5111  Storage  Oscilloscope 

TEKTRONIX  SA18N  Dual  Trace  Ampl 
TEKTRONIX  5AZ1N  Differential  Ampl 
TEKTRONIX  5B10N  Tine  Base 

Anemometer  Thermal  Systems  Inc. 

1851-Z  Monitor  and  Pouer  Supply 
10S4-A  Linearized  Anemometer 
Anemometer  uire  u/  #323  screen  guard 

Optical  Signal  Detector 

Optical  Signal  Detector  Pouer  Supply 

Optical  Source  and  Pouer  Supply 

Oscilloscope  Camera 


Listing  of  equipment  used  during  particle  sampling 
operations. 
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To  determine  the  f luidization  conditions  uithin  the  bed,  a 
set  of  manometers  uere  used.   These  manometers  measured  pressures 
uithin  the  bed,  at  heights  from  4.13  cm  <i.S2  in)  to  37.15  cm 
(14.SZ  in)  above  the  distributor,  in  Z . 54  cm  (1  in)  increments. 
The  pressure  data  corresponding  to  each  trap  position  and  bed 
velocity  is  listed  in  Appendix  I. 

To  determine  the  gas  flou  conditions  uithin  the  bed,  an 
orfice  flow  meter  uith  ID  -  1/ZD  taps  uas  located  upstream  of  the 
distributor.   The  computer  program  listed  in  Appendix  G  uas  used 
to  convert  the  pressure  tap  data  to  mean  air  velocities  uithin  the 
bed. 

Sampling  Procedure 

Data  uas  collected  for  four  (4)  mean  bed  velocities  at  six 
(6)  different  sample  heights.   The  sample  trap  uas  placed  at  a 
given  sampling  height  (measured  from  the  distributor  to  the  bottom 
of  the  trap),  and  ten  (10)  samples  uere  collected  at  each  desired 
velocity.   The  trap  position  uas  then  changed  to  a  new  height. 

During  certain  sampling  conditions,  those  uhich  involved  lou 
sampling  heights  uith  the  higher  air  velocities,  the  paddle  arms 
uould  occasionally  impact  the  sides  of  the  sample  trap  and  not 
close  the  sample  trap  completely.   It  is  assumed  that  this  occured 
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uhen  a  large  bubble  erupted  directly  under  the  paddle  arn  and 
deflected  the  paddle  arn  into  the  side  of  the  trap.   Whenever  this 
occured,  the  trap  uas  de-energized  and  the  closure  cycle  repeated. 

Each  time  the  height  or  velocity  uas  changed,  a  complete 
sampling  cycle  uas  conducted  and  this  sample  discarded.   This  uas  to 
prevent  any  accumulation  of  particles  (in  the  entrance  to  the 
vacuum  or  purge  lines  on  the  sample  trap)  from  being  added  to  the 
first  sample  at  the  neu  height  or  velocity. 

For  each  set  of  data  at  a  given  height  and  velocity  the 
follouing  information  uas  recorded^ 

1)  Fluidized  bed  height  above  the  distributor 

determined  visually  and  by  pressure  measurements. 
Z)  Pressure  upstream  of  orfice  plate  (PI) 

3)  Pressure  difference  across  orfice  plate  (delta  P) 

4)  Air  temperature  in  bed 

5)  Pressure  distribution  in  bed 

6)  Height  of  sample  trap  above  distributor 

The  follouing  procedure  uas  used  during  sampling-- 

For  each  sample  to  be  taken  uithin  a  data  set. 

ft.  I m -rial  conditions 

1)  Vacuum  air  supply  OFF 

Z)  Sweep  trigger  on  oscilloscope        RESET 

5a 


3)  O-scope  memory  ON-CLEAR 

4)  Solenoid  pouer  supply  switch  SI      OFF 

5)  Solenoid  trigger  suitch  SZ  ON 


B.  Sampling  Procedure 


1)  Trigger  oscilloscope  sueep  and  uait  until  sueep 
is  at  the  center  of  the  CRT. 

2)  Close  the  Solenoid  pouer  supply  suitch  SI.   Uhen 
SI  is  shut,  the  oscilloscope  will  shou  an 
additional  trace.   This  third  trace  is  used  to 
determine  the  closure  time  relative  to  the 
presence  of  gas  jets  and  bubble  eruptions.   An 
example  of  a  typical  oscilloscope  trace  is  shoun 
in  Fig.  14. 

3)  Open  the  solenoid  trigger  suitch  SZ.   This 
reduces  the  current  to  the  solenoid  to  prevent 
overheating.   The  maximum  alloued  time  to  let  SZ 
remain  closed  is  five  (5)  seconds. 

4)  Turn  on  the  air  supply  to  the  vacuum  system  and 

leave  on  one  (i)  minute.   The  exhaust  air  from 

the  ventun  on  the  vacuum  system  must  be  directed 

into  the  purge  line  in  an  oscillatory  manner. 

This  uill  agitate  the  particles  uithin  the  trap 
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Fig.  14     Typical  Oscilloscope  Trace  Obtained 
During  Sanpling  Operation. 
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so  as  to  move  then  into  the  vacuum  line. 
5)  Turn  off  the  vacuum  air  supply. 
S)  Turn  off  the  solenoid  power  supply  suitch  SI. 

7)  Remove  sample  from  sample  container  and  place  in 
specimen  bottle. 

8)  Photograph  trace  on  oscilloscope. 

9)  Repeat  from  ft. 


Sample  Analysis 

Sample  Ueight  Determination 

Each  sample  obtained  uas  ueighed,  using  a  Torsion  Balance  Co. 
TORBflL  scale,  to  an  accuracy  of  0.01  grams.   The  average  value  and 
standard  deviation  uas  then  determined  for  each  set  (specific 
height  and  velocity)  of  samples.   The  weight  of  particles  in  a 
completely  filled  sample  trap  uas  also  determined  for  void 
determination.   These  results  are  listed  in  Appendix  I. 
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Determination  of  Particle  Size  Distribution 

To  determine  the  particle  size  distribution  uhich  occurs  at 
various  heights  for  a  specified  gas  velocity,  three  (3)  samples 
from  each  set  of  data  at  a  specified  height  and  velocity,  and  frorc 
the  bulk  bed  material  uere  analyzed  using  a  Zeiss  Videoplan  Z 
Image  analyzer.   The  software  used  uas  the  "Image  Analysis  System 
MOP-Videoplan"  distributed  by  KONTRON  Electronics  Group. 

For  each  sample  selected,  a  microscope  slide  uas  prepared 

using  double  sided  adhesive  tape  on  uhich  a  portion  of  the 

selected  sample  uas  placed.   To  ensure  a  sharp  contrast  uas 

achieved,  each  slide  uas  backed  uith  uhite  paper.   The  sample  uas 

then  placed  under  a  Zeiss  microscope  to  uhich  the  image  analyzer 

uas  connected  via  a  vidiocon  tube.   The  magnification  used  uas 

53X  uhich  provided  an  average  vieu  of  about  8  particles  at  a 

time.   The  analyzer  uas  then  used  to  analyze  the  partices  present 

on  each  slide  of  interest.   The  mode  in  uhich  the  image  analyzer 

uas  used  uas  the  equivalent  diameter  mode.   This  mode  determined 

the  cross  sectional  area  of  each  particle  sampled  and  calculated 

the  diameter  of  a  circle  uith  the  same  area.   Uhen  a  data  set  uas 

completed,  the  data  uas  analyzed  for  particle  size  distribution. 

The  resulting  output  (Appendix  J)  consists  of  a  particle  count  vs 

particle  size  histogram,  a  cumulative  frequency  plot  and  a 

classification  data  list.   ft  gaussian  distribution  fit  for  the 

data  is  also  plotted  on  the  histogram  and  cumulative  frequency 
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plots.   The  X  axis  of  the  plots  are  normalized  uith  a  range  of 
zero  (0)  to  four  (4).   To  determine  the  actual  diameter  of  the 
particle  in  microns  for  a  given  normalized  value,  the  normalized 
value  must  be  multiplied  by  the  conversion  factor  Z48.887. 
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CHAPTER   IV 


COMPUTER  MODEL 
Introduction 

This  chapter  descibes  the  theory,  logic  flou  and  testing  of  a 
particle  trajectory  computer  model  to  predict  the  particle  loading 
in  the  freeboard  due  to  erupting  buubles.   The  program,  listed  in 
Appendix  L,  is  written  in  HP  BASIC  2.0  and  uas  run  on  a  Hewlett 
Packard  981S  microcomputer.   The  output  from  the  model  is  discussed 
in  chapter  VI  and  compared  uith  experimental  results.   An  in  depth 
analysis  of  the  program  logic  and  structure  is  given  in  Appendix  L. 

Model  Theory 

The  model  developed  here,  is  based  on  calculating  the 
trajectory  of  a  single  particle  as  it  is  ejected  from  the  bed 
surface  and  is  acted  upon  by  gravitational  and  drag  forces.   The 
drag  force  is  due  to  the  difference  in  absolute  particle  and  air 
velocities.   The  air  velocity  is  a  combination  of  the  initial  jet 
velocity  produced  uhen  a  bubble  bursts,  and  the  superficial  bed 
velocity.   To  ensure  that  the  particle  drag  is  calculated 
accurately,  the  following  drag  coefficient  correlation  given  by 
Uhite  [17]  uas  used. 
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Cd  =  Z4   +     6     +  0.4 

Re     1+  /Re-1  (3) 


where: 

Cd  =  Drag  coefficient  for  sphere 
Re  =  Reynolds  number 


Eqn  3  is  valid  over  the  range  0<Re<10  .   To  calculate  the 
particles  position  and  velocity,  the  computer  uses  a  forward 
difference  method.   Using  Neuton's  Lau  (ZF=ma)  the  acceleration  of 
the  particle  due  to  gravity  and  drag  is  determined.   Inserting 
this  acceleration  into  Eqn  4,  the  particles  neu  velocity  is 
determined. 


V  =  Vo  +  a  t  (4) 


uhere: 

V  =  Particles  neu  velocity 
Vo  =  Particles  present  velocity 
a  =  Acceleration  of  particle 
t  =  Time  increment  of  calculation 


To  determine  the  particles  neu  position,  the  velocity 

calculated  in  Eqn  4  is  inserted  into  Eqn  5. 
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H  =  Ho  +  (Vo  +  V)  t 

Z  (5) 


uhere- 

H  =  Particles  neu  height 
Ho  =  Particles  present  height 
Vo  =  Particles  present  velocity 
V  =  Particles  neu  velocity 
t  =  Tine  increment  of  calculation 


Using  Eqns  3,4, and  5,  the  trajectory  of  the  particle  is  calculated 
from  the  tine  it  initially  leaves  the  bed  until  the  time  that  it 
returns  to  the  bed. 

These  calculations  are  repeated  over  a  range  of  particle 
diameters  from  80  to  570  nicrons.   By  determining  the  residence 
time  of  each  particle  uithin  a  specified  height  increment  (A  H) 
above  the  surface  of  the  bed,  a  particle  density  distribution 
above  the  bed  is  determined.   The  height  increment  (AH)  used  in 
the  program  is  Z  cm.   At  the  end  of  each  time  step  uhen  the  height 
calculation  ( Eqn  5)  is  completed,  the  counter  representing  the 
particular  Z  cm  height  increment  which  the  particle  is  in,  is 
incremented  by  one.   Each  particle  size  has  its  oun  set  of 
counters  to  alloy  individual  particle  analysis. 

The  calculated  density  distribution  is  then  weighted  uith  the 
particle  size  distribution  of  the  bulk  bed  material  since  the 
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probability  of  a  given  size  particle  being  present  at  a  specified 
height  is  dependent  upon  the  number  of  particles  uithin  the 
system.   This  is  accomplished  by  multiplying  each  height  counter 
of  a  given  particle  size  with  the  number  of  particles  for  that 
given  size  present  in  the  input  bed  distribution.   By  summing  the 
density  values  for  each  set  of  particle  diameters  at  a  given 
height  over  the  entire  freeboard  of  the  bed,  the  overall  particle 
density  above  the  bed  surface  is  determined. 

The  model  assumes  that  all  the  particles  are  ejected 
perpendicular  to  the  surface  of  the  bed  and  are  initially  at  a 
uniform  velocity.   Because  the  model  uses  single  particles  for  the 
analysis,  the  effects  of  mul tiparticle  interactions  are  not 
included  in  the  model. 

Testing  of  Program 

To  evaluate  the  validity  of  the  program,  tuo  tests  were 
conducted.   The  first  test  compared  the  height  solution  produced 
by  the  computer  uith  a  closed  form  solution.   The  second  test 
involved  running  the  program  uith  different  particle  diameter  and 
particle  distribution  height  intervals  to  ensure  that  a  valid 
sample  size  uas  being  used. 

Closed  Form  Solution 

To  determine  a  closed  form  solution  for  particle  height  as  a 
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function  of  initial  particle  velocity  and  superficial  velocity  Uo, 
a  force  balance  uas  used.   The  forces  acting  on  a  particle  are 
gravitational  and  drag.   The  gravitational  force,  Fl,  is  simply 
the  volume  of  the  particle  multiplied  by  the  particles  density  and 
the  gravitational  acceleration,  and  can  be  uritten  as: 


Fl  =  -   Pa1T  D   g 

s  <s> 


uhere: 

Fl  -  Force  due  to  gravity 
Pp  =  Density  of  particle 
D  =  Diameter  of  particle 
g  =  gravitational  acceleration 


In  order  to  get  a  closed  form  solution  that  did  not  involve 
non-linear  differential  equations,  Stokes  flou  uas  used  for  the 
closed  form  solution  only.   The  computer  model  used  the  Stokes 
equation  only  to  compare  results  uith  the  closed  form  solution, 
afteruard,  Eqn  3  uas  used.   Using  the  Stokes  drag  coefficient 
relation,  the  drag  force  on  a  particle  can  be  determined  as: 


uhere 


FZ  =  3  (Uo  -  Up)  »    D  (7) 


FZ  =  Force  due  to  drag 
Uo  =  Superficial  bed  velocity 
Up  =  Velocity  of  particle 
l*    =  Absolute  viscosity  of  air 
D  =  Diameter  of  particle 


B0 


By    inserting   Eqns    S   and  7    into   Neuton's   Lau   (£  F=ma>    and 
simplifying,    the   following   differential    equation    is    obtained: 


uhere^ 


boundary  conditions^ 


X  +  CI  X  =  C2  (8) 

X  =  Acceleration  of  particle 

X  =  Velocity  of  particle 

CI  =  18^ 
Pp  Dz 

CZ  =  CI  Uo  -  g 


1)  t=0    X=Uo 

2)  t=0   X=0 


This  second  order  linear  differential  equation  can  be 
solved  using  the  given  boundary  conditions  uith  the  resulting 
closed  form  solution  given  as: 


Gl 


X  =  C3  +  CZ   t  -  C3  exp<  -CI  t  } 
CI 


(9) 


where- 


X  =  Height  of  particle 
C3 


Uo  -  CZ 

CI 


JL 

ci 


Using  the  initial  conditions  listed  in  table  Z,  the  solution 
obtained  using  the  computer  model  (maximum  height=  15.Z70  cm,  time 
to  maximum  height=  0.Z75  sec)  was  identical  to  three  decimal 
places  uith  the  solution  obtained  using  Eqn  9. 


TABLE  Z 


Superficial  Velocity  ( Uo ) :  G0.96  cm/s  (Z  ft/s) 

Initial  particle  velocity  ( Upo ) :  304.8  cm/s  (10  ft/s) 

Particle  diameter:  Z00 

Particle  density: 

Time  increment  (computer):  0.001  sec 

Jet  velocity:  0.0    cm/s 


microns 
5000   kg/m3 


List  of  parameters  used  to  check  computer 
calculations  against  closed  form  solution. 
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Sample  Size  Sensitivity  Test 

To  ensure  that  appropriate  sample  sizes  uere  used  to  minimize 
errors  due  to  coarse  sampling  intervals,  tuo  sensitivity  tests 
uers  run.   One  test  involved  changing  the  particle  diameter 
interval  from  10  microns  to  5  microns.   The  second  test  changed 
the  height  sampling  interval  (A  H)  from  2  cm  to  5cm  and  then  1  cm. 
For  each  of  the  tests,  the  same  initial  conditions  uere  input  into 
the  program.   Table  3  shous  the  resulting  output  from  the  program 
listing  the  initial  conditions  and  the  resulting  diameter  versus 
maximum  height  data  calculated.   Fig.   15  shous  a  plot  of  the 
calculated  maximum  height  vs  particle  diameter  data  listed  in 
table  2.   Fig.   IS  shous  the  bed  particle  distribution  used  for 
each  of  the  tests. 

Fig.   17a  shous  the  entrainment  calculation  using  a  diameter 

interval  of  10  microns  and  a  A  H  of  2  cm.   Fig.   17b  shous  the 

same  calculation  using  a  diameter  interval  of  10  microns  and  a  A  H 

of  5  cm.   The  curve  is  not  as  smooth  but  still  retains  the  same 

general  shape.   The  peak  of  the  curve  shoun  in  Fig.   17b  occurs  at 

a  height  of  about  20  cm  uhereas  the  peak  in  Fig.   17a  occurs  at 

about  IS  cm.   n  semi-log  plot  of  these  curves  uould  shou  that  the 

sloce  of  the  line  to  the  right  of  the  peak  uould  be  larger  for  the 

data  represented  by  Fig.   17b.   The  effect  of  maintaining  A  H  at  5 

cm  but  decreasing  the  diameter  interval  to  5  microns  is  shoun  in 

Fig.   17c.   There  is  no  readily  detectable  difference  betueen  Fig. 
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cm/s 
Mean  Bed  Velocity-  57.912 

Initial  Particle  Velocity3   97.2312 
Peak  Jet  Velocity3  609.5 

Ga5  Jet  Duration3  .02  s 

Maximum   Max  Ht 

Diameter   Heiaht    Tine       Diameter 
urn       cm      seconds        urn 


Maximum 

Heiaht 
cm 


Max  Ht 

Time 

seconds 


30 
100 
120 
140 
1S0 
180 
200 
220 
240 
2B0 
280 
300 
320 
340 
3G0 
380 
400 
420 
440 
460 
480 
500 
523 
540 
553 


55.1 
57.2 
55.8 
52.7 
48.8 
44.7 
40.7 
37.0 
33.7 
30.8 
28.2 
25.1 
24.2 
22.5 
21.1 
19.8 
18.7 
17.7 
16.9 
16.1 
15.4 
14.3 
14.2 
13.7 
1  -.    i 


.232 
.297 
.305 
.307 
.303 
.295 
.285 
.275 
.264 
.254 
.245 
.236 
.228 
.220 
.213 
.207 
.201 
.196 
.191 
.187 
.183 
.179 
.176 
.173 
.170 


90 
110 
130 
150 
170 
190 
210 
230 
250 
270 
290 
310 
330 
350 
370 
390 
410 
430 
450 
470 
490 
510 
530 
550 
570 


57.0 
55.8 
54.4 
50.3 
46.8 
42.7 
38.8 
35.3 
32.2 
29.5 
27.1 
25.1 
23.3 
21.3 
20.4 
19.2 


18. 

17, 
16, 
15, 


15.1 

14.5 
14.0 
13.5 
13.3 


.290 
.302 
.307 
.305 
.299 
.290 
.280 
.270 
.259 
.249 
.240 
.232 
.224 
.217 
.210 
.204 
.199 
.194 
.189 
.185 
.181 
.178 
.174 
.171 
.158 


Table  3     Listing  of  Input  and  Resulting  Maximum 
Particle  Heights  with  Time  to  Maximum 
Height.   Thes  Values  Uere  Used  [Juring 
the  Increment  Sensitivity  Tests. 
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Fig.  17     Plots  of  Relative  Particle  Density  vs  Freeboard 
Height  shouin  the  effects  of  varying  Diameter 
Interval  and  Height  Interval. 
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17c  and  Fig.   17b.   In  Fig.   17d,  the  effect  of  changing  A  H  to  1 
en  and  diameter  spacing  to  5  microns  is  shown.   The  difference 
betueen  Fig.   17a  and  Fig.   17d  is  barely  noticable  and  no 
detectable  change  in  the  slope  to  the  right  of  the  peaks  is 
present.   As  a  result  of  this  analysis,  the  program  uas  operated 
uith  a  A  H  of  2  cm  and  a  diameter  interval  of  10  microns.   This 
reduced  the  calculation  time  to  half  of  that  required  uhen  using  a 
diameter  spacing  of  5  microns  uhen  the  same  total  diameter  span 
uas  used  and,  as  uas  seen  in  Figs.   17a  and  17d,  the  difference  in 
output  does  not  require  the  finer  increment. 
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CHAPTER   V 


EXPERIMENTAL  RESULTS  FIND  DISCUSSION 
Minimum  Fluidization  Velocity 

Fig.   18  is  a  plot  of  mean  bed  velocity  (Uo)  vs  pressure  drop 
through  the  bed  (Pb).   The  velocity  values  uere  determined  using 
the  program  in  Appendix  G  and  the  pressure  data  listed  in  Appendix 
I.   From  this  plot,  the  minimum  fluidization  velocity  ( Umf )  for 
the  bed  conditions  used  during  this  study  is  determined  to  be  0.15 
m/s  (0.5  Ft/sec). 

Entrainment  Analysis 

Table  4  lists  the  averaged  sample  weights  and  their  standard 
deviations  for  the  samples  (Appendix  I)  collected  by  the  sampling 
apparatus.   Looking  at  the  standard  deviation  of  the  sample 
groups,  the  standard  deviation  is  fairly  large  compared  to  the 
average  values.   Houever,  visual  observations  of  the  fluidized  bed 
in  operation  would  suggest  that  a  larger  standard  deviation  uould 
be  expected.   The  short  sample  cycle  (2  ms )  and  the  bubble  burst 
activity  in  the  bed  are  the  main  reasons  for  this  conclusion.   The 
average  standard  deviation  is  Z9  7.    of  the  average  sample  weight 
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and  the  range  extends  from  68  Z   of  the  average  for  the  sample  set 
numbers  10S  thru  115  to  only  7  X   of  the  average  for  the  samples 
285  thru  295.   This  variation  is  believed  to  be  due  to  the  nature 
of  the  bubble  activity  uithin  the  bed.   It  uould  appear  that  uhen 
a  bubble  erupts  at  a  time  just  prior  to  closure  of  the  trap,  a 
sample  ueight  greater  than  the  average  uould  be  obtained. 
Houever ,  uhen  no  bubble  has  been  present,  the  material  caught  by 
the  trap  should  consist  mainly  of  particles  returning  to  the  bed 
and  therefore  be  belou  the  average  ueight  caught.   Table  4  also 
shous  the  density  of  each  sample  collected.   These  values  uere 
determined  by  dividing  the  averaged  sample  weights  by  the  volume 
of  the  sample  trap  uhich  is  43.02  cubic  cm. 
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TABLE  4 


Superficial  Bed  Velocity  ( Uo ) : 
Non-dimensionalised  velocity  (Uo/Umf) 


58.1  cm/s  (1.905  ft/s) 
3.81 


Sample 
Numbers 


28S-Z35 

24S-255 

11B-12S 

IBS- 175 

7B-85 

176-185 


Height 
Above  Bed 
(en) 


3. 

7. 

12. 


17.8 


22 
31, 


Average 
Ueight 
( grams ) 


7.22 
4.48 
2.41 
1.08 
1.32 
0.31 


Standard 

Average 

Deviation 

Density 

( grams ) 

(grams/cm  ) 

0.51 

1.68  (-01) 

0.72 

1.04  (-01) 

0.7B 

5.60  (-02) 

0.16 

Z.51  (-02) 

0.23 

3.07  (-02) 

0.10 

7.21  (-03) 

Superficial  Bed  Velocity  (Uo): 
Non-dimensionalised  velocity  (Uo/Umf) 


48.3  cm/s  (1.585  ft/s) 
3.17 


Sample 

Height 

Average 

Standard 

Average 

Numbers 

Above  Bed 

Ueight 

Deviation 

Oensity 

(cm) 

( grams ) 

( grams ) 

(grams /cm  ) 

276-285 

5.1 

4.35 

1.05 

1.01  (-01) 

236-245 

8.9 

1.88 

0.22 

4.37  (-02) 

96-105 

13.6 

1.36 

0.32 

3.16  (-02) 

156-165 

19.1 

0.60 

0.12 

1.39  (-02) 

56-65 

24.1 

0.92 

0.35 

2.14  (-02) 

186-195 

31.8 

0.18 

0.04 

4.18  (-03) 

Superficial  Bed  Velocity  (  Uo  )  : 
Non-dimensional ised  velocity  (Uo/Umf) 


39.5  cm/s  ( 1.297  ft/s) 
2.59 


Sample 

Height 

Average 

Standard 

Average 

Numbers 

Above  Bed 

Ueight 

Deviation 

Density 

(cm) 

( grams ) 

( grams ) 

(grams/cm  ) 

265-Z75 

5.7 

1.31 

0.37 

3.05  (-02) 

225-235 

9.5 

0.S2 

0.17 

1.44  (-02) 

126-135 

14.3 

0.28 

0.08 

5.51  (-03) 

136-145 

19.7 

0.08 

0.04 

1.86  (-03) 

86-95 

24.8 

0.09 

0.05 

2.09  (-03) 

19G-205 

32.4 

0.06 

0.02 

1.39  (-03) 
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TABLE  4  (cont) 


Superficial  Bed  Velocity  (Uo):  35.4  cm/s  (1.161  ft/s) 

Non-dimensionalised  velocity  (Uo/Umf):    2 . 3Z 


Sample 

Height 

Average 

Standard 

Average 

Numbers 

Above  Bed 

Ueight 

Deviation 

Density 

(cm) 

(grams) 

( grams ) 

(grams /cm  ) 

255-255 

S.4 

0.63 

0.13 

1.45  (-02) 

Z1S-225 

10.2 

0.32 

0.04 

7.44  (-03) 

106-115 

14.9 

0.22 

0.15 

5.11  (-03) 

146-155 

20.3 

0.02 

0.01 

4.55  (-04) 

S6-75 

26.0 

0.12 

0.04 

2.79  (-03) 

20S-215 

33.0 

0.01 

0.004 

2.32  (-04) 

List  of  experimental  data  shouing  sample  averages,  standard 
deviations,  heights,  and  velocity  conditions  measured. 
Density  values  are  calculated  by  dividing  the  average  sample 
ueight  by  the  sample  trap  volume. 


Fig.   19  shous  the  relationship  between  the  average  particle 
density  caught  in  the  sample  trap  and  the  sample  trap  height  above 
the  bed  surface  as  a  function  of  Uo/Umf. 


For  the  relatively  lou  fluidization  velocities  used  during 
the  data  measurements,  (maximum  Uo/Umf  =  3.81)  the  freeboard 
height  can  be  assumed  infinite.   Using  this  assumption  implies 
that  all  of  the  particles  return  to  the  bed  and  none  are 

elutriated,  le  comclete  reflux.   Under  these  conditions,  the 
following  equation  has  been  suggested  by  Leuis  et  al  i 33  and  Kunii 
and  Levenspiel  [11  to  model  the  particle  loading  uithin  the 
freeboard. 
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Fig.  19  Plot  of  Particle  Density 
as  a  Function  of  Uo/Unf. 
Obtained  Using  the  Particle  Sampler. 
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vs  Freeboard  Height 
This  data  was 


Pr  =  Po  expOa  1>  (  10) 


Pr  =  Particle  density  at  height  1 

Po  =  Particle  density  at  (+)bed  surface 

a  =  Characteristic  particle  decay  length 

1  =  Height  above  bed  surface 


Table  5  shous  the  values  for  the  parameters  Po  and  a  obtained 
by  linear  regression  for  the  curves  shoun  in  Fig.  13. 
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TABLE  5 


Uo/Umf        Po         a       Correlation 

Coefficient 


3.81 

0.238 

0.1097 

0.979 

■3   1*7 

0.137 

0.1023 

0.944 

Z.53 

0.040 

0.1181 

0.337 

2.32 

0.031 

0.1399 

0.8S7 

Pr  =  Po  exp-C-a  1> 

Results  of  linear  regression  analysis  for  particle 
loading  density  (grans/cn  )  vs  height  above  the  bed 
surface  ( cm) . 


The  parameter  Po  physically  represents  the  particle  loading 
density  which  uould  be  obtained  if  the  sample  uere  taken  at  the 
surface  of  the  bed.   This  is  not  necessarily  the  case  as  is 
indicated  by  the  computer  model  uhich  is  discussed  in  chapter  VI, 
but  is  oniy  a  parameter  descibing  the  particle  loading 
distribution  in  the  region  of  the  data  obtained.   The  dashed  lines 
belou  freeboard  heights  of  4  cm  indicate  the  region  in  question. 
The  parameter  a  is  a  characteristic  length  of  decay  for  the 
particle  flux. 

A  correlation  betueen  the  values  for  Po  in  table  5  and  Uc/Umf 
is  shown  in  Fig.   Z0.   This  plot  shows  that  Po  is  closely  related 
with  Uo/Umf.   Po  vanes  uith  (Uo/Umf  -  1)  approximately  to  the  2.3 
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1  )  Showing  Strong 


>o  on  Uo. 


pouer,  uhere  as  Ualsh  et  al  [10]  determined  the  coefficient  to  be 
about  2.1.   Uen  and  Chen  t5]  reported  a  correlation  for  particle  flux 
uhich  is  proportional  to  bubble  diameter  and  (Uo/Umf  -  1 >  to  the 
5/2  pouer.   These  correlations  are  listed  in  table  S. 


TABLE  S 

_  Present  uorlc    B  .»  0.13  <Uo/Umf  -  1)     kg/  m   sec 
Uen  and  Chen    B  =  1.34  (04)  ( Uo-Umf )     kg/  m   sec 
Ualsh  et  al     B  =  18  ( Uo  -  Umf )    kg/  m   sec 


Comparison  of  least  square  fit  relations  for  Po 
as  functions  of  Uo  and  Unf. 


The  differences  in  these  correlations  are  due  to  the 
different  bed  configurations  in  uhich  the  data  uas  taken  and  the 
measurement  technique  used.   The  present  uork  utilized  a  bed  uith 
a  relatively  closely  spaced  tube  configuration  and  steel  grit 
(S.G.   8.1,  median  size  230  microns)  for  the  bed  mass.   Both 
ascending  and  descending  particle  fluxes  were  captured  in  the 
sample.   The  uork  of  Ualsh  et  al  [10],  used  a  bed  uith  tuo  (2) 
uidely  spaced  horizontal  serpentine  tubes  and  Ottaua  sand  (S.G. 
2.G,  median  size  755  microns).   Also,  only  descending  particle 
flux  uas  used  in  determining  their  relations.   The  correlations  of 
Uen  and  Chen  [5]  are  a  result  of  studies  conducted  on  previous 
research  using  cylindrical  column  beds  and  lou  mass  bed  materials 

77 


(S.G.   0.8  -  2.S).   The  data  for  these  analysis  is  based  mainly  on 
pressure  measurements. 

It  has  been  observed  by  Leuis  et  al  [31  and  (Jen  and  Chen  £51 
that  a  is  not  a  strong  function  of  Uo.   Both  of  these  studies 
recommend  that  the  characteristic  particle  decay  length,  1/a, 
could  be  approximated  by  an  expression  of  the  form: 


1/a  =  C  Uo  (11) 


Table  7  is  a  list  of  correlations  obtained  by  other  studies 
and  in  the  present  uork .   The  study  by  Leuis  et  al  [31  uas 
conducted  uith  75  micron  glass  spheres  in  a  cylindrical  bed.   Fig. 
21  shows  the  relationship  between  1/a  and  Uo  in  the  present  uork. 


TABLE  7 

Present  uork     1/a  =  (0.19  +  0.03  s)Uo  m 

Leuis  et  al      1/a  =  (1.42  +  0.14  s)Uo  m 

Uen  and  Chen     1/a  =  (0.2S  +  0.09  s)Uo  m 

Ualsh  et  al      1/a  =  (0.32  +  0.05  s)Uo  m 


Comparison  of  least  square  fit  relations  for 
1/a  as  functions  of  Uo  (m/s). 
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Particle  Size  Distribution 

Tuo  methods  for  determining  particle  size  distribution  uere 
used.   The  sieve  method,  described  in  Appendix  E,  uas  used  for  bed 
material  analysis  only,  due  to  the  small  sample  sizes  obtained  from 
the  trap.   To  analyze  the  small  trap  samples,  an  image  analyzer 
uas  utilized.   The  procedure  used  uith  the  image  analyzer  is 
described  in  the  experimental  procedure  chapter. 

The  sieve  data  representing  bed  mass  as  a  function  of 
particle  size,  is  listed  in  Appendix  E.   Fig.   22a  shous  the  bed 
mass  data  converted  to  particle  number  as  a  function  of  particle 
diameter  uhere  the  interval  of  particle  diameter  between 
successive  measurements  is  10  microns.   These  plots  are  faired 
from  the  data  listed  in  Appendicies  E  and  J.   Fig.   22b  shous  the 
bed  mass  distribution.   The  conversion  from  a  mass  distribution  to 
a  particle  number  distribution  uas  calculated  assuming  that  each 
particle  uas  spherical  in  shape.   The  volume  corresponding  to  a 
given  particle  diameter  uas  multiplied  by  the  particle  density  to 
get  a  unit  particle  mass.   The  mass  fraction  of  the  sieve  analysis 
corresponding  to  the  specified  particle  diameter  uas  then  divided 
by  the  unit  particle  mass  to  obtain  the  represent i tive  particle 
number.   The  overall  resulting  particle  number  distribution  curve 
uas  then  normalized  uith  rescect  to  a  maximum  value  of  100. 

The  particle  size  distribution  data  obtained  from  the  image 
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Sieve   Analysis 

"i 


Image   Analyzer    Analysis 


\&  100  iUU  3Qid  4  00  500  b00 

Particle   Diameter    (microns) 

a)   Particle  Distributions 


130 


230 


1330 


Fig.  22 


Particle  Diameter  (microns) 

b)  Mass  Distribution 

Particle  Size  vs  Mass  distribution 
and  Particle  Number  of  Bed  Material. 
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analyzer  is  listed  in  Appendix  J.   The  image  analyzer  uas  used  to 
analyze  the  bed  nass  and  sample  data  from  the  Uo/Umf  =  3.81  data 
set.   The  bar  graphs  and  the  cumulative  percentage  plots  shou  the 
number  of  particles  vieued  by  the  image  analyzer  plotted  as  a 
function  of  normalized  particle  diameter. 

Fig.   22a  shows  a  plot  of  the  particle  number  distribution  as 
a  function  of  particle  diameter  as  determined  by  the  image 
analyzer.   The  data  is  normalized  with  respect  to  100  and  is 
compared  uith  the  data  as  determined  by  the  sieve  analysis.   The 
descrepancy  between  the  tuo  plots  can  be  explained  by  the  methods 
used  to  determine  the  respective  data.   The  analyzer  first 
determines  the  cross  sectional  area  of  the  vieued  particle.   n 
circle,  having  the  same  cross  sectional  area,  is  then  calculated. 
This  results  in  averaging  the  smaller  minimum  diameter  uith  the 
larger  maximum  diameter  of  all  particles  uith  long  cylindrical  or 
ellipsoidal  shapes.   In  the  sieve  however,  a  large  portion  of  the 
particles  will  pass  through  the  sieve  screen  by  means  of  the  small 
cross  sectional  area  presented  by  their  longitudinal  direction. 
To  be  consistant  in  the  following  sections,  the  image  analyzer 
data  will  be  used  to  correlate  all  particle  density  distributions. 

The  particle  size  distribution  as  a  function  of  height  above 
the  bed  surface  uas  evaluated  using  the  data  collected  at  Uo/Umf  = 
3.81.   The  image  anaiyser  was  used  to  determine  the  distribution 
by  viewing  three  (3)  random  samples  from  each  of  the  six  (G) 
height  positions.   Table  8  lists  the  various  statistical  values 


for  the  particle  size  distribution  data  obtained  (Appendix  J! 


TABLE  8 


eedoard 

Average 

Median 

Mode 

Di 

ameter  at 

leight 

Diameter 

Diameter 

Diameter 

Cu 
Di 

502  of 
nulative 
Number 
stribution 

(en) 

( microns ) 

( microns ) 

( microns ) 

( 

microns ) 

3.8 

199 

19Z 

14Z 

159 

7.5 

219 

Z14 

187 

203 

1Z.3 

204 

189 

1SZ 

157 

17.8 

179 

1S4 

154 

179 

22.2 

197 

179 

15Z 

177 

31.1 

Z1Z 

197 

194 

190 

Staistical  values  for  particle  number  distribution 
as  a  function  of  freeboard  height.   A  complete 
listing  of  the  data  is  given  in  Appendix  J. 


Fig.   23  shous  these  values  plotted  against  the  bed  height  at 
which  they  uere  taken.   The  excursion  of  points  at  the  7.5  cm 
height  is  assumed  to  be  due  to  analysis  error.   The  linear 
regression  lines  for  the  average  and  median  values  shou  that  they 
are   ueak  functions  of  collection  height.   For  the  particle 
diameter  values  representing  the  mode  and  the  502  point  on  the 
cumulative  percentage  plot,  the  linear  regression  lines  shou  a 
stronger  dependence  on  collection  height. 
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Oscilloscope  Trace  Analysis 

Figs.   24  thru  27  shou  examples  of  the  oscilloscope  traces 
given  in  Appendix  K.   Figs.   24  and  25  represent  the  traces 
obtained  during  lou  Uo  coditions.   Under  these  conditions,  the 
normal  bubble  probe  output  is  high  (  +  5  volts  dc).   Uhen  a  bubble 
errupts,  particles  are  throun  from  the  surface  of  the  bed  and 
eclipse  the  light  path  at  the  tip  of  the  probe  (trace  goes  to 
zero).   Figs.   2S  and  27  represent  the  traces  obtained  at  higher 
Uo  conditions  uhen  the  bubble  probe  is  normally  eclipsed  (bubble 
probe  output  is  lou,  zero)  by  particles.   The  prese.nce  of  a  bubble 
is  detected  by  a  high  output  from  the  bubble  probe  due  to  the 
bubble  creating  a  void  through  uhich  the  light  bean  can  pass.   The 
anemometer  output  just  above  the  bubble  probe  is  used  to  determine 
the  velocity  of  the  gas  jet  leaving  the  erupting  bubble.   As  seen 
in  the  figures,  there  is  not  aluays  a  bubble  associated  uith  a  gas 
jet  and  vice  versa.   This  is  due  to  the  anemometer  probe  being 
slightly  offset  from  the  bubble  probe  as  described  in  the 
experimental  procedure  section. 

The  gas  jet  in  Fig.   24  is  delayed  5  ms  and  the  gas  jet  in 
Fig.   25  is  delayed  35  ms  from  the  point  uhere  the  bubble  traces 
begin  their  excursion  to  the  zero  (eclipsed)  condition.   Since  the 
seperation  betueen  the  bubble  probe  and  the  anemometer  is  1  cm, 
this  results  in  an  estimated  jet  velocity  of  200  cm/s  (Fig.   Z4) 
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Bed  Surface  Eclipses  Pro 


Ret i vat  ion  of  Solenoid 


Fig.  24    Oscilloscope  Trace  of  Bubble  Probe, 
Anemometer  Probe,  and  Solenoid 
Actuation  at  Lou  Uo. 


Bed  Surface  Eclipses  Prob 


Activation  of  Solenoid 


Fig.  Z5     Oscilloscope  Trace  of  Bubble  Probe, 
Anemometer  Probe,  and  Solenoid 
Actuation  at  Lou  Uo. 


8G 


Activation  of  Solenoid 


Fig.  ZS    Oscilloscope  Trace  of  Bubble  Probe, 
Anemometer  Probe,  and  Solenoid 
Actuation  at  Higher  Uo. 


Activation  of  Solenoid 


Fig.  11 


Oscilloscope  Trace  of  Bubole  Probe, 
Anemometer  Probe,  and  Solenoid 
Actuation  at  Higher  Uo. 
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and  28  cn/s  (Fig.   25).   Using  the  calibration  curve  in  Appendix 
M,  the  velocities  are  determined  to  be  518  cn/s  for  the  jet  in 
Fig.   24  and  792  cn/s  for  the  jet  in  Fig.  25.   The  large 
difference  in  these  values  could  be  due  to  a  bulge  in  the  bed 
surface  from  the  bubble  underneath  eclipsing  the  bubble  probe. 
If  this  is  the  case,  there  mould  be  a  delay  betueen  the  tine  at 
uhich  the  bubble  probe  uas  eclipsed  and  the  eruption  of  the  bubble 
from  the  bed  surface. 

An  analysis  of  Figs.  2S  and  27  results  in  the  sane 
discrepancies  betueen  velocities  calculated  from  delay  tines  and 
measured  by  the  anenometer  probe.   Fig.  27  shows  another  phenomena 
uhich  occurs  quite  frequently.   The  gas  jet  appears  before  the 
bubble  is  detected.   There  are  several  possibilities  uhich  can 
explain  these  occurances.   First,  the  bubble  occurs  off  center 
from  the  bubble  probe.   Under  this  condition,  the  bubble  may  erupt 
and  initiate  a  jet  uhich  is  registered  by  the  anemometer.   The 
bubble  then  continues  to  rise  and  the  side  of  the  bubble  is 
registered  by  the  bubble  probe.   This  explanation  can  be  altered 
to  include  bubbles  coalescing  belou  the  surface.   Under  this 
condition,  a  bubble  may  be  pulled  into  the  vortex  of  an  already 
erupting  bubole  and  it  eclipses  the  bubble  probe. 

An  analysis  of  the  jet  velocities  and  durations  for  all  the 

oscilloscace  traces  in  Appendix  K,  resulted  in  the  data  listed  in 
table  9.   In  later  calculations,  the  values  to  be  used  for  the  jet 
velocity  and  duration  uill  be  B09 . S  cn/s  (20  ft/s)  and  20  ms 
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respectively. 


TABLE  9 


Average  Stnd  Dev 

Jet  Velocity      S39     cm/s  533    cm/s 

2Z.5   ft/s  18.    ft/s 

Jet  duration      21.25   ms  13.22  ms 


Average  and  standard  deviation  of  jet  velocity 
determined  from  oscilloscope  traces  in 
Appendix  K. 


Sample  Ueioht  Versus  Bed  Activity  Correlation 

The  main  purpose  of  inserting  the  bubble  and  anemometer  probe 
beneath  the  sampling  apparatus  uas  to  determine  uhether  a 
correlation  exists  betueen  the  sample  ueight  collected  and  the 
presence  of  bubble  eruptions  and  gas  jets.   To  evaluate  the 
photographs  taken  of  the  oscilloscope  traces  (Appendix  K),  the 
following  information  uas  required^ 

A.  The  average  velocity  of  the  particles  as  they 
travel  from  the  bed  surface  to  the  trap. 

B.  The  distance  the  particle  must  travel  to  reach 

the  trap. 

C.  The  closure  time  of  the  trap  relative  to  the 
bubble  eruptions  and  -gas  jets. 
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The  average  particle  velocity  uas  calculated  by  determining 
the  average  particle  size  and  using  the  height  output  from  the 
computer  trajectory  model.   The  average  particle  size  as 
determined  by  the  image  analyzer  uas  180  micons.   To  determine  the 
average  velocity  using  the  model,  the  following  initial  conditions 
were  input  to  the  program: 

Suoerficial  velocity  ( Uo )  =  57.3  cm/s   (1.9  ft/s) 

This  uas  the  actual  velocity 
measured  during  sampling. 


Jet  velocity  (  Uj  )  =  S09 . S  cm/s   (20  ft/s) 

This  value  uas  determined  from 
the  anemometer  data. 


Jet  duration  (t  jet)  =  20  ms 


This  value  uas  determined  from 
the  oscilloscope  traces. 


The  height  attained  by  a  180  micron  particle,  as  determined 
by  the  mocel  uith  the  above  conditions  is  44.7  cm  in  0.295 
seconds.   This  results  in  an  average  velocity  of  44.7/0.295  = 
151.5  cm/s  (4.97  ft/s). 

The  distance  a  particle  must  travel  to  reach  the  center  of 
the  trap  from  the  bee  surface  is  obtained  from  the  data  listed  in 
P.ppendix  I.   Since  trao  height  is  measured  from  the  bed  surface  to 
the  bottom  of  the  trap,  4  cm  (0.13  ft)  must  be  added  to  the  trap 
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heights  to  obtain  the  distance  to  the  center  of  the  trap.   The 
total  tine  required  for  the  particles  to  leave  the  bed  surface  and 
arrive  at  the  center  of  the  trap  can  nou  be  determined.   Table  10 
shous  the  transit  tine  for  each  set  of  data  for  uhich  the 
oscilloscope  traces  uere  photographed  (Appendix  K). 


TABLE  10 


Data      Trap       Trap  Total  Time 

set      Height    Height  Transit  Before 

Number     3ottom     Center  Time  Actuation 

(cm)       (cm)  ( ms )  ( ms ) 

27.13  179  136 

30.03  128  155 

26.23  173  130 

28.76  190  147 

17.65  IIS  73 


List  of  transit  times  for  particles  traveling  from 
the  bed  surface  to  the  center  of  the  trap.   The  time 
prior  to  actuation  of  the  sample  trap  is  also  shown 
(Total  time  -  42 . S  ms ) . 


56-65 

24.13 

66-75 

26.03 

76-85 

22.23 

86-95 

24.76 

36-105 

13.55 

The  total  closure  time  of  the  sampling  apparatus  (closure  of 
solenoid  pouer  suoply  suitch  SI  until  trap  is  shut)  is  determined 
in  Appendix  D  to  be  42.6  ms .   This  value  is  subtracted  from  the 
particle  transit  time  to  obtain  the  time  before  actuation  in  uhich 
particles  leaving  the  bed  surface  will  be  caught  in  the  trap. 

To  determine  whether  or  not  a  correlation  exists  between  the 
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ueight  of  a  sample  and  the  amount  of  bed  activity  present  in  the 
bed,  a  weighted  analysis  uas  used.   Each  sample  uas  evaluated 
three  uays.   First,  if  the  sample  ueight  uas  less  than  the  average 
sample  ueight,  a  ueight  factor  (U)  of  -1  uas  assigned  to  it.   If 
the  sample  ueight  uas  greater  than  the  average,  a  ueight  factor  of 
+1  uas  assigned.   Second,  the  oscilloscope  trace  uas  analyzed  at 
the  point  corresponding  to  the  "time  before  actuation"  listed  in 
Table  10.   This  position  is  the  time  before  the  solenoid  pouer 
supply  trace  is  present  (Fig.   Z8>.   If  the  trace  shoued  signs  of 
a  bubble  eruption  or  a  gas  jet  at  this  point,  the  position  factor 
(p)  uas  assigned  a  value  of  +1.   If  no  activity  uas  present,  the 
position  factor  uas  assigned  a  value  of  -1.   Third,  each  trace  uas 
analyzed  again  at  the  "before  actuation  time"  but  bed  activity 
uithin  a  +/-  25  ms  region  uas  counted.   If  there  uas  bed  activity 
uithin  this  region,  the  region  factor  (r)  uas  assigned  the  value 
of  +1.   If  no  activity  uas  present,  the  region  factor  uas  assigned 
the  value  of  -1.   The  follouing  equation  uas  used  to  assign  a 
correlation  factor  to  each  sample'- 


Q  =  U  p  +  U  p  (10) 

uhere: 

Q  =  Correlation  factor 

U  =  +1  if  sample  is  >  average  of  data  set 
-1  if  sample  is  (  average  of  data  set 

p  =  -1  if  no  activity  is  present  in  bed 

+1  if  bubble  or  gas  jet  activity  present 

07 


Time  ( ns ) 
200  300 


500 


Activation  of  Solenoid 


Trap  Closes 
at  4Z  ms 


Fig.  28    Typical  Oscilloscope  Trace  During 
Sampling  Procedure. 
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i"  :  -1  if  no  activity  is  present  +/-  ZS  ns 
+1  if  activity  is  present  +/-  ZS  ns 


Table  11  shows  the  correlation  factor  for  each  sample.   The 
resulting  average  correlation  factor  for  all  samples  is  0.3Z. 
This  value  suggests  that  no  correlation  can  be  made  between  the 
sample  ueight  and  bed  activity.   There  are  several  explanations 
for  this.   Babbles  not  detected  by  the  bubble  probe  or  debris  from 
previous  bubbles  returning  to  the  bed  surface  can  influence  the 
sample  ueight  by  increasing  the  amount  of  particles  collected. 
Also,  bubbles  of  smaller  size  have  a  smaller  velocity  and  uould 
therefore  produce  a  particle  stream  uhich  is  either  slower  or  just 
delayed  in  leaving  the  bed  surface.   This  uould  result  in  louer 
sample  weights  than  expected. 
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Sample 
Nunber 


TABLE  11 


55 
57 
58 
53 
50 
52 
S3 
54 

56 
67 
63 
78 

77 
73 
83 
84 

86 
88 
83 
30 
32 
33 
34 
35 


2 

-2 

2 

2 

-2 

-2 

2 

2 

0 

2 

-2 

-2 

2 

-2 
2 
2 

0 
2 
•2 
0 
0 
2 
0 
2 


36 

37 

38 

33 

1C0 

101 

102 


Average  =    0. 32 


List  of  samples  and  their  correlation  parameters, 
resulting  average  value  for  Q  indicates  that  no 
correlation  can  be  made  from  the  data  obtained  to 
indicate  by  sample  weight  whether  or  not  any  bed 
activity  occured  belou  the  sample  trap. 


The 


35 


CHAPTER   VI 

TRAJECTORY  MODR.  RESULTS  AND  DISCWSSTQN 
Selection  of  Baseline  Parameters 

The  input  to  the  node!  consists  of  the  following  five  (5) 
parameters : 

1)  Superficial  bed  velocity  (Uo) 

2)  Initial  particle  velocity  (Upo) 

3)  Peak  gas  jet  velocity  ( U j ) 

4)  Gas  jet  duration  (tj) 

5)  Particle  distribution  of  the  bed  mass 

The  baseline  values  for  each  of  these  inputs  uas  determined 
to  be  as  close  to  the  actual  experimental  bed  conditions  as 
possible. 

Superficial  Bed  Velocity 

The  superficial  bed  velocity  uas  determined  directly  from  the 
experinental  data.   For  the  bed  conditions  discussed  in  this 
section,  the  superficial  velocity  used  is  the  same  as  the  highest 
velocity  condition  under  uhich  the  particle  sampler  uas  used. 
The  superficial  velocity  uithin  the  bed  for  these  samples  uas 
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calculated  to  be  57. 9  cn/s  (1.9  ft/s). 
Initial  Particle  Velocity 

To  determine  the  initial  velocity  of  the  particle,  it  uas 
assumed  that  the  particle  uas  located  at  the  nose  of  a  bubble  and 
would  therefore  have  the  bubbles  velocity.   To  determine  the 
bubble  velocity,  Kunii  and  Levenspiel  ill    give  the  following 
equation^ 


Ub  =  0.711  [g  QbJ  1+  Uo  -  Umf  (11) 


where: 


Ub  =  Bubble  velocity 

Uo  =  Superficial  bed  velocity 

Umf  =  Minimum  fluidization  velocity 

Db  =  Bubble  diameter 

-  -  -  rg - =  graviatational  acceleration 


Observations  of  the  bed  material  during  the  sampling 
operation  suggests  that  the  average  bubble  diameter  present  in  the 
bed  is  approximately  B  cm  (Z.4  in).   For  Uo  =  57.9  cn/s,  Umf  = 
15. Z  cm/s  and  Db  =  Bern,  the  bubble  velocity  is  calculated  to  be 
97.2  cm/s  (3.19  ft/s).   This  value  was  used  as  the  initial 
particle  velocity  for  the  base  line  data. 
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Peak  Gas  Jet  Velocity 

To  determine  an  average  peak  jet  velocity,  the  anemometer 
output  on  the  oscilloscope  traces  (Appendix  K>  uere  analyzed.   To 
determine  the  relationship  betueen  the  voltage  output  from  the 
anemometer  (uhich  is  displayed  on  the  oscilloscope)  and  the  actual 
gas  jet  velocity,  the  anemometer  uas  calibrated  in  a  uind  tunnel. 
The  resulting  calibration  curve  for  the  anemometer  output,  as 
displayed  on  the  oscilloscope  trace,  is  given  in  Appendix  M.   The 
average  peak  velocity  of  the  anemometer  traces  analyzed  in  chapter 
V  uas  determined  to  be  approximately  609. S  cn/s  ( Z0  ft/s). 

Sas  Jet  Duration 

To  determine  an  average  gas: jet  duration  time  to  use  as  a 
baseline,  the  anemometer  traces  in  Appendix  K  uere  used.   The  time 
duration  of  each  gas  jet  analyzed  is  determined  directly  from  the 
oscilloscope  trace.   The  average  duration  of  a  gas  jet  uas 
determined  to  be  approximately  Z0  ms  in  chapter  V. 

Particle  Distribution  of  the  Bed  Mass 

Tuo  particle  distributions  for  the  bed  mass  uere  available 
for  use  as  the  baseline  data.   Both  particle  distributions  uere 
determined  using  the  same  sample  but  analyzed  using  different 
analysis  techniques.   The  tuo  methods  used  are  the  sieve  method 
and  the  image  analyzer  method,  both  of  uhich  are  described  earlier 
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in  chapter  V.   The  tuo  distributions  are  shoun  in  Fig.   Z3. 
The  bed  distribution  determined  by  the  image  analyzer  uas  used  to 
be  consistent  uith  the  particle  distributions  measured  at  various 
freeboard  heights. 


Typical  Output  Using  Baseline  Parameters 

Fig.   30  shous  the  resulting  maximum  particle  height  obtained 
as  a  function  of  particle  diameter  for  the  baseline  conditions. 
The  general  shape  of  the  curve  is  determined  by  the  dominating 
force  acting  on  the  particle.   For  large  particles,  the  dominating 
force  at  these  air  velocities  is  gravity.   Therefore,  momentum 
(initial  particle  velocity)  is  the  controlling  factor  determining 
the  maximum  height  attained  by  the  particle.   As  the  particle  size 
decreases,  the  proportion  of  drag  force  to  gravitational  force 
becomes  larger.   For  the  smaller  particles,  the  drag  force  becomes 
the  dominate  force.   The  result  is  that  the  maximum  height 
attained  by  a  particle  continues  to  increase  as  particle  size 
decreases.   This  continues  until  the  value  of  Uo  approaches  the 
terminal  velocity  of  the  smaller  particles.   At  this  point,  the 
particle  is  totally  dominated  by  the  air  flow  in  the  bed.   This  is 
seen  in  the  small  decrease  in  height  attained  by  particles  less 
than  110  microns  as  a  result  of  the  particles  rapid  deacceleration 
tc  Uo  after  the  gas  jet  has  stocned  while  the  larger  particles 
continue  a  little  higher  due  to  momentum.   The  particle  height 
continues  to  decrease  for  a  short  time  until  Uo  becomes  equal  to 
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Faired  Data  From  Appendix  E  and  J.  Interval  of  10  microns. 
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or  greater  than  the  terminal  velocity  for  the  remaining  particles. 
These  particles  would  be  elutriated  from  the  system. 

Figs.   31  thru  3B  shou  the  individual  particle  densities  at 
different  heights  above  the  bed.   fill  of  the  heights  shou  an 
increase  in  particle  density  for  the  smaller  diameter  particles. 
This  is  the  result  of  the  smaller  particles  falling  at  their 
respective  terminal  velocities  uhich  is  slower  than  for  the  larger 
particles.   The  result  is  that  the  smaller  the  particle  velocity, 
the  longer  the  particle  exists  uithin  a  given  height  region  and 
thus,  the  larger  the  particle  density.   It  is  also  due  to  the 
larger  number  of  smaller  particles  present  uithin  the  system. 

Fig.   37  shows  the  particle  density  distribution  above  the 
bed.   The  two  curves  represent  the  effect  of  varying  the  range  of 
particle-  diameters  used  in  the  model.   The  two  ranges  are  80  -  570 
microns  and  80  -  1070  microns.   The  peak  of  the  80  -  570  micron 
distribution  is  at  about  18  cm  of  bed  height  whereas  the  other 
distribution  peaks  at  about  11  cm.   The  region  shown  on  the  graph 
below  these  heights  is  the  area  known  as  the  splash  zone.   The 
shape  of  this  part  of  the  curve  is  due  to  the  initial  acceleration 
of  the  particles  leaving  the  surface  of  the  bed  by  the  gas  jet. 
The  particles  then  begin  to  slow  down  due  to  drag  and  gravity 
resulting  in  an  increase  in  particle  density.   The  peak  on  this 
curve  coincides  with  the  maximum  height  attained  by  the  largest 
particles,  and  therefore,  it  can  be  assumed  that  the  change  in 
slope  is  due  to  the  loss  of  particles  as  they  return  to  the  bed. 
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Fig.   38  is  a  semi-ln  plot  of  the  sane  data  as  plotted  in  Fig. 
37.   The  sharp  drop  at  the  right  end  of  the  plot  is  due  to 
analyzing  particles  uith  a  minimum  diameter  of  80  microns.   Had 
smaller  particles  been  included,  the  slope  uould  have  approached 
zero  instead  of  infinity  uhich  uould  model  the  elutriation  of 
particles.   The  slope  of  the  line  to  the  right  of  the  peaks  shown 
in  Fig.   38,  uhich  are  located  at  11  and  18  cm  of  height  above  the 
bed  surface  are  the  same  and  are  approximately  -0.117  grams/cm  . 
To  decrease  computation  time,  the  80  -  570  micron  particle 
distribution  uas  used  in  the  following  sensitivity  analysis  as  the 
decay  slopes  are  equal. 

Model  Sensitivity  Analysis 

--~    An  analysis  of  the  change  in  bed  characteristics  due  to  a 
change  in  one  of  the  five  model  input  parameters  uas  conducted. 
The  parameters  changed  and  the  values  used  are  as  follous- 

Superficial  velocity  (Uo): 

Baseline  Uo:     57.9  cm/s  (1.9  ft/s) 

1.  30.5  cm/s  ( 1  ft/s) 

2.  91.4  cm/s  (3  ft/s) 

Initial  particle  velocity  ( Upo ) : 

Baseline  Upo:    119.2  cm/s  (3.19  ft/s) 

1.  61 .0  cm/s  (2  ft/s) 

2.  152.4  cm/s  (5  ft/s) 
10G 


3.   305. 0  cn/s  (10  ft/s) 

Gas  Jet  velocity  ( U j ) : 

8aseline  Uj :     B09 . S  cn/s  (20  ft/s> 

1.   305.0  cn/s  (10  ft/s) 
Z.   457. Z  cn/s  (15  ft/s) 

Gas  jet  duration  (tj): 

Baseline  t j -  Z0  ns 

1.   10  ns 
Z.   30  ns 

Particle  distribution  of  bed  mass: 

.  Baseline  distribution-'   Inage  Analyzer 

.--•.,  ...  ,    1.   Sieve  analysis  data 

The  bed  characteristics  uhich  are  evaluated  against  the 
baseline  characteristics  are: 

1.   The  naxinun  height  attained  by  each  particle  fron  80  to 
570  microns. 

Z.   The  density  of  each  particle  size  (80  -  570  microns)  at 

freeboard  heights  of: 

a.  4  en  ( 1 . 5  in) 

b.  8  en  ( 3. 1  in) 
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c.  1Z  en  (4.7  in) 

d.  18  cm  (7.1  in) 

e.  ZZ  cm  (8.7  in) 

f.  31  en  (1Z.Z  in) 

3.   The  particle  density  distribution  in  the  freeboard. 


Variation  of  Superficial  Velocity  (Uo) 

Fig.   39  shous  the  influence  of  Uo  on  the  maximum  height 
attained  by  the  particles.   As  uas  described  in  the  previous 
section,  the  larger  particles  are  dominated  by  momentum  and  not 
drag.   This  is  shoun  in  Fig.   39  by  the  small  change  in  maximum 
height  attained  by  the  large  particles  due  to  a  change  in  Uo.   As 
.  the  particle  size  decreases,  the  drag  force  becomes  the  dominate 
force  and  the  effect  of  Uo  on  particle  height  increases.   As  seen 
in  Fig.   39,  a  change  in  Uo  produces  a  moderate  change  in  the 
maximum  height  attained  by  the  small  particles. 

Figs.   40  thru  45  shou  the  effect  of  Uo  on  individual 

particle  densities  at  different  heights  above  the  bed.   All  of  the 

heights  shou  an  increase  in  particle  density  for  the  smaller 

diameter  particles  as  Uo  is  increased.   This  is  the  result  of  the 

smaller  particles  falling  at  their  respective  relative  terminal 

velocities  uhich  is  slower  for  smaller  particles.   A  decrease  in 

the  maximum  particle  diameter  present  at  a  given  freeboard  height 

as  Uo  is  decreased  can  be  observed  in  Figs.   40  thru  45.   This 
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Relative  Particle  Number  vs  Particle 
Diameter  as  a  Function  of  Uo. 
Freeboard  Height  of  22  cm. 
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decrease  is  due  to  the  decrease  in  lift  given  to  the  ascending 
particle  as  Uo  is  decreased.   As  a  result,  all  particles  achieve  a 
louer  maximum  height  uhen  Uo  is  decreased. 

Fig.   4S  shows  the  effect  of  Uo  on  the  particle  density 
distribution  above  the  bed.   There  is  a  slight  increase  in  overall 
particle  density  uith  an  increase  in  Uo.   fts  expected,  there  is 
also  an  increase  in  the  maximum  height  attained  by  the  particles 
uith  an  increase  in  Uo.   Fig.   47  is  a  semi-ln  plot  of  the  same 
data  plotted  in  Fig.   46.   Table  12  lists  the  slopes  of  the  lines 
to  the  right  of  the  peak  shoun  in  Fig.   47  uhich  is  located  at 
about  18  cm  of  height  above  the  bed  surface.   This  table  shous  an 
increase  in  the  particle  distribution  slope  as  Uo  is  decreased. 


TABLE 

12 

Uo 

Slope 

Intercept 

cm/s 

gms/cm 

gms/cm 

30.5 

-0.136 

7.036 

57.9 

-0.117 

6.758 

91.4 

-0.098 

6.473 

Effect  of  Uo  on  the  slope  of  the 
particle  density  distribution  as  a 
function  of  height  for  the 
distributions  shoun  in  Fig.  47. 
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Variation  of  Initial  Particle  Velocity  ( Upo ) 

Fig.   48  shous  the  relation  between  initial  particle  velocity 
(Upo)  and  the  maximum  height  attained  by  particles  of  different 
size.   This  figure  also  shous  that  the  larger  particles  are  mainly 
momentum  dependent  uhile  the  smaller  particles  are  drag  dependent. 
If  Upo  uere  increased  further,  the  larger  particles  uould  continue 
to  increase  their  maximum  height.   The  smaller  particles  uould 
approach  a  height  limit  uhich  is  dependent  upon  the  superficial 
velocity  in  the  bed. 

Figs.   49  thru  54  shou  the  effects  of  increasing  Upo  on  the 
individual  particle  densities  for  increasing  heights.   In  general, 
the  individual  particle  distributions  undergo  the  same  relative 
changes  from  the  bed  surface  to  the  maximum  height  position.   The 
difference  being  the  height  above  the  bed  surface  at  uhich  the 
particular  distribution  is  present. 

Fig.   55  shous  the  effect  of  varying  Upo  on  the  particle 
density  distribution  in  the  freeboard.   As  uas  noted  previously, 
the  peak  density  occurs  at    the  point  uhere  the  largest  particles 

attain  their  maximum  height  above  the  bed.   As  Upo  is  increased, 
the  curve  to  the  right  of  the  peak  becomes  shorter  and  steeper. 
This  trend  uill  continue  until  the  the  effect  of  the  smaller 
particles  returning  to  the  bed  cause  the  slope  to  decrease.   Fig. 
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SB  is  a  semi-ln  plot  of  the  data  in  Fig.   55.   Table  13  confirms 
that  the  slope  is  becoming  steeper  as  Upo  is  increased  although 
the  influence  of  initial  particle  velocity  on  the  slope  of  the 
density  vs  height  curve  is  moderate. 


TABLE  13 


Upo        Slope     Intercept 
cm/s       gms/cm     gms/cm 

SI. 0  -0.118  6.058 

97.2  -0.117  S.758 

15Z.4  -0.125  8.500 

304.8  -0.200  18.20 

Effect  of  Upo  on  the  slope  of  the 
particle  density  distribution  as  a 
-function  of  height  for  the 
distributions  shown  in  Fig.  56. 


Variation  of  Jet  Velocity  <  Uj  ) 

Fig.   57  shows  the  effect  on  maximum  particle  height  uhen  Uj 
is  varied.   As  in  the  two  previous  analysis',  the  influence  of 
momentum  and  drag  on  the  different  particle  sizes  is  apparent.   A 
decrease  in  jet  velocity  from  the  baseline  value  results  in  a 
large  drop  in  small  particle  height  and  a  small  change  in  large 
particle  height.   This  suggests  that  jet  velocity  has  a  strong 
effect  on  the  aiscersion,  or  seperation  of  particles  of  different 
sizes  at  increasing  freeboard  heights.   Increasing  Uj  results  in  a 
stretching  effect  of  the  particle  size  distributions  above  the 
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bed.   One  interesting  point  to  note  is  the  loss  of  the  rounded 
edge  at  the  heights  corresponding  to  the  smaller  particles.   This 
shous  the  extremely  large  effect  that  drag  plays  uihen  the  air 
velocity  is  very  close  to  the  terminal  velocity  of  the  particle. 

Figs.  58  thru  S3  shou  the  effect  of  varying  Uj  on  the 
individual  partical  densities  at  varying  heights  above  the  bed. 
Again,  as  in  the  previous  parameter  analysis,  the  same 
distribution  shapes  can  be  seen  for  each  value  of  Uj ,  the  only 
difference  being  the  height  at  uhich  it  occurs. 

Fig.   64  shous  the  change  in  particle  density  distribution  as 
_Uj  J.s_varied^___F.ig-  _^65-  is_aj:semi-ln_plot  of  this  data  shouing 
that  as  Uj  is  decreased,  the  slope  of  the  distribution  changes 
substantially.   Table  14  also  shous  this  effect.' 


TABLE  14 


Uj         Slope      Intercept 
cm/s       gms/cm     gms/cm 

S09.6  -0.117  S.7S8 
457. 2  -0.207  7.380 
304.8      -0.400      7.800 

Effect  of  Uj  on  the  slope  of  the 
particle  density  distribution  as  a 
function  of  height  for  the 
distributions  snoun  in  Fig.  BS. 
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Variation  of  Jet  (Duration  ( t  j  ) 

Fig.   SB  is  a  graph  shouing  the  effect  of  varying  jet 
duration  on  the  maximum  heights  attained  by  the  particles.   It  is 
easily  seen  that  a  small  change  in  jet  duration  results  in  a 
substantial  change  in  particle  height.   This  effect  indicates  that 
the  distribution  present  in  an  actual  fluidized  bed  is  probably  a 
statistical  average  of  a  rapidly  fluctuating  particle  distribution 
uhich  is  controled  by  the  durations  of  the  jets  from  the 
neighboring  bubble  eruptions  as  well  as  local  eruptions. 


Figs.   G7  thru  72  shou  the  effects  on  the  individual  particle 
density  distributions,  at  different  freeboard  heights,  as  a 
function  of  jet  duration.   These  figures  also  shou  the  drastic 
change  in  density  distributions  caused  by  changes  in  jet  duration. 

Fig.   73  shous  the  particle  density  distributions  above  the 
bed  as  affected  by  jet  duration.   Fig.   74,  uhich  is  a  semi-ln 
plot  of  the  data,  shous  the  drastic  changes  in  the  decay  slopes. 
Table  15  also  shous  this  drastic  variation  in  slope. 
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TABLE  IS 


tj        Slope     Intercept 
cm/s       gms/cm     gms/cm 

10.0  -1.071  12.50 

20.0  -0.117  6.7S8 

30.0  -0.020  S.234 

Effect  of  tj  on  the  slope  of  the 
particle  density  distribution  as  a 
function  of  height  for  the 
distributions  shoun  in  Fig.  74. 


Variation  of  Particle  Distribution  in  Bed  (lass 


Fig.-  75, shows  the  two  bed  particle  distributions  used  in 


this  analysis.   The  determination  of  these  distributions  is 
descibed  earlier  in  this  uork . 

Figs.   76  thru  81  shou  the  effects  of  varying  the  bed 
particle  distribution  on  the  individual  particle  distributions  at 
different  heights  above  the  bed  surface.   Using  the  sieve  particle 
distribution  instead  of  the  image  analysis  distribution  produces  a 
shift  in  the  mean  particle  distribution  towards  the  smaller 
particles.   The  distributions  resulting  from  the  sieve  data  do  not 
shift  as  much  as  the  image  analyzer  data  when  the  freeboard  height 
is  increased.   The  small  shift  in  riean  diameter  exhibited  by  these 
plots  is  consistant  uith  the  distributions  listed  in  Appendix  J. 

Fig.   8Z  shous  the  particle  density  distributions  in  the 
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freeboard  for  the  tuo  bed  mass  conditions.   The  effect  of  the 
sieve  distribution  (which  has  a  much  lower  average  particle 
diameter)  is  to  decrease  the  initial  rate  at  which  particle  mass 
returns  to  the  bed.   This  is  due  to  the  smaller  particles  having  a 
much  smaller  terminal  velocity.   This  will  not  only  raise  the  bulk 
of  the  particles  to  a  higher  height,  but  uill  also  increase  the 
time  required  for  the  particles  to  return  to  the  bed. 


Comparison  of  Model  with  Experimental  Results 

As  discussed  in  the  previous  section,  the  baseline  parameters 
were  selected  as  the  experimental  conditions  present  in  the  bed 


uhen  sampling  at  the  highest  Uo  setting.   Table  IS  lists  these 
baseline  parameters  again  for  review. 


TABLE  IS 


Uo   =  57. 9  cm/s  (1.9  ft/s) 
Upo  =  97.2  cm/s  (3.19  ft/s) 
Uj   =  609. S  cm/s  (20  ft/s) 
tj   =  20  ms 

Baseline  parameters  used  in  computer  model 


Fig.   93  shows  the  particle  density  distribution  predicted  by 
the  model  using  the  baseline  conditions.   Fig.   34  is  a  semi-log 
plot  of  this  same  data.   A  comparison  of  the  slope  of  the 
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distribution  line  after  the  peak  value  in  Fig.   84  with  the  slope 
obtained  from  the  expe^i mental  data  shows  very  good  agreement 
between  the  tuo.   This  data  is  also  given  in  Table  17  along  uith 
the  data  obtained  for  the  other  experimental  data  sets  and  is 
shoun  plotted  in  Fig.   85.   The  slopes  for  the  model  data  show  a 
steady  increase  in  slope  as  Uo/Umf  is  decreased.   This  trend  is 
also  followed  by  the  experimental  data  but  is  not  as  smooth.   The 
slope  of  the  experimental  data  is  also  changing  faster  than  the 
slope  of  the  model  data.   One  reason  for  this  is  that  all  the  model 
calculations  uere  performed  using  the  same  value  for  jet  velocity 
and  duration  of  609 .  S  cm/s  and  20  ms  resc-ectively .   The  magnitudes 
of  these  values  should  decrease  uith  decreaseing  Uo/Umf,  but  the 
values  to  be  used  for  these  different  conditions  could  not  be 
determined. 


TABLE  17 

Uo/Umf       Experimental  Model  Diff 

Slope  Slope  j 

gms/cm  gms/cm 

3.81          -0.1097  -0.117  S.G 

3.17          -0.10Z9  -0.1Z3  19. S 

Z.S9          -0.1181  -0.1Z9  9.Z 

Z.3Z          -0.1399  -0.13Z  5.S 


Comparison  of  slooes  for  the  particle  density 
distributions  acove  the  bed  as  derived  from 
the  experimental  cata  and  the  computer  model. 
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The  slope  shoun  in  Fig.   84  for  heights  between  approximately 
4  and  18  en  is  totally  different  from  the  slopes  listed  in  table 
17.   This  is  also  the  louer  range  in  uhich  the  experimental  data 
uas  collected.   As  uas  shoun  earlier,  when  the  particle 
distribution  analyzed  included  particles  from  80  -  1070  microns, 
the  peak  uas  located  at  11  cm.   This  is  the  height  uhere  the 
largest  particles  attain  their  maximum  height  and  begin  to  fall 
back  to  the  bed.   It  does  not  seem  likely  that  including  the  feu 
particles  present  above  1070  microns  uould  totally  account  for 
this  discrepancy.   Other  factors  uhich  the  model  has  not  accounted 
for  uhich  may  explain  these  discrepancies  is  the  effect  of  varying 
jet  velocities  and  durations  uhich  uere  shoun  in  chapter  V  to  vary 
greatly,  and  the  effect  of  particles  uhose  velocity  vector  is  not 
perpindicular  to  the  bed  surface. 


I  O.A. 


CHAPTER   VII 


CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 

Several  important  conclusions  can  be  made  concerning  the 
particle  trajectory  computer  model  and  the  experimental  data. 


1.   ft  particle  trajectory  model,  even  as  simple  as  the  one 
discussed  in  this  uork,  produces  results  uhich  closely  predict 


several  aspects  of  particle  activity  uithin  a  fluidized  bed.   The 
more  important  predictions  include  the  particle  density 
distribution  in  the  freeboard  and  the  height  distribution  of 
particles.   The  results  obtained  uith  the  current  model  indicate 
that  further  uork  on  improving  the  particle  trajectory  model  to 
include  particle-particle  interactions  and  paticle  velocities  not 
perpendicular  to  the  bed  surface  is  very  desirable. 

Z.   The  sensitivity  analysis  indicates  that  the  jet  duration 

and  jet  velocity  are  critical  parameters  in  determining  particle 

loading  conditions  in  the  freeboard.   Since  these  sr~   constantly 

changing  from  one  bubble  eruption  to  the  next,  a  statistical 

distribution  will  be  required  to  accuarately  model  freeboard 

particle  activity.   The  developement  of  a  basic  model  to  predict 
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the  jet  duration  and  velocity  as  a  function  of  bubble  size,  Uo,  Umf , 
etc  is  therefore  needed. 

3.   Particle  distribution  analysis  of  experimental  samples  at 
increasing  heights  above  the  bed,  shou  the  presence  of  large 
particles.   Many  of  these  are  above  the  maximum  calculated 
trajectory  heights  for  these  particle  sizes.   This  indicates  that 
either  the  jet  velocity  or  jet  duration  carried  these  particles  to 
these  heights,  or,  particle-particle  collisions  are  present  in 
numbers  great  enough  to  be  important  in  the  analysis. 


-4.   The  sampling  apparatus  designed  and  built  for  this  study 
appears  to  operate  satisfactory.   The  data  obtained  correlates 
uith  uork  done  by  other  researchers  and  uith  the  computer  model. 
Houever,  additional  uork  is  needed  to  correlate  sample  size  and 
particle  size  distribution  uith  freeboard  height. 

Recommendations 

1.   Further  experimental  data  is  needed  at  higher  velocities 
and  higher  heights  above  the  bed  using  the  sampling  apparatus 
designed  in  this  uork.   This  uill  provide  additional  data  to 
evaluate  the  operation  of  the  sampling  apparatus  and  the  computer 
model.   niso,  a  scale  cacable  of  measuring  quantities  of  samples 
less  than  3.31  grams  ana  a  particle  removal  system  uhich  removes 
the  particles  uith  less  error  from  the  sample  trap  is  necessary 

for  accurate  uork . 
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2-   A  sampling  device  uhich  has  a  shorter  trap  height  should 
be  designed.   The  present  sample  trap  is  very  direction  oriented 
and  samples  particles  traveling  only  in  a  narrow  range  from  the 
vertical.   This  would  allou  more  accurate  uork  to  be  done  in  the 
splash  zone  uhere  particles  are  more  likely  to  be  traveling  in 
directions  other  than  vertical. 

3.   The  computer  model  needs  to  have  encorporated  in  it,  a 
statistical  distribution  model  for  jet  velocity  and  jet  duration. 
This  uill  allou  the  determination  of  the  effects  of  varying  jet 


velocity  and  duration  on  particle  distributions  in  the  freeboard. 

4.  A  correlation  for  jet  velocity  and  duration  as  functions 
of  bubble  diameter,  Uo,  and  Umf  should  be  determined  and  included 
in  the  model. 
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APPENDIX      A 

Moment    of   Inertia  Calculations   for  Paddles 

The  calculation  of   the   total   Moment   of    inertia  (I)    of   the 
paddles   can  be  broken  up   into   three   seperate   calculations.      First, 
the   moment   of  inertia.  ( II  >    of   the   tuo   aluminum   cylinders   used   to 
mount    the   paddle   arms   to   the   solenoid   shaft    is   calculated.      The 
second   calculation   (IZ),    accounts   for    the   moment    of    inertia  of    the 
paddles    themselves,    uhich   are   constructed   of    a  foam,    bassuood   and 
epoxy   laminate    (Fig.      ft-1).      Final ly,:    the   moment    of    inertia   (13)    of 


the   harduood   mounting   ends   en   the   paddles    is   calculated. 

The  following   equations   uere  used   to   calculate   the   mass 
moments    of    inertia- 


II    =    P       H    r* 


(fll) 


12    =      L    U   (L2    j   UZ)    (PI    tl    *   PI    t2) 

12  <A2> 


13    =   ZLUtP(L2    +    U2) 

1Z  (A3) 

Table  ft.l  is  a  listing  of  parameters  recuired  for  the  moment 
calculations. 
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TABLE  fl.l 


Component 

Aluminum  Cylinder 

Density  Al  [Pi- 
Radius  [r] 
Height  [H] 

Paddles 

Density  Foam  CPU 
Density  3assuood  IPZ1 
Thickness  Foam  [tl] 
Thickness  Bassuood  [t21 
Length  Foam  [LI] 
Length  Bassuood  [L2] 
Uidth  Foam  LU11 
Uidth  Bassuood  CUZ] 


MKS  units 


Z780 

kg/m 

0.0171 

m 

0.0381 

m 

English  units 


0.1004   lbf/in' 
0.S7S    in 
1.5     in 


50 

kg/m 

0.00181 

lbf/ 

in 

400 

kg/m 

0.0144 

Ibf/ 

in3 

0.010Z 

m 
m 

0.4 
0.0GZ5 

in 
in 

0.00159 

-  ■ 

0.584Z 

m 

Z3.0 

in 

0.584Z 

m 

Z3.0 

in 

0.0305 

m 

1.2 

in 

0.0305 

m 

l.Z 

in 

Harduood 

Oensity  Maple  [PI 
Thickness  Maple  [t] 
Length  Maple  [L] 
Uidth  Maple   [U] 


650  kg/m 

0.010Z  m 

0.043Z  m 

0.0348  m 


0.0Z35  lbf/in 

0.4  in 

1.7  in 

l.Z  in 


Listing  of  paddle  components  and  parameters 


15Z 


Equation  (A.l)  is  used  to  calculate  the  mass  moment  of 
inertia  for  a  right  circular  cylinder  rotating  about  its  Z  axis 
(Fig.   A-Z).   The  height  H,  accounts  for  both  the  upper  and  the 
lower  solenoid  shaft  cylinders.   Equations  (ft. 2)  and  (A. 3) 
determine  the  mass  moment  of  inertia  for  a  rectangular  prism 
rotating  about  its  X-Y"  centrcidal  axis  (Fig.   A-3).   The  length  L 
~in  equation  AZ  accounts  for  both  the  upper  and  louer  paddle 
lengths.   The  factor  of  Z  in  equation  (A. 3)  is  because  both  the 
upper  and  louer  hardwood  sections  ar&   identical  and  can  be 
combined.   By  inserting  the  respective  values  from  Table  A.l  into 
equations  (A.l  )T  (A.Z),  and  (A. 3),-  the  value  for  each  inertia 


component  can  be  calculated.  Adding  these  components,  the  total 
mass  moment  of  inertia  applied  to  the  solenoid  by  the  paddles  is 
determined.   These  results  are  listed  in  Table  A.Z. 
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Fig-    A-Z        Diagram   for    Moment    of    Inertia  Calculation 
Used   for    Cylinder    About    Z    Axis. 


Fig.    A- 3        Diagram    for    Moment    of    Inertia   Calculation 
Used    for    Rectangular   Pnsim   About    X    Axis. 

1S4 


TABLE  A. 2 
Moment  of  inertia  components^ 
MKS  units  English  units 

11  =  1.4Z2  (-05)  kg  m2      4.911  (-02)  lbf  in2 

12  =  5.823  C-04)  kg  m*      1.982  (00)   lbf  in2 

13  —2. 555  (-06)  kg  m2      S.919  (-03)  lbf  in 

Total  moment:     I  =11  +12  +13 
I  =  5.991  (-04)  kg  rT       2.037  (00)   lbf  in 


Calculated  values  for  mass  moment  of  inertia  of 
paddles  and  paddle  components.   II  is  inertia 
value  for  the  aluminum  cylinders.   12  is  the 
inertia  value  for  the  paddles  themselves  and  13 
is  for  the  harduood  end  pieces. 
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APPENDIX   B 


Error  Determination  of  Vacuum  Collection  System 


Tests  of  the  vacuum  sample  removal  system  indicated  that  some 
particles  remained  in- the  sample  trap.  These  particles  were 
located  in  the  corners  of  the  trap  where  the  equalizing  air  stream 
could  not  agitate  them  enough  to  move  them  into  the  vacuum  strean. 
As  a  result,  it  became  necessary  to  determine  to  what  extent  these 
remaining^  particles  affected  the  accuracy  of  the  sample  attained. 


The  procedure  used  involved  placing  samples  of  known  weight 
and  particle  size  distribution  (equivalent  to  the  bed  material) 
inside  the  sample  trap.   Sample  sizes  of  10,  IS,  20,  and  25  grams 
uere  used.   The  sample  particles  were  then  vacuumed  out  and  their 
weight  determined.   The  difference  in  weights  of  the  samples  were 
then  calculated  along  with  the  percentage  differences.   These 
results  are  listed  in  Table  8.1. 

An  analysis  of  the  results  listed  in  Table  B. 1  indicates  that 
the  average  difference  between  the  two  sample  weights  is  only  0.52 
I  while  the  maximum  difference  observed  was  3.93  %.   It  can 
therefore  be  concluded  that  the  particles  remaining  in  the  sample 
trap  do  not  significantly  affect  the  accuracy  of  the  sample 
attained. 
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TABLE  9.1 

InitiaL 

Sample 

Ueight 

Ueight 

SAMPLE 

Ueight 

Ueight 

Difference 

Dif f erenc 

No. 

grans 

grans 

grams 

Percent 

1 

10 

9.95 

0.05 

0.50 

2 

10 

9.92 

0.08 

0.80 

3 

10 

9.95 

_ 

0.50 

0.05 

4 

10 

IS 

9.98" 
15.03" 

0.02 
-0.03 

0.20 

S 

-0.20 

G 

IS 

15.02 

-0.02 

-0.13 

7 

IS 

14.86 

0.14 

0.93 

8 

IS 

14.87 

0.13 

0.87 

9 

IS 

14.97 

0.03 

0.20 

ia   . 

20 

19.93 

0.07 

0.35 

ll  .; . . 

20  .\ '. 

19.35  ~ 

0.14 

0.70 

12 

20   . 1 1 
20 

19.90  . 

0.10 

0.50 

13 

19.93 

0.07 

0.35 

14   . 

20 

19.94 

0.06 

0.30 

IS 

25 

24.83 

0.17 

0.58 

IS 

2S 

24.80 

0.20 

0.80 

17 

25 

24.78 

0.22 

0.88 

18 

Z5 

Z4.83 

0.17 

0.58 

13 

Z5 

Z4.77 

0.Z3 

0.93 

RV6= 

0.10 

0.5Z 

STNUV= 

0.08 

0.34 

Results  of  vacunn  sample  removal  test, 
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APPENDIX   C 


Sample  Trap  Closure  Tine  Test 


Due  to  the  high  velocity  of  the  particles  with  respect  to  the 


size  of  the  sample  trap,  it  uas  important  to  ensure  that  the 
sample  trap  uas  closed  in  a  very  short  period  of  time.   To  measure 
the  closure  time,  the  set-up  shoun  in  Fig.   C-l  uas  used. 


fin  electric  eye,  consisting  of  a  light  emitting  diode  and  a 
photo  transistor,  uas  used.   The  diode  uas  pouered  by  a  G  volt 
battery  and  the  photo  transistor  uas  uired  directly  to  an 
oscilloscope.   The  electric  eye  circuit  uas  then  placed  at  a 
distance  of  0.ZS7  m  (10.5  in)  from  the  pivot  point  of  the  paddles, 
fit  this  position,  the  closing  time  calculation  involves  only  a 
simple  proportion  relationship  betueen  the  paddle  and  the  sample 
trap  uidths.   The  eye  uas  also  positioned  as  close  to  the  trailing 
edge  of  the  fully  closed  paddle  as  possible  to  minimize  error. 

The  oscilloscope  sueep  uas  set  to  trigger  off  of  the 

initial  change  of  state  from  the  photo  transistor  uhen  the  paddle 

first  eclipsed  the  light  beam.   The  resultant  traces  uere  then 

analysed  to  determine  closure  time.   Time  t=0  uas  set  equal  to  the 

initiation  of  the  trace  on  the  oscilloscope.   Time  t=ti  uas 

defined  as  the  paint  at  uhich  the  trace  begins  its  excursion  back 
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Fig.  C-l    Diagram  of  Closure  Time  Determination  Set  Up. 
ftl 1  Dimensions  in  cm. 
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to  the  base  voltage  as  shoun  in  Fig.   C-2.   By  using  the  ratios 


Ui  =  U2 

tl    t2  (C.l) 


Uhere: 


■ 


Ul  =  uidth  of  the  paddle 
U2  =  uidth  of  sample  trap    ----- 
—  -—     —  tl  - --   tine  delay  measured  on  oscilloscope 
t2  =  closure  time  of  sample  trap 


the  closure  time  of  the  sample  trap  can  be  determined.   Several 
trials  uere  run  with  the  resulting  data  listed  in  Table  C-l.   All 
of -the  trials  uere  very  consistant .:  -The  resultant  determination 


for  the  average  closure  time  uas  1.44  ms. 


TA8LE  C-l 


Ul    = 

0. 

0305   n 

(1.2 

in) 

uz  = 

3. 

,0127   m 

(0.5 

in) 

Trial 

tl 

tz 

No. 

ins 

ms 

1 

3.4 

1.42 

2 

3.4 

1.42 

3 

3.55 

1.48 

4 

3.5 

1.4G 

5 

3.5 

1  .45 

flvg        3.47       1.44 


Samel e  trap  closure  data:  Ui  is  the  uidth  of 
the  paddle,  U2  is  the  uidth  of  the  sample  trap 
and  time  tl  is  eclipse  time  of  paddle  through 
light  beam.   Time  t2  is  closure  time  of  Sample 
trap  as  calculated  using  equation  (C.l). 
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APPENDIX   D 


Solenoid  Torque  and  Dynamic  Analysis 


The  goal  of  this  design  uas  to  shut  the  sample  trap  in  about 
1  millisecond  <  ms ) .   At  this  speed,  a  particle  traveling  at  10 
meters  per  second  (m/s)  uould  travel  1  cm  or  approximately  11  X  of 
the  sample  trap  length.   A  velocity  of  10  m/s  is  at  the  upper 
limit  of  the  particle  velocity^ distribution  and  uould  account  for 
only. a  small  percentage  of  the  particles  ejected  by  the  bed.   The 
majority  of  the  particles,  for  the  fluidization  conditions  used  in 
the  bed,  have  an  average  ejection  velocity  of  1-2  m/s,  based  on 
data  from  George  and  Grace  [81. 

TORQUE  RNALYSIS 

To  determine  the  minimum  required  torque  output  of  the 
solenoid,  tuo  analysis  uere  performed.   The  first  analysis 
calculates  an  average  torque  required  to  produce  the  desired 
velocity.   The  second  analysis  uses  the  torque  data  from  the 
selected  solenoid  and  calculates  the  expected  closure  time  and 
suing  time  for  the  paddle  arms. 
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First  Analysis 

To  close  the  sample  trap  in  1  ms,  the  angular  velocity  of  the 
paddles  at  closure  can  be  calculated  as- 


u  =   U2 

R  tZ     --.-■  (D.l) 


uhere-  u      =   Angular    velocity   of    paddles 

-:.  J*Z   =   Uidth   of    sample   trap 

R     =   Paddle   pivot  to   sample   trap   distance 
+2  =  ~61osure   time   of    sample   trap 


Using   the    values    of    UZ=    0.00137   m,    R=    0.0413   m,    and   t2=    1    ms, 
equation   (D.l)    results    in   an   angular   velocity   of   47.70   rad/sec. 
The    radial    acceleration    required    to    attain    this    velocity    through    a 
deflection   of ^74   radians    (45   Deg),    is   calculated   by: 


a  =      u2  (D-2) 


2    8 


uhere:     a  =  Angular  acceleration  of  paddles 
u  =  Angular  velocity  of  paddles 
8  =  Angular  deflection  of  paddles 


Using  the  result  of  equation  Dl  and  8="74  radians,  equation  D2 
gives  the  required  angular  acceleration  of  the  paddles  as  1448.5 
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rad/s/s.   Containing  this  result  with  the  result  for  the  total  mass 
moment  of  inertia  from  Appendix  A,  the  required  average  torque 
output  of  the  solenoid  can  be  calculated  as- 


T  =  I  a  (D3) 


uhere;     a  =  Angular  acceleration  of  paddles 
I  =  Mass  moment  of  inertia 
T  =  Required  solenoid  torque 


The  resultant  value  for  the  torque  (T)  is  0.87  N  m  (7.7  Lbf  in). 
A  geometric  average  of  the  selecred  solenoid's  torque  output,  as 
shoun  in  Fig.   D-la,"  indicates  an  average  of  1.0  N  m  (8.8  Lbf  in) 
Therefore,  the  selected  solenoid  has  sufficient  torque  output  to 
achieve  the  desired  paddle  velocity. 


SECOND  ANALYSIS 

For  the  second  analysis,  it  uill  be  assumed  that  the  torque 
output  of  the  selected  solenoid  can  be  Modeled  as  a  linear  spring. 
Fig.   D-lb  shous  the  torque  output  of  the  solenoid  as  a  function 
of  angular  displacement .   Fron  Fig.   D-lb,  a  spring  constant  of  k: 
0.72  N  m/rad  (0.11  Lbf  in/deg)  can  bs  used  to  approximate  the 
tsr^je  curve.   The  spring  constant  is  derived  from  the  slops  of 
the  tor que- angular  displacement  curve. 


LS4 


12 


10- 


o 

S3 


o 


4  - 


Goemetric  Average 


End  of  Stroke 


Triangle 
Approximation 
of  Torque 
Curve 


10 


Z0 


30 


40 


Angular  Displacement  (Degrees) 
Fig.  D-la   Torque  Output  of  Rotary  Solenoid 
Shouing  Triangle  Approximation. 

1Z  , 


10- 


Slope    of    Approximation 


•a 

S3 


a 

c 
o 


0  10  -  Z3  30  40 

Angular  Displacement  (Degrees) 

Fig.  D-lb   Torque  Output  of  Rotary  Solenoid 
Showing  Slope  of  Approximation 
Curve  Used  to  Determine  Spring 
Constant . 
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The  system  can  then  be  modeled  as  a  simple  rotational  mass  - 
spring  system.   The  differential  equation  of  notion,  general 
solution,  and  boundary  conditions  are: 


Diff  Eqn:         I  8  +  k  8  =  0  (D.4a) 


Gen  Soln--         8  =  CI  sin(m  t)  +  C2  cos (  pi  t)         (D.4b> 
uhere:    m   =   k 


b.c.:         1)    t  =  0:        8=0 

Z)    t  =  0:       8  =  1.S5  rad 


Solving   equation   (0.4b)   using   the   given   boundary   equations,    gives 
the   following   expressions   for   any   angular   displacement    (8)    and 
angular    velocity    (8): 


8   =    1.S5   cos(m   t)  (D.5a) 

8    =    -1.65   m    sin(m    t)  (D.5b) 


Knouing  that  the  solenoid  stroke  is  '^74  radians,  the  total  solenoid 
actuation  time  can  be  obtained  from  equation  (D.Sa)  by  setting  0  = 
( 1 . S5  -  '^74 ).   The  total  actuation  time  uas  determined  to  be  Z9.2 
ns.   Inserting  this  value  for  tine  into  equation  (0.5b),  a 
rotational  velocity  of  48. S  rad/s  uas  found.   This  velocity  is 
very  close  to  the  desirea  velocity  obtained  from  equation  (D.l). 
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DYNAMIC  ANALYSIS 

The  desire  to  determine  uhether  the  debris  of  a  specific 
bubble  has  been  trapped  by  the  sampling  device  requires  the 
determination  of  the  total  actuation  time  for  the  apparatus. 
-    Knouing.this  value,  the  delay  between  the  detection  of  a  bubble 
and  the  closure  of  the  sample  trap  can  be  determined,  and  hence, 

. uhether  the  particles  from  the  detected  bubble  uere  uithin  the 

vicinity  of  the  trap. 


=-.e--  i-_  Appendix  £L  determined  that  the-  actual  closure  time  for  the 
sample  trap  uas  1.44  ms.   Using  this  value  in  equation  (0.1),  a  = 
rotational  velocity  at  the  end  of  the  stroke  of  only  33.12  rad/sec 
-is  calculated.   Because  the  mass  moment  of  inertia  for  the 
solenoid  uas  unknown  in  the  beginning,  the  analysis  ignored  it. 
Uith  the  determintion  of  the  actual  closure  time  in  Appendix  C, 
the  total  moment  of  inertia,  including  the  solenoid  can  be 
approximated.   By  iterating  equations  (0.5a)  and  (0.5b)  uith  8  = 
33. 1Z  radians/sec,  a  neu  inertia  value  of  1.27  (-03)  kg  m  is 
determined.   Using  this  inertia  value  in  equation  (0.5b),  the 
total  actuation  time  is  42 . S  ms.   This  increase  in  actuation  time 
of  13.5  ms  is  relativily  close  to  the  manufactures  quoted 
actuation  time  for  the  solenoid  of  12  ms  and  suggests  that  the 
values  are  uithin  reason. 

Table  0.1  provides  a  summary  of  the  results  obtained  from  the 

dynamic  analysis. 
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TABLE  D.l 


Mass  moment  of  inertia- 


Paddles  (I)  5.991  (-04)  kg  n* 

2.037   (00)    lb   in1 


Paddles   +   Solenoid   (It)  1.27   (-03)    kg  m 

4. 32  (00)  lb  in2" 


Closure  time: 

Sample  trap  (t2)  1.44  ms 

Total  suing  <t3)  42. S  ms 

Summary  of  dynamic  analysis  results. 
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APPENDIX   E 

Bed  Particle  Size  Distribution  Analysis 

A  sample  of  particles  uas  obtained  from  the  fluidized  bed 
after  the  bed  had  been  operating  for  several  hours.   Sample  sizes, 
ranging  from  1.2  to  2.8  kg,  uere  then  taken  from  the  central  area 
of  the  bed.   Each  sample  uas  then  sieved  through  a  series  of  14  US 
standard  uire  mesh  sieves  using  a  Tyler  Industrial  Products  Model 
RX-24  portable  sieve  shaker  for  30  minuites.   The  contents  of  each 
sieve  uas  then  ueighed,  using  a  Torsion  Balance  Co  TORBAL  scale, 
to  an  accuracy  of  0.01  gram.   The  resulting  average  particle  size 
distribution  is  listed  in  Table  E.l  and  shoun  in  fig.   E-l. 
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TABLE  E.i 


US 

Seive 

Ueight 

Percent 

Cumulative 

Standard 

Size 

In  Seive 

of  Total 

Percent 

Seive  No. 

upi 

grans 

Sample 

18 

1000 

8.49 

0.31 

0.31 

za 

350 

Z4.S0 

0.90 

1.21 

30 

S00 

101. 9S 

3.73 

4.35 

35 

500 

110.85 

4.0G 

9.01 

45 

355 

17S.3S 

S.46 

15.46 

Ed 

Z97 

430.77 

15.78 

31.  Z4 

S3 

Z50 

385.59 

14. 1Z 

45.35 

.78 

Z12 

~oc  32 

14.43 

53.34 

80 

130 

4S3.Z0 

17.13 

77.03 

100 

149 

339.93 

1Z.4S 

89.43 

1Z0 

1ZS 

125.01 

4.58   - 

34.06 

140 

106 

109.45 

4.01 

98.36 

170 

90 

4G.S7 

1.71 

99.77 

Z00 

75 

5.91 

0.ZZ 

39.99 

<7S 

0.30 

0.01 

1 00 . 00 

Particle  Specific  Gravity  =8.1 

Average  particle  size  distribution  of  bed  material  in  grams 
and  percentage  of  total  ueight  using  sieve  analysis. 
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APPENDIX   F" 

Mean  Bed  Flou  Velocity  Determination 

To  determine  the  mean  velocity  of  air  flowing  through  the 
bed,  the  ASME  report  on  fluid  meters  [17]  was  used.   Based  on  this 
papery  equation  (F.l)  uas  used  for  determining  the  mass  flou  rate 
of  air  through  the  bed. 


2.  h. 

U  =   0.0S370Z    (KYFd      )    (  6    dP    )  (F.l) 


uhere: 


U  =  Mass  flou  rate  (lbm/sec) 
d  =  Orfice  diameter  (inches) 
K  =  Flou  coefficient 

F  =  Thermal  expansion  factor  (1  for  air) 
Y  =  Expansion  factor 
dP  =  Differential  pressure  (inches  uater) 
£  =  Density  of  air  (ahead  of  orfice) 


The  flou  coefficient  (K)  is  a  function  of  Reynold's  number,  the 
orfice  diameter  (d),  and  the  pipe  diameter  (D).   For  the  pressure 
tap  configuration  used  in  the  MIT  atmospheric  fluidized  bed  (1-0  / 
1/2-D)  the  flou  coefficient  is  calculated  using  eauation  (F.2). 


uhere: 


K  =  Ko  +  1000  b 

VB  Re'  (F.2) 


Re  =  Reynolds  No. 

B  =  d/0  =  (orfice  dia  /  pipe  dia) 

Ko  =  (0.S014  -  0.01352  D  H  ) 

+  (0.37S0  +  0.07257  D  '*) 

(      0.00025      +  B*   +  1.5  b'4  ) 
DZ  BZ+0. 00025  D 


b  =  (0.0002  +  0.0011  ) 
D 

+  (0.0038  +  0.0004) 


(  8*  +  (IS. 5  +  5  0)  b'   ) 


The  expansion  factor  (Y)  is  determined  from  equation  (F.3)  and  is 
a  function  of  the  diameter  ratio  (B),  the  ratio  of  specific  heats 
(s),  and  the  ratio  of  the  differential  pressure  to  inlet  pressure. 


Y  =  1  -  (0.410  +  0.350  BH    )   dP 

PI  s  (F.3) 

uhere: 

PI  =  Inlet  pressure 

s   =  Ratio  of  specific  heats  (1.4  for  air) 
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Mass   flou   rate    is   a   function   of    velocity,    and   therefore,    by 
definition   a   function   of    Reynold's    number.      Equation    (F.4)    is  used 
to   determine   mass    flou   rate   as   a   function   of    viscosity   and 
Reynold's    number. 


U   =   u    02 

C   Re  (F.4) 


uhere 


u  =  Absolute  viscosity  ( lbm/f t  sec) 
C  =  Constant  ( 15.28) 


This  requires  that  the  solutions  of  equation  (F.l)  and  equation 
(F.4)  be  iterated  until  the  Reynold's  numbers  converge.   The 
determination  of  mean  air  velocity  is  obtained  from  the  mass  flou 
rate  (U)  using  equation  (F.5). 


V  =   U 


uhere: 


e  a  (F.s) 


C    -    Density  of  air  ( lbm/f t  > 
A  =  Area  of  bed  (ftZ) 
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To  perform  these  calculations,  a  computer  program  uas  used.   The 
program  (APPENDIX  G)  is  uritten  in  HP  BASIC  2.0  for  running  on  an 
HP  981S  series  200  micro-computer.   Convergence  usually  required 
tour  to  seven  iterations  and  should  be  correct  to  uithin  0.5Z. 
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APPENDIX      G 

This    Appendix    lists    the   BASIC    computer   program  used   to 
calculate   mean   bed   velocity   as   discussed    in   Appendix    F.       It    is 
uritten    in   HP    BASIC   Z.0. 
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10 

2a 

30 
40 

50 
SO 
70 

30 
30 

100 

110 

,?atn 

123 

130 

HO 

150 

150 

170 

180 

190 

200 

210 

220 

230 

240 

250 

ISO 

270 

230 

230 

300 

310 

320 

330 

3*0 

350 

3S0 

370 

330 

3=0 

400 

410 

420 

430 

4.40 
450 
450 

470 
430 

4S0 

500 
510 
520 
530 
542 
550 
550 
370 
530 
530 


MAIN  Program 

Metering  of  gases  oy  means  of  the  ASME  sauare-eaced  orfice  with 
1-0  1/2-0  taps.   Reference  Fluid  Meiers  Their  Theory  ana 
Application,   ASME  Report  5th  edition,  ASME  New  York,  NY.  1371 
Program  must  oe  altered  for  Orfics  diameter  other  than  7.071  in 
anc  pipe  diameter  otner  than  10.02  in.  (see  line  130) 


INTE-ER  Answer 

REAL  Dl.D2.U\X,Yl,?:,?la,7l.S,F_aei.5el3,Area_bed.vmu,ReynoIc-s_nG,Rhc,Rhol 

Velocity 

DIM  HomeSC2:,ClearSC21 

DATA  7.371,10.02.1.4,11.511.50000. ,75. 

READ  02,Dl,S,Area_ped,Reynolas_no,Patm 

CiearS«CHRs<255 >&CHRS(7S) 

Hone=«CHRS( 255 !iCHRS(  84 ) 


!  CLEAR  screen 
!  HOME  screen 


Input  variables 


INPUT  'Enter  static  pressure  PI  (cm  Hg > :  ",P1 

INPUT  'Enter  pressure  orop  d?  (  incnes  water):  "  ,.=_del 

INPUT  "Enter  air  temperature  (Oegrees  ?)■■    *,Tl 


Compute  data 

OUTPUT  2: HomeS ; 

OUTPUT  2; Clears; 

3eta-02/Dl 

Pia-?l-!-?atm 

Pl-(Pla)/2.54 

T2-T1+4SS.S7 

Rho«FN0ensity<  Fatm,T2 ) 

Rhol-FN0en«ity(Pla.T2  > 

Umu"FNUl 5COSlty(T2 ) 

Yl-FNExoan(Beta,P1.3,P_del  ) 

K«FNFIow_cpef f (Reynolas_no.3eta.Dl  ) 

U-FNMass_f low(D2,K,Yl  ,Rhoi,?_del  ) 

Ro-FNReyn(U,vmu,D2 ) 

1-    r:5S(  Rd-Reynoi  cs_no  )>  1 .  THEN 

Reynolds_no»Rd 

SOTO  350 
END  IF 

velocity-wY(Rho*Area_bed) 
PRINT  "Velocity-  ^Velocity!"   ft/s* 
INPUT  "Enter  ( 1  )  to  continue,  (0)  to 
IF  Ansuier-0  THEN  STOP 
IF  AnsuerOl  THEN 

PRINT  "Value  must  be  either  1  or  0;  try  again 

50TC  450 
END  IF 
SOTO  150 
END 


Home  display 
Clear  display 
Calculate  beta 
Absolute  pressure 
Change  to  inches 
Convert  to  R 
Air  density  bed 
Air  density  upstream 
Air  viscosity 
Expansion  factor  Y 
Flow  coefficient  K 
Mass  flow  rate 
Reynolds  2  calculated 
Check  accuracy 


!  Ueiocity  calculation 
!  Output  results 
stop" , Answer 


MAIN. 


.aegin  FUNCTIONS 


i 

2E-  "NExcani seta  .  =  1 . £ . =_cei ) 
Holri-r*_dei/< Pi ♦5*13. 52=5  ) 

Y»l-< .4i-.33*5eta*4)*Hoid 
RETURN  V 


rector 
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503  FNEND 

513  ! 

520  ! 

£50  DEF  FNFiou_coef f ( Reynolds_no,Beta,01 )          !  Flow  coefficient  K 

543  Ko«.6B14-.ai352*01*(-.2S)+< .575-r  . 37237*01 *( -.2S ) )•< .30025/C 01*01 *Se:a*Se 

ta+.3002S*Ol >+Beta*»+l.S*eeta"lB ) 

5=3  3«.3002+.00ii/Dl+< .0038+. 3004/01 >»<Beta*8eta+< IS-S+S*01 >*8eta~lG > 

553  K*KoTl000*3/SQP.(3eta*neynoIds_no  > 

570  RETURN  K 

530  FNENO 

5=0  ! 

733  ! 

713  OEF  FHflass_flow<02,K,Yl,Rhol,P.jJaI )           !  Mass  flow  rate  U 

723  y«.3SS732*02*C2*K*Yl*5GR(Rhoi*P_dei ) 

730  RETURN  U 

740  FNENO 

7E3  ! 

753  ! 

770  OEF  FNReyn(U,Vnu,02>                          !  Reynolds  $  calculated 

730  Rd»iE.2S*U/(Umu*02) 

733  RETURN  Sa 

303  FNENO 

313  ! 

323  ! 

333  DEF  FN0ensity(Pet.-n.T2 )                        !  Calculate  air  density 

343  Rho*.522406*Patn/T2 

3E3  RETURN  Rho 

333  FNENO 

3T0  ! 

383  ! 

390  DEF  FNUiscosity< T2 )                            I  Calculate  air  viscosity 

333  U«u«7.3B3E-7*T2*1.5/<T2+lS8.G > 

310  RETURN  Vnu 

323  FNENO 
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APPENDIX   H 

This  Appendix  contains  the  parts  list  for  the  sampling  system 
in  three  tables.   Table  H.l  is  a  list  of  components  for  the  sampling 
apparatus.   Table  H.2  is  a  list  of  components  for  the  electrical 
circuit  uhile  Table  H.3  is  for  the  vacuum  system. 
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Table   H.l 


Item   No.  Component  Description 

1  Solenoid  LEDEX   Size   SS   45  Degree   Right 

Hand   Stroke   Rotary   Solenoid 
Part   No.    S-8204-0Z9 
LEDEX   Inc. 

301    Scholz   Or,    P.O.    Box    427 
Vandal i a,    Oh     45377 
(513)    838-3GZ1 

Z  Paddle  Foam,    Bassuood   epoxy   laminate 

with   Maple   mounting   blocks 
0.5  X   1.2   X   11. S   inches 

3  Interface  Aluminum  Cylinders   with   set    screu 

fastener.      R=   0.S75    in   H=    1.5    in 

4  Solenoid  Aluminum   mounting   plate 
Mounting  1/4   X   Z.8   X   Z.8   inches 

5  Extension  Aluminum   bar 

Bar  1/4  X  3/4  X  1Z.G  inches 

5         Trap  Aluminum  Plate 

Mount  1/4  X  1.1  X  1.5  inches 

7  Sample  1/1S  inch  Aluminum  plate 
Trap             Inside  Dimensions^ 

0.5  X  1.5  X  3.5  inches 

8  Vertical         Aluminum  bar 
Mounting         1/4  X  3/4  X  24  inches 

Slide  Adjustment  holes  drilled  every 

0.5  inch 

9  Base  Tripod  Aluminum  Structure 
Structure 


List  of  components  for  sampling  accaratus. 
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TABLE  H.Z 

Item  No.       Component  Description 

1  Bl      Bridge  rectifier  assembly,  Silicon, 

LEDEX  Part  No.  121011-001,  includes 
arc  suppressor  in  unit. 

Z  Dl      Arc  Suppressor,  not  needed  if  above 

rectifier  assembly  used. 
LEDEX  Part  No.  1ZZSSS-001 

3  Ml  Solenoid  (see  Table  H.l) 

4  Rl  Resistor,  Z  Mohm,  IX,  Z  Uatt 

5  RZ  Resistor,  100  kohm,  IX,  1  Uatt 
S  R3  Resistor,  ZS0  ohm,  10X,  50  Uatt 

7  SI      Suitch,  SPST,  Push  Button, 

10  A,  250  V 

8  SZ      Suitch,  SPST,  Toggle,  10  A,  ZS0  v 
List  of  components  for  solenoid  pouer  supply. 
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TABLE  H.3 

Item  No.  Component       Description 

1  VI  3/8  inch  Ball  Valve 

Z  PI  Eductor 

3  Fl  Screen  Filter,  3Z0  un  mesh 

4  CI      Sample  Container,  Snail  Plastic 

Bottle,  1  Pt 

5  CZ      Sample  Trap  (see  Table  H.l) 
S  LI       1/4-  inch  Polyflou  Tubing 

List  of  components  for  vacuum  system. 
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APPENOI X 


This   Appendix    contains    a   complete   listing   of    all    data 
obtained   during   this    study.       The   data   is   arranged 
acording   to   the   distance   of    the   trap   height    above   the 
distributor. 
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DATA    for   Urn    DETERMINATION 


Hb 

PI 

Pb 

dP 

T 

( inches ) 

(cm  Hg) 

(en  H20) 

(n  H20) 

(F) 

10.75 

7.2 

51.9 

1.805 

52.0 

10.75 

5.8 

50.9 

1.517 

52.0 

10.25 

6.2 

49.1 

0.932 

62.0 

9.75 

5.0 

48.8 

0.772 

52.0 

9.50 

5.5 

47.5 

0.515 

62.0 

9.25 

5.1 

44.9 

0.335 

62.5 

8.75 

4.8 

42.9 

0.255 

63.0 

8.50 

4.5 

"39.9 

0.177 

63.0 

8.50 

4.0 

39.5 

0.125 

63.5 

8.50 

3.7 

33.9 

0.107 

64.0 

8.50  S_ 

3.5-  j_ 

32.0 

0.094 

64.0 

8.50 

3.  3 

30.1 

0.084 

64.0 

8.50 

2.4 

22.3 

0.043 

64.5 

8.50 

2.0 

18.5 

0.027 

65.0 

8.50 

1.5 

14.5 

0.017 

65.0 
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Trap   Height  above  Distributor:  en  (in) 

Bed   Height  above  Distributor"  en  (in) 

Trap  Height  above  bed  surface:  en  (in) 

Pi:  en  Hg   (in  Hg) 

dP:  en  uater   (in  uater) 

Temp:   C  (F) 


■z. 

30. 9G 

(12.19) 

= 

25. 4 

(10) 

- 

5.56 

(2.19) 

- 

5.8 

(2.23) 

■=. 

1.234 

(0.486) 

- 

19.  S 

(S7.0) 

SAMP   It 

Ueight    of    Sanple 

( grans ) 

1 

3. IS 

2 

2.36 

3 

1.88 

4 

2.43 

5 

2.48 

6 

2.10 

7 

1.84 

8 

2.S7 

9 

2. 57 

10 

2. IS 

11 

2. IS 

12 

1.03 

13 

2.46 

14 

1.13 

IS 

1.21 

flVG    =  2.09   grans 
STN  DV  =  0.S0  grans 


Trap   Height  above  Distributor:  en  (in) 

Bed   Height  above  Distributor:  en  (in) 

Trap  Height  above  bed  surface-"  en  (in) 

Pi:  en  Hg   ( in  Hg) 

dP:  en  uater   (in  uater) 

Tenp:   C  (F) 


30.96 

(12.19) 

26.67 

(10.5) 

4.29 

(1.69) 

6.8 

(2.68) 

1.887 

(0.743) 

18 

(84.5) 

SAMP   * 

Ueight    of    Sanple 

( grans ) 

16 

7.72 

17 

9.12 

18 

3.83 

19 

4.04 

20 

5.66 

21 

3.08 

22 

8.86 

23 

4.50 

24 

6.55        fWG 

ZS 

5.84        STN 

=  5.96   grans 
STN  DV  =  2.12   grans 


i  no. 


0.52 

0.55 

0.54 

0.54 

0.76 

0.53 

0.61 

0.47 

0.85 

AVG    =  0.5Z  grans 

0.61 

STN  OV  =  0.11  grans 

Trap   Height  above  Distributor:  cm  (in)  =  37.31  (14.69) 

Bed   Height  above  Distributor:  cn  (in)  =  25.4  (10) 

Trap  Height  above  bed  surface:  cn  (in)  =  11.91  (4.59) 

Pi:  cn  Hg   (in  Hg)  =  5.8  (2.28) 

dP-    cn  water   (in  uater)  =  1.201  (0.473) 

Tenp:   C  (F)  =  20.5  (69.0) 


SAMP  ft     Weight  of  Sanple 
(grans ) 


25 
27 
28 
29 
30 
31 
32 
33 
34 
35 


Trap   Height  above  Distributor:  cn  (in)   =  37.31   (14.69) 

Bed   Height  above  Distributor: 

Trap  Height  above  bed  surface: 

P 1 :  cn  Hg   (  i  n  Hg ) 

aP:  cn  water   (in  water) 

Tenp:   C  (F) 


SAMP  ft  weight  of  Sanple 
( grans ) 

36  2 . 87 

37  1 . 42 

38  1 . 89 

39  1 . 34 

40  1.18 

41  1.52 

42  1 . 57 

43  1.53 

44  1.56   AVG    =  1.6S   grans 

45  1.57    STN  OV  =  0.46  grans 
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cn  (in) 

s 

25.57 

(10.5) 

cn  (in) 

= 

10.54 

(4.19) 

= 

5.3 

(2.28) 

= 

1.897 

(0.747) 

= 

19 

(66.0) 

Trap   Height  above  Distributor-'  en  (in)  =  37.31  (14.S9) 

Bed   Height  above  Distributor:  en  (in)  =  27.94  (11.0) 

Trap  Height  above  bed  surface:  en  (in)  =  9.37  (3. S3) 

Pi:  en  Hg   (in  Hg)  =  S.8  (2.S8) 

dP:  en  uater   (in  uater)  =  1.887  (0.743) 

Tenp:   C  (F)  =  13  (So.0) 


SAMP  Jt  Weight  of  Sanple 
( grans ) 

4G  4.3Z 

47  2 . 48 

48  3.12 
43  2.83 

50  4.4S 

51  2.37 

52  4.46 

53  4.95 

54  4.15   AVG    =  3. 51   grans 

55  2.9G   STN  DV  =  0.95  grans 


L8S 


TRAP  HEIGHT  =  20.25  in 


Height    of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

for 

for 

for 

Tap   above 

Samples 

Samples 

Samples 

Samples 

Distributor 

55-75 

35-95 

55-55 

76-85 

( i  nches > 

(en   H20) 

(ca   H20) 

(cm   H20) 

(cm   H20) 

1.6 

52.9 

54.0 

55. 1 

58.5 

2.5 

48.8 

49.3 

51.7 

52.8 

=  -  -~ 

3.5 

40.5  .__ 

41.3 

43.8 

44.7 

4.5 

34.9 

35.8 

38.9 

40.0 

5.5 

27.8 

29.7 

32.5 

33.6 

' 

S.S 

21.2 

22.3 

25.5 

27.5 

— 

7.5 

14.5 

15.1 

19.1 

20.7 

8.5 

8.3 

9.7... 

12.5 

14.8 

1    ■■■ 

3«y 

2.2 

3.5     E 

5.4 

8.5 

10.5 

0.0 

0.1    : 

1.5 

2.4 

.  .11.5 

0.0 

0.0 

0.0 

0.3 

12.5 

0.0 

0.0 

0.0 

0.0 

13.5 

0.0 

0.0    , 

0.0 

0.0 

14.5 

0.0 

0.0 

0.0 

0.0 

TRAP  HEIGHT  =  23.25  in 


Height    of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

for 

for 

for 

Tap   above 

Samples 

Samples 

Samples 

Samples 

Distributor 

206-215 

196-205 

186-195 

176-185 

( inches ) 

(cm   H20) 

(cm   H20) 

(cm   H20) 

(cm   H20) 

1.5 

52.0 

53.4 

56.3 

57.7 

2.5 

48.9 

49.5 

51.8 

52.2 

3.5 

40.7 

41.5 

43.6 

44.7 

4.5 

34.7 

35.8 

38.8 

40.0 

5.5 

28.3 

29.3 

32.0 

33.3 

6.5 

21.5 

22.9 

25.6 

26.6 

7.5 

14.9 

16.6 

18.8 

20.3 

8.5 

8.4 

9.8 

12.9 

13.8 

9.5 

2.2 

3.3 

6.2 

7.5 

10.5 

0.0 

0.  1 

1.3 

2.9 

11.6 

0.0 

0.0 

0.0 

0.3 

12.5 

0.0 

0.0 

0.0 

0.0 

13.5 

0.0 

0.0 

0.0 

0.0 

14.5 

0.0 

0.0 

0.0 

0.0 
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51.44 

(20.25) 

29.11 

(ll.S) 

22  .23 

(8.75) 

22.55 

(8.88) 

8.3 

(3.21) 

5.985 

( 2 . 75 ) 

IS. 5 

(  52 . 0 ) 

Ueight  of  Sample 

(grams) 

1.47 

1.35 

1.24 

1.29 

0.74 

1.38 

1.38 

1.47 

1.58 

AVG 

=  1.32 

grams 

1.33 

STN 

DV  =  0.23 

grams 

Trap   Height  above  Distributor:  cm  (in) 

Bed   Height  above  Distributor:  cm  (in) 

Trap  Height  above  bed  surface'-  cm  (in) 

Static  Bed  Height:  cm  (in) 

Pi:  cm  Hg   (in  Hg) 

dP'    en   uater   (in  uater) 

Temp:   C  (F) 


SAMP  * 


75 
77 
78 
79 
30 
81 
82 
83 
84 
85 


Trap   Height  above  Distributer:  cm  (in) 

Bed   Height  above  Distributor:  cm  (in) 

Trap  Height  above  bed  surface:  cn   (in) 

Static  Bed  Height:  cm  (in) 

Pi:  cm  Hg   (in  Hg) 

dP :  cm  uater   (in  uater) 

Temp:   C  (F) 


SAMP  #  Ueight  of  Sample 
( grams ) 

56  0.95 

57  0-51 

58  1.83 

59  0.84 

50  1 . 02 

51  0.75 

52  0 . 57 

53  1.05 

64  0.60   AVG    =0.32   grams 

55  0.88    STN  DV  =  0.36  grams 
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51.44 

(20.25) 

27.31 

(10.75) 

24.13 

(9.5) 

22.56 

(8.88) 

7.3 

(  2 .  87 ) 

4.910 

(1.933) 

16.5 

(62.0) 

Trap   Height  above  Distributor:  cm  (in) 

Bed   Height  above  Distributor:  en  (in) 

Trap  Height  above  bed  surface:  cm  (in) 

Static  Sed  Height:  cm  (in) 

Pi:  cm  Hg   (in  Hg) 

dp:  cm  uater   (in  water) 

Temp:   C  (F) 


SAMP  It 


86 
87 
38 
89 
90 
91 
9Z 
93 
94 
95 


51.44 

(Z0.Z5) 

ZS.S7 

(10.5) 

Z4.7S 

(9.75) 

ZZ.56 

(8.88) 

5.5 

(Z.56) 

3.ZS1 

(1.Z84) 

16 

(51.0) 

._ . 

Weight  of  Sample 

(grams ) 

0.07 

0.08 

0.17 

0.04 

0.08 

0.04 

0.07 

0.13 

0.07 

AVG     =  0.09 

grams 

0.17 

STN  DV  =  0.05 

grams 

Trap   Height  above  Distributor:  cm  (in) 

Bed   Height  above  Distributor:  cm  (in) 

Trap  Height  above  bed  surface-"  cm  (in) 

Static  Bed  Height:  cm  (in) 

P 1 :  cm  Hg   (  i  n  Hg ) 

dP:  cm  uater   (in  uater) 

Temp:   C  (F) 


SAMP  8 


66 
67 
68 
53 
70 
71 
7Z 
73 
74 
75 


51.44 

(Z0.Z5) 

Z5.4 

(10.0) 

Z6.03 

(10. Z5) 

ZZ.56 

(8.88) 

6.Z 

(  Z . 44 ) 

Z.581 

(1.016) 

17 

(63.0) 

ight  of  Sample 

( grams ) 

0.05 

0.15 

0.14 

3.07 

0.11 

0.10 

0.18 

0.08 

0.11 

AVG 

=  0.1Z 

grams 

0.16 

STN 

DV  =  0.04 

grans 

!  Hq 


TRAP  HEIGHT  =  IS. 12  in 


Height    of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

for 

for 

for 

Tap   aoove 

Samples 

Samples 

Samples 

Samples 

Distributor 

106-115 

126-135 

95-105 

116-125 

( inches) 

(cm  H20) 

(cm  H20) 

(cm   H20) 

(cm   H20) 

1.6 

52.5 

54.4 

55.5 

58.5 

2.S 

48.7 

49 . 9 

51.3 

52.9 

3.S 

40.3 

41.8 

.43.3 

44.4 

4.5 

34.5 

35.9 

37.9 

39.  S 

S.5 

27.9 

29.2 

31.3 

33.9 

S.S 

20. S 

22.8 

25.2    -- 

27.1 

7.S 

13.8 

15.5 

18.9 

20.5 

8.5 

7.8 

9.4 

12.3 

14.5 

9.S 

1.8 

3.2 

5.  3 

7.S 

13.5 

0.0 

0.3 

0.9 

2.5 

11.5 

0.0 

0.0 

0.0 

0.3 

12.5 

0.0 

0.0 

0.0 

0.0 

13.5 

0.0 

0.0 

0.0 

0.0 

14.5 

0.0 

0.0 

0.0 

0.0 

TRAP  HEIGHT  =  18.25  in 


Height    of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

for 

for 

for 

Tap   above 

Samples 

Samples 

Samples 

Samples 

Distributor 

146-155 

136-145 

155-165 

166-175 

( inches ) 

(cm  H20) 

(cm   HZO) 

(cm   H20) 

(cm   H20) 

1.5 

52.7 

54.3 

55.8 

57.9 

2.6 

49.0 

49.9 

51.5 

52.6 

3.5 

41.3 

41.6 

43.3 

44.9 

4.6 

34.5 

36.1 

38.1 

39.9 

5.6 

28.0 

29.3 

31.5 

33.5 

6.6 

21.5 

22.5 

Z4.9 

26.6 

7.6 

14.9 

15.8 

18.8 

20.7 

8.5 

8.5 

9.9 

12.2 

14.5 

9.5 

Z.4 

3.4 

5.0 

7.8 

13.5 

0.0 

0.1 

3.9 

Z.5 

11.5 

0.3 

0.3 

3.3 

0.3 

12.5 

0.0 

0.0 

0.0 

0.3 

13.5 

0.0 

0.0 

0.0 

3.3 

14.5 

0.3 

0.3 

0.0 

3.0 

133 


Trap   Height  above  Distributor:  cm  (in)   =  40.95  (1G.1Z) 

Bed   Height  above  Distributor:  en  (in)    =  28.58  (11.Z5) 

Trap  Height  above  bed  surface:  en  (in)    =  12.37  (4.88) 

Static  Bed  Height:  cm  (in)  =  22. SB  (8.88) 

Pi:  cm  Hg   (in  Hg)  =  8.0  (3.15) 

dP:  cm  uater   (in  uater)  =  S.8S8  (2.7) 

Temc:   C  (F)  =  IS  (B1.0) 


SAMP  ft     Ueight  of  Sample 
( grams) 


US  4.18 

117  3.03 

118  2.3Z 

119  1.47 

120  2.34 

121  2. 54 

122  1.73 

123  2.12 

124  1.96   flVG    =2.41   grams 

125  2.38   STN  DV  =  0.76  grams 


= 

40.95 

(16.12) 

= 

27.30 

(10.75) 

= 

13.  65 

(5.38) 

= 

22.56 

(8.88) 

= 

7.4 

(2.91) 

= 

4.808 

(1.893) 

= 

16 

(61.0) 

Trap   Height  above  Distributor-"  cm   (in) 

Bed   Height  above  Distributor-'  cm  (in) 

Trap  Height  above  bed  surface:  cm  (in) 

Static  Bed  Height:  cm  (in) 

Pi:  cm  Hg   (in  Hg) 

dP :  cm  uater   (in  uater) 

Temp:   C  (F) 


SAMP  3  Ueight  of  Sample 
( grams ) 

96  1.33 

97  0.81 

98  1.37 

99  1.56 

100  1.42 

101  1-43 

102  1.27 

103  2.05 

104  1.21    AVG    =1.36  grams 

105  1-11    STN  DV  =  0.32   grams 


191 


Trap   Height  above  Distributor:  en  (in)   =  40.35  (IS. 12) 

Bed  Height  above  Distributor:  en  (in)    =  2S.E7  (10.5) 

Trap  Height  above  bed  surface:  en  (in)    =  14. Z8  (5.62) 

Static  Bed  Height:  en  (in)  =  22.55  (8.88) 

Pi:  en  Hg   (in  Hg)  =  S.S  (2.S0) 

dP:  en  uater   (in  water)  =  3.252  (1.283) 

Tenp:   C  (F)  =16  (51.0) 


SAMP  * 


125 
127 
128 
129 
130 
131 
132 
133 
134 
135 


weight  of 

Sanple 

(grans 

) 

0.36 

0.31 

0.21 

0.27 

0.18 

0.37 

0.34 

0.29 

0.34 

AVG     =  0.28 

grans 

0.16 

STN  DV  =  0.08 

grans 

Trap   Height  above  Distributor:  en  (in) 

Bed   Height  above  Distributor:  en  (in) 

Trap  Height  above  bed  surface:  en  (in) 

Static  Bed  Height:  en  (in) 

Pi:  en  Hg   (in  Hg) 

dP :  en  uater   (in  uater) 

Tenp:   C  (F) 


SAMP  » 


106 
107 
138 
109 
110 
Hi 
112 
113 

114 
115 


40  95 

(16.12) 

26.04 

(10.25) 

14.91 

(5.88) 

Z2.56 

(8.88) 

6.3 

(  2 . 48 ) 

2.591 

(1.020) 

16 

(81.0) 

Ueight  of 

Sanple 

( grans 

) 

0.38 

0.12 

0.15 

0.17 

0.13 

0.18 

0.56 

0.05 

0.18 

AVG     =  0.22 

grans 

0.24 

STN  OV  =  0.15 

grans 
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Trap   Height  above  Distributor:  cm  (in) 

Bed   Height  above  Distributor-"  cn  (in) 

Trap  Height  above  bed  surface:  cn  (in) 

Static  3ed  Height:  cn  (in) 

Pi:  cn  Hg   (in  Hg) 

cP :  cn  uater   (in  water) 

Tenp:   C  (F) 


4S.3S 

(18.25) 

28.58 

(11.25) 

17.78 

(7.00) 

22.55 

(8.88) 

8.1 

(3.19) 

5.858 

(2.7) 

IS 

(51.0) 

SAMP  3 

Ue 

ight  of 

Sanple 

-  -- 

"--"• 

(grans 

) 

. -  - 

166 

- 

1.13 

---——-—. 

167 

1.18 

.  -  -.         '--■- 

158 

1.34 

165 

0.93 

170 

0.88 

171 

1.01 

172 

1.20 

173 

1.15 

174 

1.08 

AVG    =1.08 

grans 

175 

0.85 

STN  DV  =  0.16 

grans 

Trap   Height  above  Distributor:  cn  (in) 

Bed   Height  above  Distributor:  cn  (in) 

Trap  Height  above  bed  surface-'  cn  (in) 

Static  Bed  Height:  cn  (in) 

PI:  cn  Hg   (in  Hg) 

dP:  cn  uater   (in  uater) 

Tenp:   C  (F) 


SAMP  * 


156 
157 
158 

153 
1S0 
151 
162 
163 
164 
165 


46.36 

(18.25) 

27.30 

(10.75) 

19.06 

(7.5) 

22.56 

(8.88) 

7.3 

(  2  .  87 ) 

4.808 

(1.893) 

15 

(51.0) 

Ueight  of 

Sanple 

( grans 

) 

0.64 

0.48 

0.78 

0.58 

0.44 

0.58 

0.83 

0.57 

0.52 

AVG 

=  0.50 

grans 

0.62 

STN 

DV  =  0.12 

grans 
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Trap   Height  above  Distributor:  en  (in) 

Bed   Height  above  Distributor:  en  (in) 

Trap  Height  above  bed  surface:  en  (in) 

Static  Bed  Height:  cm  (in) 

PI:  en  Hg   (in  Hg) 

dP:  en  uater   (in  water) 

Tenp:   C  (F) 


SAMP  ft 


136 
137 
138 
13S 
140 
141 
142 
143 
144 
145 


4S.3G 

(18.25) 

2S.S7 

(10.5) 

19. S3 

(7.75) 

22.  SS 

(8.88) 

6.S 

(  2 . S0 ) 

3.244 

(  1 . 277  ) 

IS.  5 

(S2.0) 

Weight  of 

S, 

anple 

(grans 

) 

- 

-  0.03 

0.13 

0.0G 

0.05 

0.07 

0.06 

0.05 

0.05 

0.11 

AVG    =  0.08 

grans 

0.14 

STN  DV  =  0.04 

grans 

46.36 

(18.25) 

26.04 

(10.25) 

20.32 

(8.00) 

22.56 

(8.88) 

6.3 

(2.48) 

2.631 

(1.036) 

16.5 

(62.0) 

Trap   Height  above  Distributor:  en  (in) 

Bed   Height  above  Distributor:  en  (in) 

Trap  Height  above  bed  surface:  en  ( in) 

Static  8ed  Height:  en  (in) 

Pi:  en  Hg   (in  Hg) 

dP:  en  uater   (in  uater) 

Temp:   C  (F) 


SAMP  *  weight  of  Sanple 
( grans ) 

146  0.01 

147  0.01 

148  0.02 
1-19  0.02 

150  0.03 

151  0-32 

152  0.01 

153  0.01 

154  0.02    AVG    =  0.02   grans 

155  0-05    STN  DV  =  0.01   grans 
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Trap   Height  above  Distributor:  cm  (in) 

Bed   Height  above  Distributor:  cm  (in) 

Trap  Height  above  bed  surface:  cm  (in) 

Static  Bed  Height:  cm  (in) 

PI:  cm  Hg   (in  Hg) 

dP~    cm   uater      (in  uater) 

Temp:       C    (F) 


SfiW3    3 


59.06 

(23.25) 

Z7.34 

(11.00) 

31.12 

(12.25) 

Z2.56 

(8.88) 

8.2 

(3.23) 

S.8S8 

(2.7) 

IS. 5 

(SZ.0) 

Ueight  of 

S, 

ample 

(grams 
0.32 

) 

— —   : - 

—  '  - 

0.22 

-  - 

0.24 

""* 

0.52 

0.31 

0.42 

0.34 

0.31 

0.23 

nVG    =0.31 

grams 

0.21 

STN  DV  =  0.10 

grams 

176 
177 
178 
179 
130 
181 
182 
183 
134 
185 


Trap   Height  above  Distributor:  cm  (in) 

Bed   Height  above  Distributor:  cm  (in) 

Trap  Height  above  bed  surface:  cm  (in) 

Static  Bed  Height:  cm  (in) 

Pi:  cm  Hg   (in  Hg) 

dP :  cm  uater   (in  uater) 

Temp:   C  (F) 


SAMP  S  Ueight  of  Sample 
( grams ) 

.  186  0.25 

187  0.17 

188  0.14 

189  0.17 

190  0.12 

191  0.23 

192  0.16 

193  0.25 

194  0.18    AVG     =  0.18   grams 

195  0.17    STN  DV  =  0.04   grams 


19S 


59.06 

(23.25) 

27.18 

(10.70) 

31.76 

(12.50) 

22.56 

(8.88) 

7.3 

(  2 . 87 ) 

4.813 

(1.895) 

16 

(61.0) 

Trap   Height  above  Distributor:  en  (in)   =  53.06   (23.25) 

Bed   Height  above  Distributor:  en  (in) 

Trap  Height  above  bed  surface:  Cn  (in) 

Static  Bed  Height:  cn  (in) 

Pi:  cn  Hg   (in  Hg) 

dP:  cn  uater   (in  uater) 

Tenp:   C  (F) 


26.67 

( 10.50) 

32.39 

(12.75) 

22.56 

(8.88) 

6.5 

(2.56) 

3.236 

( 1 . 274 ) 

15 

(61.0) 

SAMP  * 

Uei 

.ght  of 

Sanple 

196 

- 

( grans 
0.05 

) 

■■<  '-— 

197 

0.07 

198 

0.07 

199 

0.04 

200 

0.06 

201 

0.09 

202 

0.04 

203 

0.07 

204 

0.06 

AVG    =  0.06 

grans 

205 

0.06 

STN  DV  =  0.02 

grans 

Trap   Height  above  Distributor: 
Bed   Height  above  Distributor: 
Trap  Height  above  bed  surface: 
Static  Bed  Height:  cn  (in) 
Pi:  cn  Hg   (in  Hg) 
dP :  cn   uater   (in  uater) 
Temp:   C  (F) 


SAMP  tt     Ueight  of  Sanple 
( grans ) 


cn  (in) 

=  59.06 

(23.25) 

cn  (in) 

=  26.04 

(10.25) 

cn  (in) 

=  33.02 

(13.00) 

=  22.56 

(8.88) 

=  6.3 

(2.48) 

=  2.519 

( 1.031) 

=  16 

(61.0) 

206  0.01 

207  0.02 

208  0.01 

209  0.01 

210  0.01 

211  0.01 

212  0.01 

213  0.02 

214  0.01     AVG     =  0.01   grans 

215  0.01     STN  DV  =  0.004  grans 
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TRAP  HEISHT  =  1Z.7S  in 


Height    of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

for 

for 

for 

Tap   above 

Samples 

Samples 

Samples 

Samples 

Distributor 

ZS6-2S5 

ZSS-Z75 

Z7S-Z8S 

Z8G-Z95 

(  i  riches  > 

(cn   HZO) 

(cn   HZO) 

(cn   HZO) 

(cm   HZO) 

l.S 

53.7 

54.8 

57.  Z 

58.8 

Z.S 

48.5 

49.9 

51.7 

53. Z 

3.S 

40.8      . 

4Z.1 

44.4 

45.3 

4.S 

34.3 

35.7 

38.7 

40.4 

s.s      ~   - 

Z8.1    __ 

Z9.7 

3Z.5   Z~- 

34.0 

5.S 

Zl.S 

Z3.8 

ZS.l 

Z7.8 

-    7.B 

14.7 

17. Z 

13.1 

Zl.S 

8.5 

8.1 

10.7 

13.5 

1S.Z 

9.5 

Z.l 

4.5 

5.8 

9.5 

10.6 

0.0 

0.5 

l.S 

3.7 

11.5 

0.0 

0.0 

0.0 

0.7 

IZ.S 

0.0 

0.0 

0.0 

0.0 

13.5 

0.0 

0.0 

0.0 

0.0 

14. S 

0.0 

0.0 

0.0 

0.0 

TRAP  HEIGHT  =  14. ZS  in 


Height    of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

for 

for 

for 

Tap    above 

Samples 

Samples 

Samples 

Samples 

Distributor 

Z16-ZZS 

ZZ5-Z3S 

Z36-Z4S 

Z46-Z55 

( inches ) 

(cm   HZO) 

(cm   HZO) 

(cm   HZO) 

(cm   HZO) 

l.S 

5Z.Z 

53.4 

56.4 

58.0 

Z.S 

48.7 

49.7 

51.  6 

5Z.7 

3.S 

40.9 

4Z.3 

43.9 

45.3 

4.S 

34.7 

3G.1 

38. Z 

40.  Z 

5.S 

Z7.3 

Z3.S 

31. Z 

34.1 

S.S 

Z0.6 

Z3.0 

ZS.7 

11 .1 

7.S 

14.9 

16.9 

18.9 

Z1.3 

8.S 

8.6 

10.7 

13.4 

14.4 

9.6 

Z.5 

3.7 

6.5 

8.9 

10. S 

0.3 

0.3 

l.Z 

3.1 

11.6 

0.0 

0.0 

0.3 

0.0 

IZ.S 

0.3 

0.0 

0.3 

3.3 

13. S 

0.0 

0.0 

3.3 

3.3 

14.6 

0.0 

0.0 

3.3 

3.0 
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Weight    of 

S, 

ample 

( grans 

) 

3.34 

5.3G 

4.S4 

5.71 

4.17 

4.33 

4.55 

4.33 

•8.44 

AVS           =   4.48 

grans 

3.30 

STN   DV    =   13. 72 

grains 

Trap   Height  above  Distributor:  en  (in)   =  36.20  (14.25) 

Bed   Height  above  Distributor:  en  (in)    =  28.58  (11.25) 

Trap  Height  above  bed  surface:  en  (in)    =7.62  (3.00) 

Static  Bed  Height:  en  (in)  =  22.56  (3.88) 

Pi:  en  Hg   (in  Hg)  =8.2  (3.23) 

dP:  en  uater   (in  uater)  =  6.38  (2.75) 

Tenp:   C  (F)  =  16.5  (62.0) 


SAMP  * 


246 
247 
248 
243 
250 
251 
252 
253 
254 
255 
*  value  not  used  in  conputing  average. 


Trap   Height  above  Distributor:  en  (in) 

Bed   Height  above  Distributor:  en  (in) 

Trap  Height  above  bed  surface:  en  (in) 

Static  Bed  Height:  en  (in) 

Pi:  en  Hg   (in  Hg ) 

dP:  en  uater   (in  water) 

Tenp:   C  (F) 


SAMP  8  Weight  of  Sanple 
(grans ) 

236  1.55 

237  1.88 

238  1.32 
233  1.80 

240  1.63 

241  2.35 

242  1.84 
2*3  1-80 

244  1.87    AVG     =1.88   grans 

245  2.38    STN  DV  =  0.22   grans 


198 


36.20 

(14.25) 

27.30 

(10.75) 

8.90 

(3.50) 

22.56 

(8.88) 

7.3 

(  2 . 87 ) 

4.768 

(  1 . 877  ) 

16.5 

(62.0) 

Trap   Height  above  Distributor:  cm  (in) 

Bed   Height  above  Distributor:  en  (in) 

Trap  Height  above  bed  surface:  en  (in) 

Static  Bed  Height:  en  (in) 

PI:  en  Hg   (in  Hg) 

dP:  en  uater   (in  uater) 

Tenp:   C  (F) 


SAMP  ft  Weighs  of  Sanple 
( grans ) 

22S  0.93 

227  0.68 

228  0.61 

229  0.5S 

230  0.81 

231  0.52 

232  0.38 

233  0.57 


= 

3G.20 

(14.25) 

s 

26.  S7 

(10.50) 

= 

9.53 

(3.75) 

s 

22.56 

(8.88) 

s 

6.6 

(2.60) 

= 

3.218 

(1.267) 

= 

15.5 

(  52 . 0 ) 

234  0.75   AVG    =  0.52   grans 

235  0.44   STN  DV  =  0.17  grans 


36.20 

(14.25) 

25.04 

(10.25) 

10.16 

(4.00) 

22.56 

(8.88) 

6.3 

(  2 . 48 ) 

2.624 

(1.033) 

16.5 

(62.0) 

Trap   Height  above  Distributor:  en  (in) 

Bed   Height  above  Distributor:  en  (in) 

Trap  Height  above  bed  surface:  en  (in) 

Static  Bed  Height:  Cn  (in) 

Pi:  cn  Hg   (in  Hg) 

dP:  cn  uater   (in  uater) 

Tenp:   C  (F) 


SAMP  *  Ueight  of  Sanple 
( grans ) 

216  0.32 

217  0.35 

218  0.36 

219  0.30 

220  0.32 

221  0-35 

222  0.24 

223  0.30 

Z24  0.34     AVG     =  0.32   grans 

Z25  0-30     STN  DV  =  0.04   grans 


i  qq 


Trap   Height  above  Distributor:  cm  <in)   =  3Z . 39  (12.75) 

Bed   Height  above  Distributor:  en  (in)    =  28.58  (11.Z5) 

Trap  Height  above  bed  surface-'  en  (in)    =3.81  (1.50) 

Static  Bed  Height:  cm  (in)  =  Z2.5S  (8.88) 

Pi:  en  Hg   (in  Hg)  =  8.1  (3.13) 

dP:  en  uater   (in  uater)  =  6.385  (2.75) 

Tenp:   C  (F)  =17  (53.0) 


samp  & 

Ueight  of  Sanple 

( grans ) 

285 

5. 55 

287 

7.51 

288 

7.46 

289 

8.05 

290 

S.9Z 

291 

7.38 

232 

7.23 

293 

7. 55 

234 

7.04    AVG 

295 

5.38    STN 

=  7.22   grans 
STN  DV  =  0.51   grans 


Trap   Height  above  Distributor 
Bed   Height  above  Distributor: 
Trap  Height  above  bed  surface: 
Static  Bed  Height:  en  (in) 
Pi:  en  Hg   (in  Hg) 
dP:  en  uater   (in  uater) 
Tenp:   C  (F) 


SAMP  tt  Ueight  of  Sample 
( grans ) 

276  2.44 

277  3 . 77 

278  3.53 
273  4.04 
Z80  4.55 

231  3.81 

232  4.55 
293  5. 50 


en  (in) 

=  32.39 

(12.75) 

en  (in) 

=  27.30 

(10.75) 

en  (in) 

=  5.09 

(2.00) 

=  22.55 

(8.88) 

=  7.3 

(  2 . 87 ) 

=  4.836 

(1.304) 

=  17 

(53.0) 

234  5.75   AVG    =  4.35  grams 

295  5.43    STN  DV  =  1.05  grams 
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Trap   Height  above  Distributor: 
8ed   Height  above  Distributor: 
Trap  Height  above  bed  surface: 
Static  Bed  Height:  en  (in) 
Pi:  cm  Hg   (in  Hg) 
dP:  cm  uater   (in  water) 
Temp:   C  (F) 


cm   (in) 

= 

3Z.39 

(12.75) 

cm   (in) 

= 

2B.G7 

(10.5) 

cm   (in) 

= 

5.72 

(2.25) 

= 

Z2.5G 

(8.88) 

= 

6.5 

(2.50) 

= 

3.259 

(1.287) 

= 

15.5 

(  52 . 0 ) 

MP   * 

Weight   of    Si 
(grams ) 

ample 

ZSG 

0.74 

257 

1.80 

2  58 

1.17 

253 

1.17 

270 

1.87 

271 

1.15 

272 

1.58 

273 

1.03 

274 

1.00 

AVG 

275 

1.55 

STN 

=  1.31   grams 
=  0.37  grams 


Trap   Height  above  Distributor:  cm  (in)   =  32.39  (12.75) 

Bed   Height  above  Distributor:  cm  (in)    =  26.04  (10.25) 

Trap  Height  above  bed  surface:  cm  (in)    =  6.35  (2.50) 

Static  Bed  Height:  cm  (in)  =  22.56  (8.88) 

Pi:  cm  Hg   ( in  Hg)  =6.3  (2.48) 

dP"-  cm  uater   (in  uater)  =  2.604  (1.025) 

Temp:   C  (F)  =  16.5  (62.0) 


SAMP  «  Weight  of  Sample 
( grams ) 

ZS5  0.50 

257  0.76 

Z58  0.67 

Z59  0.51 

260  0.54 

251  0.53 

Z62  0.47 

263  0.65 

264  0.61     AVG     =  0.63   grams 

265  0.87     STN  DV  =  0.13   grams 
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APPENDIX  J 


This  Appendix  contains  the  output  fron  the  inage  analyzer 
The  output  is  arranged  in  the  following  order. 


1.  Bed  distribution 

Z.  4  en  freeboard  height 

3.  8  en  freeboard  height 

4.  12  en  freeboard  height 

5.  18  en  freeboard  height 
B.  ZZ  en  freeboard  height 
7.  31  en  freeboard  height 


Each  of  the  distributions  at  a  given  bed  height  applies  only 
to  the  Uo/Unf  =  3.81  condition.   For  each  of  the  seven  inage 
analyzer  outputs  listed,  the  follouing  three  fornats  are  used. 

1.       Histogran  of  absolute  particle  frequency 
vs  particle  dianeter. 

Z.       Cunulative  percentage  plot  of  particle 
size  distribution. 

3.       Table  listing  of  the  above  data. 


Z0Z 


The  dashed  line  on  both  the  bar  graphs  and  on  the  cumulative 
percentage  plots  represent  the  Gaussian  distributions  uhich  fit 
the  given  set  of  data.   This  Gaussian  distribution  should  be  used 
only  as  a  rough  estimate  of  the  data  because  the  analysis  included 
"particles"  less  than  70  microns.   By  vieuing  a  blank  slide  with 
only  the  tape  applied,  these  "particles"  uere  confirmed  to  be 
bubbles  and  dirt  entrapped  in  the  adhesive  on  the  tape  used  to 
hold  the  sample  particles  in  place. 
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APPENDIX   K 


This  Appendix  contains  the  photographs  of  the  oscilloscope 
traces  obtained  uhile  sampling.   Each  of  the  pictures  is  labeled 
uith  the  sample  number  for  uhich  it  represents. 
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Oscilloscope  Trace  for  Sample  Number  56 


Oscilloscope  Trace  for  Sample  Number  57 
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Oscilloscope  Trace  for  Sample  Number  58 


Oscilloscope  Trace  for  Sample  Number  59 
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Oscilloscope  Trace  for  Sample  Number  B0 


Oscilloscope  Trace  for  Sample  Number  G2 
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Oscilloscope  Trace  for  Sample  Number  S3 


Oscilloscope  Trace  for  Sanple  Number  B4 


Oscilloscope  Trace  for  Sample  Number  SG 


Oscilloscope  Trace  for  Sample  Number  S7 
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Oscilloscope  Trace  for  Sample  Number  S3 


Oscilloscope  Trace  for  Sample  Number  70 


Oscilloscope  irace  for  Sample  Number  77 


"gygrV?-: 

:-1 ' 

=? 

■- 

?: 

>                            * 

M 

? 

. 

-. 

g 

. 

C 

S 

i 

- 

«.  - 

-i 

T- 

=* 

4 

■^ 

a 

i 

" 

1          : 

& , 


>a^ 


-  i 


■■-'•"•53 


MB  iir  — 


y 


Oscilloscope  Trace  for  Sample  Number  79 
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Oscilloscope  Trace  for  Sample  Number  83 


Oscilloscope  Trace  for  Sample  Number  84 
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Oscilloscope    Trace    for    Sample    Number    8S 


Jscilloscope    Trace    for    Sample    Number    88 
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Oscilloscope  Trace  for  Sample  Number  89 


Oscilloscope  Trace  for  Sample  Number  90 
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Oscilloscope  Trace  for  Sample  Number  92 


Oscilloscope  Trace  for  Sample  Number  9  3 


Oscilloscope  Trace  for  Sample  Number  34 


Oscilloscope  Trace  for  Sample  Number  95 
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Oscilloscope  Trace  for  Sample  Number  96 


Oscilloscope  Trace  for  Sample  Number  97 
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Oscilloscope  Trace  for  Sample  Number  98 


Oscilloscope  Trace  for  Samole  Number  99 
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Oscilloscope  Trace  for  Sample  Number  100 


Oscilloscope  Trace  for  Samole  Number  101 
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Oscilloscope  Trace  for  Sample  Nunber  102 
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APPENDIX   l_ 


This  Appendix  contains  a  complete  listing  of  the  particle 
trajectory  computer  model  used  in  this  research.   In  addition,  an 
in  depth  program  description  and  flow  diagram  is  included. 
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Model  Computer  Program 

This  section  describes  the  computer  program  listed  in 
Appendix  L.   To  assist  in  understanding  the  program  logic,  the 
program  itself  contains  numerous  comment  statements.   The  program 
consists  of  a  main  program  and  seven  (7)  subroutines.   The  main 
program  controls  the  input  of  parameters,  trajectory  calculations, 
and  output  selections.   The  subroutines  control  the  actual  data 
output  in  either  table  form  or  graphics.   The  following  is  an  in- 
depth  description  of  the  program: 


Lines  Description 


50-150        This  section  is  used  to  explicitly  define 

the  major  variables  used  uithin  the  program. 
The  large  arrays  are   defined  in  common  block 
form  to  save  memory.   Constants  used  in  the 
program  are  also  defined. 

133-S10       This  section  is  used  to  input  variable  data 
to  the  program.   Three  options  exist  for  the 
bed  particle  distribution  input: 

1)  The  default  condition  sets  the  quantity 
of  each  particle  to  unity.  This  option  is 
generally  used  uhen  a  height  determination 
is  required  or  being  sought  after. 

Z)  The  particle  number  option  allows  the 
entry  of  bed  distribution  by  the  number  of 
particles  present  in  each  diameter  range. 
This  is  used  uhen  aata  from  the  image 
analyzer  uas  being  used. 

3)  The  last  option  is  used  uhen  the  particle 
size  distribution  is  determined  from  a  sieve 
analysis  and  the  data  is  measured  in  grams 
mass.   The  program  will  then  determine  the 
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particle  number  density  based  on  the 
assumption  that  the  particls  are  spherical 


550-590       These  lines  convert  the  velocity  input 
values  from  Ft/s  to  m/s. 

700-710       Determine  the  slope  of  the  triangular  jet 
using  the  input  amplitude  and  duration. 

720-800       Initialize  the  data  arrays  to  zero  (0). 

810-1550      This  section  calculates  the  trajectory  for 

each  diameter  particle.   This  UHILE  condition 
contains  the  following  8  subsections. 

820-840     Initialize  the  height  of  the  particle  to 

pass  the  first  UHILE  statement.   Calculate 
the  particle  diameter  to  be  used. 

850-1340    Calculate  the  particles  trajectory 

parameters  while  the  particle  is  above  the 
bed  surface.   This  DO  loop  contains  the 
following  5  subsections. 

860-940  If  the  elapsed  time  since  the  particle  left 
the  bed  is  less  than  the  jet  duration  time, 
the  add  the  jet  velocity  to  Uo. 

950-990     Calculates  the  relative  velocity  of  the  air 
uith  respect  to  the  particle.   Determine  the 
sign  of  the  drag  force. 

1010-1060    Calculate  Reynold's  number. 

1100-1170   Calculate  the  drag  on  the  particle  and  then 
determine  the  particles  neu  velocity  and 
position. 

1180-1Z20    If  the  particles  velocity  is  positive, 

record  the  neu  max  height  and  elapsed  time. 

1230-1320.   Add  1  count  to  the  probability  array  in  the 
storage  position  representing  the  particles 
height/2.   By  using  the  DIV  statement,  a 
height  window  ( dH )  of  2  cm  is  created. 

1350-1370      Change  ell  heights  from  Ft/s  to  m/s 

1380-1430     Multiply  probability  distribution  by  bed 
distribution  weighting  factor. 

1440-1500      Determine  the  maximum  value  and  normalize 

Z37 


1510-1550 

1570-1560 

1700-1820 
1830-2210 

2220 
2250-2520 

2550-32  30 

2540-3120 
3130-3240 
3250-3830 

3840-3920 

3980-4130 

4180-4280 

4330-4460 

4510-4640 

4630-4770 


the   data   such   that    the   maximum   is    100. 

Calculate   the   particle   entrainment    by 
summing    the    total    volume   of    each   particle 
distribution   at    a  given   height. 

Normalize    the   entrainment    data   to    100    and 
calculate   the   In   value   if    not    equal    to   zero. 

Display   output    menu. 

Controls    the   selected   output. 

End   of   MAIN   program 

Sub  Display_Data  is  used  to  list  input 
parameters  used  and  the  maximum  height 
attained   by   each   diameter   particle. 

Sub   Oisplay_graphics   controls    the   processing 
and   output    of   graphic   information.      This 
subprogram   contains   the   follouing  4   sections. 

Performs    scaling   for   determining  graphic 
dimensions    and   limits. 

Creates   graphic   display   and   labels    the   X   and 
Y   axis. 

Controls    the   plotting   of    selected  graphic 
ou  t  pu  t . 

Prints   hard   copy. 

Function   routine   for   determining   scaling 
factor    for    graphic   display. 

Finds    the    limit    of    data   uithin   an   INTEGER 
array   to    limit    the   X   axis   on   a   plot    to   the 
range    in   which    the    Y    values    are    non-zero. 

Determines    the    maximum   value    stored    in   an 
INTEGER   array. 

Same    as    4330-4460    except    for    REAL    array. 

Same    as    4130-4230    except    for    RE3L    array. 

END   of   PROGRAM 
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13 
20 
53 
40 

5a 

53 

70 

30 

ec_?art 

90 

100 

us 

120 
130 

140 

150 
150 
170 
1S0 
190 
230 
Z10 

22a 

230 

240 

250 


.MAIN  Program. 


0PTI0N  3ASE  1 

INTESER  Answer  ,  Soeci  fy  .Point ,  Uiew,3ed_min,  I 

RE3L  Uo.Uoo.Uj .OeI_t,?art_heignt,SIooe, Jet_t, it 

CQfl  INTE5ER  Distribution<53,30),Distrib_density<53,33),REAL  Height< 2,53  ),3 


cle<50  ),S.Tt.-ai.nment<2,30  > 
DIM  Her.eSC2],21sar5C2] 
C:eer5-CH?.5<255)&CHRS(75) 
HomeS-CHRSC 255  )&CHRS( 34 ) 
Viscosi  ty_km*l . 4855-5 
Oensity_air«l .201 
Qen«ity_part»7.35E+3 
Srav:ty»9. 30555 

Incut  variaoles 


key 


prints; 

INPUT  " 

INPUT 

INPUT 

INPUT 

INPUT 


IS  1 
Sin  t  er 
5nter 
Enter 
r 
inter 


CLEAR  SCR 
HOME  key 
air,  m*«i/s 
k_/n*n*n 

kg/m*m*m 
m/s*s 


!  Output  to  CRT 


nean  air  velocity  Ua  (Ft/s)--  *,Uo 
ampiituce  of  jet  velocity  (Ft/s):  *iUj 
initial  particle  velocity  Upo  (Ft/s):  ".Uool 
duration  of  jet  (s):  ",Jet_i 
tine  increment  for  iteration  (3):  ",Del_t 


INPUT  "Do  you  desire  to  input  the  oed  particle  distribution  (1  yes, 3  no):" 


.Answer 


253 
270 
230 
220 
300 
510 
523 
530 
540 
353 
353 
370 
330 
3S0 
100 

410 
420 

430 

440 
450 

463 


IF  Answer«3~THEN  = 

FOR  1-1  TO  53 

3ed_particie(  I )-l. 

NEXT  I 
ELSE 

IF  Answer* 1  THEN 


-  -      .  .-       - 


Bed  dist  not  wanted 
Set  bed  dist  to  1 


!  Sed  dist  wanted 


INPUT  "Input  will  be;  weight  in  grans  (  1  )  or  number  <2):  ".Soecify 
IF  Soecify-1  OR  Soecify-2  THEN 

PRINT  "Enter  data  for  each  diameter"         !  Enter  gram  or  t 
Min_bed»l 

FOR  1-1  TO  53  !  Each  diameter 

PRINT  "D-  '; 1*10+70;"  urn" 
INPUT  3ed_particie(  I ) 

IF  3ed_particle( I XBed_particle(Min_bed)  THEN  Hin_bea-I 
NEXT  I 

IF  Soecify-1  THEN  !  If  weight  entry 

FOR  1-1  TO  S3  !  Calculate  5 

8ed_part icle( I )«5*3ed_part icle( I  )/(PI*< 1*10+70 >"3*0ensity_part ) 
IF  3ed_aart icie< I K5ea_part icie( Min_bea )  THEN  Min_bed=I 
NEXT  I 
END  IF 


490 

=30 
510 
523 
530 

340 

333 
353 


bulU  i=d 

ENO  IF 
END  IF 
Max_bea= 1 
FOR  !  =  '.  TG  33  '  Find  largest 

IF  Bed_part  icie1:  I  >:  Bea_par-.  :c  le(Max_bed  )  THEN  Haxjje-'"! 
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=32 

ss: 

333 
312 
5Z3 
333 
5-3 
3=3 
553 

T  *™  ■ 
-        <- 

533 
5=3 
"33 

"13 
723 
"33 

~-z 

"3 

7=3 
770 
780 

730 
333 

3;a 

IZZ 
333 
3-i0 
3=3 

353 

171 

383 

390 

=30 

=  10 

523 

=33 

3-13 

==3 

=53 

"Z 

=33 

=  30 

1030 

1013 

1SZ0 

1333 

1343 

1353 

1353 

.z~z 

1330 
1330 
1130 

1113 

1121 

.133 
11-3 
1158 


next  : 

?1ax_vaiue"130/3ec_3art  :c  is'.  "a.v_oed  )  '    detsr«iine    scaler 

FOR    -  =  -    "C    53  !    Scale    aeta    ne>.=  i30 

5ec_particiei  '.  ;=5ec_par-.  i ;  ie(  I  ^»nax_vaiue 
NEXT    I 

Conaute    aata 

Ul-Uo 
Up  i«'„'oal 
Uo-Jo-.33-iS 
Uool»'Jooi*  .  23-i3 
'J;l«'Jj*.30~8 
IF    Jet_:«0    THEN   ie*_t«i 
5  icae«2 .  *'J  j  1/  Jet_t 
FOR    :»1    70    =3 
Hei;r.t(  1,1  )-0. 
HeigntiZ, I )-3. 
FOR    J-l    TO   33 

QisiributionCI.J  I»0 
£n;riinnen:( 1 , J  !»0. 
Er.:riinnenx(Z,  J  )»0. 
NEXT   J 
NEXT    I 
FOR    I-i    TO   =3 

=  ar-._iei3nt*.  300001 
Uoo»Uool 

Oiane'.srj:ar:"(  1+7  )♦!  .0E-= 
'JHILE   Part_height>0 
it*Iet_t/Z. 

I .  -;_iou«4j:5(i:  (  3  .  I  >»QeI_t 
IF    Tine_now<Jet_t    THEN 

IF   Height(2,I)«0ei_t<Jt   THEN 

Uo»Uo-5iooe*Heignt ( 2 , I  >*Oei_t 
ELSE 

Uc«Uo+Uji-5iooe*Heignt(2,I )*Oei_t 
END  IF 
ENO  IF 

Reiat  i  ve^el'Uo-Uoo 
3 ign«i 
IF  Reiat ive_vei<3  THEN 

- I5n=-1 
ENO  IF 

Calculate  Re 
Reynoias_no»Reiat i  ve_ve  1*0 1  an-  ter_p  art  /Viscos  i  ty_k  in 
Reyncia5_no-'aiSS;  Reynoias_no  ) 
1=  Reynoids_no»0 .  "THEN 
-i:ei»rat ion«-9ravity 
SOTO  1153 
END  IF 

Calculate  dr 

!  Uses  aata  correlation  for  drag  coefficient  good  for  Re<lE5 
i 

I-2g_ccef  f  =  2i/Seynoias_no+5/(  l-^-SQRC  Reynoids_nc  ))-*-.  4 

3.-25=3  ign»Qrag_coe  f  f  *uer.s  i  ty_air  *Reiat  i  ve_vei  *Reiat  i ve_ve  i*3I*Qianeter 
►Ciafietar_part/8 

1  Calculate  acceleration 

-;cele_ation-3rac/(Censi  iy_cart  *=l  ♦2iarete'_^ar*.  "3/5  !— jrav:  ty 
!  Ca-cuiate  ael  velocity 

veioci  ty-'Jco--iCcelerat  ion»Qei_t 


!  Change  to  n/s 


!  Triangie  jet  pulse 

!  Zero  heigni  array 
■  Zero  nunoer  array 

!  Zero  dist  array 
1  Zero  entrain  array 


For  each  particle 
5et  for  UHILE  state. 
Initial  U  constant 
Particle  Die  un 


Tine  of  part  flignt 
part  in  jet? 
1st  half  of  jet 


!  2nd  half  of  jet 


Rei  v  seen  by  part 
Set  sign  positive 
Is  U  negative? 
Set  sign  negative 


ag 


240 


150  !                                              Calculate  new  position 

170  P2rt_heignt»p3rt_heignt+< Veiocity+Upo )»0el_t/2 

180  IF  Upo>0  "HEN                                !  Part  still  rising? 

120  HeigntC 1 , 1  )-Part_heignt                     !  Save  position 

200  Heignt<2, 1  )»Heignt( 2, 1 )+i                   !  Inc  t  for  max  rise 

210  ENO  IF 

220  Uoo»velocity                                 !  Set  new  U  for  next  inc 

230  ?oint*<  100»Part_.neignt  )  QIV  2                !  2  en  wide  storage  bins 

2*0  IF  Paint>72  THEN                            !  Set  default  for  fatal 

2S3  Pomt-?S 

250  ENO  IF 

2"0  IF   ?omt<«0  THEN                            !  Set  default  for  fatal 

230  ?omt«0 

2=0  ENO  IF 

330  PRINT  I,Point                               !  Indicate  coma  working 

313  !                                    Save  *  tines  part  in  height  bin 

323  Oistribution< I,?oint  +  l >«Qistribution( I ,Point  +  l  )+l 

330  END  WHILE 

3-3  NEXT  I 

3=0  FOR  1-1  TC  E3 

353  HeigntC i , I  )*Heignt( 1, 1  )*180                    !  Change  to  en 

370  NEXT  I 

384  FOR  1-1  TO  E0 

330  Ualue«3ec_?art icie(  I )            !  Part  sirs  weight  factor  from  bed  dist 

400  FOR  J»l  TO  30 

£10  Oistrib_aensity( I , J )*0istr ibut ion< 1 , J  WJalue   !  weight  dist  values 

£23  NEXT  J 

•130  NEXT  I 

■i-10  h"ax_dist=FNttax_int  (Oistr  ib_density(  *), 53,30  >       !  Find  nax  value 

453  Factor-130./Max_dist                              !'  Scale  for  L00  nax 

450  FOR  1-1  TO  33 

470  FOR  J-i  TO  30 

-80  Qistrib_densi ty< I , J )»Oistr ib_aensity( I , J )»Factor  I  Scale  values 

420  NEXT  J 

=33  NEXT  I 

=  10  FOR  1*1  TO  =3     !  Mass  density/unit  area   at  height  above  bed 

E20  Volune-PW  (  1-7 )/ 1000  )"3/S                       !  Volune  Cu  en 

E30  FOR  J-l  TO  30 

5-0  Entrainnent (  1 ,  J  >»Entrainnent  (  1 ,  J  )+yolune"0is  tr  ib_densi  ty(  I ,  J  ) 

333  NEXT  J 

553  NEXT  I 

E73  f1ax_entra  in»l 

533  FOR  1-1  TO  30                                   !  Find  nax i nun 

530  IF  Entrainnent  (  1 . 1  )>E.ntrainnent  I  1  ,Max_entram  )  THEN  Max_»ntram=I 

S30  NEXT  I 

510  Factor»i00/E.ntramnent  (  L  ,Max_entr3in  )               '  Normalize  to  100 

520  FOR  1-1  TO  30 

533  Entrainnent ( 1 , 1  )=Entr ainnent ( 1 , I  )*F actor 

5-iB  IF  Entrainnent ( 1,1  )<-0.  THEN  1S53 

1333  Entrainnent(2,I  )='_0G(  Entrainnent  (1,1)) 

1333  NEXT  : 

isTa  ' 

.3:3  '  Output  Control  Section 

1333  ! 

1704  =RINTEF  IS  1                                     '  Output  to  CRT 

171«  PRINT  USING  "9,3/" 

1720  PRINT  "  1  )      Oisolay  heignt  vs  aianeter  aata" 

1733  PRINT  "2)     Oisolay  neignt  vs  diameter  graon" 

L740  PRINT  'Z'             Oisolay  density  vs  neignt  as  function  of  dia  graon" 

1753  PRINT  '4      Oisolay  density  vs  diane'.ar  as  a  function  of  heignt  grapn" 

2^1 


I7S0  PRINT  "S)     Oisolay  sane  as  3  but  with  bed  density" 

1770  PRINT  "3)      Display  sane  as  4  but  with  bed  density" 

1780  PRINT  "7)     Oisolay  density  vs  diameter  of  bed  mass* 

1730  PRINT  "3)      Oisolay  entramment  density  aoove  bed" 

1800  PRINT  "3)     Oisolay  Ln  entramment  density  above  bed" 

1810  PRINT  "10)    EXIT  PROGRAM" 

1820  INPUT  "Enter  numoer  of  aesired  aispiay;  " .Answer 

1330  SELECT  Ansuer 

13^0  CASE  -1 

13E3  CALL  0isolay_data(Ul,Ucl,Uj,Jet_t,Dei_t,Hone5, Clears) 

1350  CASE  -2 

1S70  CALL  0ispiay_3raon( Clears, HoneS, Ansuer .View) 

1380  CASE  «3 

1333  INPUT  "Enter  particle  s;:s  to  be  viewed  (30-370  un)(0  for  all):  ".View 

1200  IF  v"ieui«0  THEN  1220 

1310  IF  Uieu<80  OR  yieu>E70  THEN  1330 

1320  View-CUieu  OIU  10 )-7 

1330  CALL  0isolay_3raoh( Clears, HoneS, Ansuer ,Uiau ) 

13A0  CASE  -4 

1SE3  INPUT  "Enter  desires  heignt  above  bed  surface  (0-1E3  en):  ",Uiau 

12S3  IF  vieu<0  OR  ViewMSS  THEN  1353 

1370  Uieu«<View  OIU  Z )+l 

1380  CALL  0isolay_3raon(ClearS, HoneS, Ansuer, Uieu  ) 

1250  CASE  *5 

2000  INPUT  "Enter  particle  size  to  be  viewed  (30-370  um)(0  for  ail):  ".View 

2010  IF  v-ieu-0  THEN  2030 

2320_  IF  Vieu<80  OR  Vieu>570  THEN  2000 

2030  <Jieu-(Vieu  OIU  10  )-7 

2040  CALL  Qisplay_graoh(ClearS, HoneS, Ansuer , View ) 

2053  CASE  -5  ~ ~" 

20S3  INPUT  "Enter  desirea  heignt  above  bed  surface  (0-153  cm)--  ",Vieu 

2070  IF  View<0  OR  Vieu>153  THEN  20S0 

2080  Uieu»(Uieu  OIU  2 )+l 

2030  CALL  0ispiay_3raomClearS,  HoneS,  Ansuer  ,  View  ) 

2100  CASE  -7 

2110  CALL  Qisolay_3raon(  Clears .HoneS , Answer ,Uiew ) 

2120  CASE  -8 

2130  CALL  Oisolay  graph ( Clears .HoneS. Ansuer .View  ) 

2140  CASE  "3 

21=3  CALL  Oisolay  araon( Clears , Howe S, Answer . Uieu  ) 

21S0  CASE  -10 

2170  STOP 

2180  CASE  ELSE 

2120  GOTO  1820 

2200  ENO  SELECT 

2213  GOTO  1700 

2220  END 

2230  ! 

2240  ! 

2253  ! 

22S3  :c:  jsed  to  list  incut  carareters  and  .la*  neigr.t  at  t  me  t  ser  diar.etar 

22  ~  3  Sun  3:j2iay_aata<  Ul .  Uol  ,'Jj  ,.'*:_:  .  2ei_: ,  HoneS,  CI  ear  3  ) 

2232  COM  INTEGER  Distribui ion( * ) . Oistr ;c_cansity( *  ■ , REAL  Height( *  )  ,3ed_ particle 

I  •  )  .Entramment ( *  ) 

2220  OUTPUT  2:Home5;                                   '  Hone  and  Clear  =creen 

2300  OUTPUT  2  .Clears. 

2310  1                                                    Print  output 

2323  PRINT  "                           cm/s" 

2330  °RirJT  "riean  3eo  Velocity-         ";U1»30.48 

Z"-0  PRINT  "Initiai  Dar,.;:;e  Velocity-  ";Uol*30.48 


2253  PRINT  "Peak  Jet  Velocity-         ";Uj»30.48 

23S0  PRINT  "Gas  Jet  Duration-          ";Jet_t 

2370  PRINT 

2380  PRINT  "   Diameter","    Heignf,"    Tine","     Diameter";  "    Height";" 

Tine" 

2320  PRINT  "      urn","     en   " , "   seconds',"        un","       en','           seconds 

2400  PRINT 

2H0  FOR  I- 1  TO  =0  STEP  2 

2-120  PRINT  USIN6  *2<SX,30,SX,30.D,5X,0.Q00,SX  )"  ;  C 1+7 >• 18, Height < 1,1  >,Height(Z 

,I)«Oel_t,(I*8)«l0,Heigni(l,I+l ), Height ( 2, 1+1 >»Gei_-. 

2430  NEXT  I 

24A0  PRINTER  IS  1 

2453  INPUT  "Print  hara  coay?  ( 1  >»  yes,  (0)"  no-     ", Answer 

24S3  I?  Ansuer»i  THEN 

2470  PRINTER  IS  701 

2480  SOTO  2320 

2480  ELSE 

2530  IF  Ansuer<>0  THEN  2453 

2510  ENO  IF 

2520  SUEENO 

2530 

■™c  * "» 

2553 

2553  !  Sub  used  to  control  graohics  outaut  of  data 

2573  SUB  0isoiay_3raon< Clears , Homes , INTEGER  Qata_set ,Uieu ) 

2530  COM  INTEGER  Distribution**  ),Qistrib_denaityt»  ), REAL  Height (*  ),8ed_particie 

( ♦  ) , En trainmen t  (  ♦ ) 

2530  REAL   Xmax,Ymax,Xtick,Ytick,Xmin,Ynm 

2S30  OUTPUT  2: HomeS ; 

2510  OUTPUT  2;ClearS; 

2523  SINIT                                         !  Initialize  graohics 

2530  GRAPHICS  ON 

25-0  SELECT  Oata_set                                  !  Scale  plot  routines 

2S53  CASE  -2 

2553  Xmax-533. 

2570  Ymax»FNKax_real(Heignt( ♦ >, 1 ,50 ) 

2S30  Xtick-53. 

2530  Ytic!c=»FNScale(Ynax  ) 

2730  CASE  -3 

2710  Ynax-'rNnax_Lnt(Oistnbution<  ♦  ),S3,30) 

2720  Xmax-PNOata.limiMDistributionC  ♦  ),53,30  ) 

2730  Xtick-FNScaiaUmax  ) 

2740  Yt:ck-FNScale<Ymax  ) 

2750  CASE   -4 

3753  Xmax-530. 

2770  Ynax-FNMax_Lnt(0i5tnbution<  ♦), 50,30  > 

2780  Xt:ck=53. 

2730  Ytick»FNScaie< Ymax  ) 

2300  CASE  »5 

2910  Xmax--N0ata_Linit(0istrib_'3ensity(  *  )  ,53,30  ) 

2323  Ynax»!20. 

2333  Xtick*FNScaie(Xnax  ) 

2340  Vtick«FNScale<Ynax) 

2953  CASE   »G 

23G3  Xnax=G30. 

2973  Y«ax-FNftax_int<Oisirib_density<  •  ),  53,30  ) 

2980  Xtick-53. 

2390  Ytick»FNScaie< Ymax ) 

2  390  CASE    -7 

24j 


2913  Xnax«500. 

2323  Ynax«100. 

2=30  Xtick-50. 

2940  Ytick-10 

2SS0  CASE  -e 

2550  Xnax»FNOata_linitr( Entrainment ( * ) ,  1 , 30  ) 

2370  Ynax-100. 

2580  Xtick«FNScaie(Xmax > 

2530  Ytick-10. 

3000  CASE  -3 

3013  Xnax--NCata_ii.*!itr(EntrainHerit  (♦),!,  33) 

20Z3  Ynax»=. 

3030  Ynm— 3. 

33*3  XticK»-NScale(Xnax  ) 

30S0  Ytick-.S 

30S0  Y_axia— 2.5 

3070  50T0  3110 

3080  5N0  SELECT 

3053  Y_axis»0. 

3100  Ymn— Ytick 

3110  Xnm— 2.*Xtick 

3i20  X_axis«0. 

3130  uInOC'J  Xnin,l.l*Xmax,Y!nin,i.l»Ynax 

3U0  AXE5  Xtick,Ytick,X_axis,Y_axis 

31E0  LORG  5 

3150  FOR  I-X_axis  TO  Xnax  STEP  Xtick»2 

3170  flOVE  I,Y_axi5 

3180  LABEL  I 

3130  NEXT  I 

3200  L0R6  3 

3210  FOR  I-Y_axis  TO  Ymax  STEP  Ytick 

3223  HOVE  X_axxs,I 

3230  LABEL  I 

3ZA0  next  I 

3250  SELECT  Oata_set 

32=0  CASE  "2 

3273  .-.OWE  80, Height (1,1) 

3Z33  FOR  1-2  TO  50 

3253  ORAU  ( 1+7  )•  10. Height  U,  I ) 

3300  NEXT  I 

3310  CASE  »3 

3323  3egm«Vieu 

3330  Finisn-'Jiew 

3340  IF  Vieu<3  THEN 

3353  3egin-l 

33S0  F:nish»53 

3370  END  IF 

3380  FOR  I»3egin  TO  Finish 

3333  MOVE  B,3istribution( I,i J 

3433  FOR  J-L  TO  Xnax/ 2 

3*1 3  ORAU  Z*J,Oistributicn( I, J) 

3-2  2  ':=:<-   J 


Set  graphics  scale 
Set  graoh  axis 
Label  X-axis 


!  Label  Y-axis 


!  Plot  data 


!  Plot  all  if<3 


34-48 

3*sa 

2453 
3473 
3480 
3453 
3530 


CASE  «* 

."!C<JE  B0,Distribution(  l.Viaw) 

FOR  >2  TO  53 

ORAU  ( 1+7)*  10, Distribution  I  ,'Jieu  ) 
NEXT  I 
CASE  =5 

3eq :  n=')  isu 


Z44 


3S10  Finish-View 

3523  IF  View<0  THEN                              !  Plot  all  if  <3 

3533  2egm*l 

354-3  Finisn-53 

35=3  END  IF 

3S53  FOR    I»9egin    TO    Finisn 

3573  ,10UE   3,0istrib_density<  1,1  ) 

3533  FOR   J-i    TO   Xmax/2 

35=3  J3RAU  2»J,J3i3trib_aensity<  I,  J  ) 

3533  NEXT  J 

3513  NEXT  I 

33Z3  CASE  -5 

3533  MOVE  30,Oistrib_density<  1,'Jiew) 

35-3  FOR  1-2  TO  53 

3553  DRAW  < 1+7 )♦ 13 ,Distr ib_densi ty<  I ,V Leu ) 

3553  NEXT  I 

3570  CASE  -7 

3533  nOUE  30.3ed_aarticle<  1  > 

3533  FOR  1-2  TO  53 

3733  ORAU  (I+7)*lQ,8sd_pariicls<I) 

3710  NEXT  I 

3723  CASE  -8 

3733  ttOUE  3,£niramment<  1,1  3 

374.3  FOR  I«2  TO  33 

3753  ORAU    [•2,£ntraiiwient(l,D 

37S3  NEX7   I 

3773  CASE   -9 

3733  HOVE   3,Eritrammeni(2,l 3 

3730  FOR    1-2   TO   30 

3303  IF  Entrainment<2,I 3-0.  THEN  38*3 

3313  ORAU  I*2,Entrainment(2,I 3 

3320  NEXT  I 

3830  ENO  SELECT 

384.0  INPUT  "Print  graoh  (13  yes,  (0)  no:  '.Answer 

3853  IF  nnsuer'i  THEN 

3353  DUMP  DEVICE  IS  701                             !  Duma  to  printer 

3873  DUMP  GRAPHICS 

3583  SCLEAR 

3833  ELSE 

3803  IF  Ansuer<>0  THEN  38*0 

3310  ENO  IF 

3320  3CLEAR 

3330  SU8ENO 

334.0  ! 

3553  ! 

3=50  i 

3370  i  Function  determines  graohics  axis  scaleing 

3330  OEF  FNScaleiYnax 3 


130  AND  fmex>20  THEN 


3330 

IF  Ymav;-  100  THEN 

4.000 

7ick=Z3. 

4010 

c!*«  c 

4020 

IF  /.-.ax'-lOC  ml 

-233 

TicK-ia. 

4.04.0 

—  -1Z. 

4.053 

IF  'nax->Z3 

4.353 

7  i  c  y.  —  Z 

4073 

ELSE 

4083 

" :  c  k  =  .  S 

4030 

EMO  IF 

4  130 

EMO  IF 

Z4S 


1110 

4120 
4130 
4140 

tisa 

^150 

4170 

A130 

•1120 
4200 

4210 

4223 

4230 

4240 

42S0 
4250 
4270 
4290 
4250 
4500 

4310 

4320 
4330 

4340 
4350 
4350 
4370 
4380 
4320 

14.00 
4.110 
4.120 
4430 
1110 
4150 
44S0 
4170 
4180 

449fl 

4500 
4510 

4520 
4530 

4540 

4550 

4550 
4570 
4530 

4  =  50 

—  0  XJ  V 

-D  1  Z 


10 

:0 

4550 
4670 
4530 
4S30 
4702 


5ND  IF 
RE7URN 
FNENO 


"ick 


1  Function  determines  max  extent  of  data  for  X-axis  limit 
DEF  FNQata_limit( IN7ESER  Qata_l( * ),Row_max ,Col_max ) 
Cai-1 
REPEAT 

Sum-0. 

FOR  1-1  70  Rou_max 

Sum-ium-rOata_l<  I. Col  ) 

NEXT    I 

Col-Col+1 
UNTIL  Sum-0  OR  Col-Col_max-l 
RETURN  Co W 
FNENO 


Function  determines  max  value  in  integer  array 
OEF  FNC!ax_int(  INTE5ER  0a:a_l(  ♦  )  ,Rou_max  ,Col_max  ) 
R_max»L 
C_m«x*l 

FOR  I«l  70  Rou_max 
FOR  J-i  70  Col_max 

IF  0ata_l(  I,J)>0ata_l(R_max,C_max)  7HEN 
R_max»I 
C_max-J 
END  IF 
NEX7  J 
NEX7  I 

Max-«0ata_l<  R_max , C_max ) 
RE7URN  flax 
FNENO 


Function  determines  maximum  value  in  real  array 
OEF  FNMax_r-al(Data_l< • ),IN7E5ER  Row_max ,Col_max ) 
R_max»l 
C_max"l 

FOR  1*1  70  Roui_max 
FOR  J-i  70  Col_max 

IF  Oata.K  I.J  )>0ata_l(R_max,C_max>  7HEN 
R_max-I 
C_mex»J 
END  IF 
NEXT  J 
NEXT  I 

.**.£>:=0a*a_l  <  ?._nax  ,  C_max  ) 
=E~URN  Rax 
=  UEN0 


1  Function  determines  ma/,  extent  of  data  for  X-axis  limit 
OEF  rNCata_l  imitr(Daia_l(  *  )  ,  Row,  Col_ma>.  ) 
Col  =  l 


Z4G 


4713 

IF  Data_l(Row,Col X 

4728 

Lim it-Co  1*2 

-1733 

RETURN  Limit 

1743 

ENO  IF 

4753 

Coi-Coi-*-! 

4753 

SOTO  4713 

4773 

FNENQ 

■8    OR   Col -Col    max    THEN 


:-i7 


APPENDIX   M 


This  Appendix  contains  the  calibration  data  for  the 
anemometer  probe.   The  calibration  uas  conducted  in  a  small  uind 
tunnel  using  a  pitot  tube  connected  to  a  micromanometer  capabl a  of 
measuring  pressures  to  uithin  0.801  ins.   of  water. 
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Fig.  M-l    Calibration  of  Anemoneter  Probe. 

Oscilloscope  Voltage  vs  Air  Velocity, 
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